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Abstract 

An overview of the random network model invented by Chalker and Coddington, 
and its generalizations, is provided. After a short introduction into the physics of the 
Integer Quantum Hall Effect, which historically has been the motivation for intro- 
ducing the network model, the percolation model for electrons in spatial dimension 2 
in a strong perpendicular magnetic field and a spatially correlated random potential 
is described. Based on this, the network model is established, using the concepts of 
percolating probability amplitude and tunneling. Its localization properties and its 
behavior at the critical point are discussed including a short survey on the statistics 
of energy levels and wave function amplitudes. Magneto-transport is reviewed with 
emphasis on some new results on conductance distributions. Generalizations are 
performed by establishing equivalent Hamiltonians. In particular, the significance 
of mappings to the Dirac model and the two dimensional Ising model are discussed. 
A description of renormalization group treatments is given. The classification of 
two dimensional random systems according to their symmetries is outlined. This 
provides access to the complete set of quantum phase transitions like the thermal 
Hall transition and the spin quantum Hall transition in two dimension. The su- 
persymmetric effective field theory for the critical properties of network models is 
formulated. The network model is extended to higher dimensions including remarks 
on the chiral metal phase at the surface of a multi-layer quantum Hall system. 
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1 Introduction 



At the end of the seventies of the past century, many thought that solid state 
physics had matured to such a degree that practically everything important 
had been discovered. Superconductivity seemed to be well understood with the 
transition temperatures ceasing to increase further. Semiconductor physics 
had developed almost into an engineering discipline. No significant further 
progress of the field of solid state physics was predicted for the foreseeable 
future. It was widely believed in the community that the development more 
or less had come to an end. In this situation, a completely new and by no one 
foreseen phenomenon was discovered in the magneto-transport properties of 
a commercial electronic device, the Silicon MOSFET. The Integer Quantum 
Hall Effect was found. 

This experimental discovery, together with several others that came roughly at 
the same time, opened a completely new area of solid state research and during 
the forthcoming years initiated a truly novel view on the field of condensed 
matter. This concerns the quantum mechanical properties of disordered and 
interacting electronic solid state systems on mesoscopic scales of which the 
Quantum Hall Effect is only one example, though a very prominent one. Until 
today, researchers in the field of the Quantum Hall Effect continue to produce 
new surprises, perhaps not on daily, but certainly on monthly time scales. In 
many cases, these concern only at first glance the smaller area of the quantum 
Hall phenomenon. Often, as in the case of the fractional Quantum Hall Effect, 
the discovered phenomena later turn out to be of much wider importance than 
foreseen at the time of their discovery. 

In this review article 2 we want to describe one example of such a development, 
namely the discovery and further development of a theoretical model which 
originally was designed to describe a special aspect of the integer Quantum 
Hall Effect, namely the localization-delocalization quantum phase transition 
in a Landau band in the limit of long-range randomness. Later, the model - 
although at first glance very specialized and restricted — has been shown to 
be able to account for many phenomena in a much wider class of systems, 
namely the disorder-induced quantum phase transitions in seemingly all of 
the different presently considered disordered systems. The model, invented by 
Chalker and Coddington in 1988, describes a random network of currents. 

2 In contrast to the short overview that has been published earlier [B. Kramer, S. 
Kettemann, T. Ohtsuki, Physica E 20, 172 (2003); cond-mat/0309115 , the present 
article is supposed to provide a much more self-contained — as much as it is possible 
- overview of the derivation of the network model and its quantum mechanical 
properties, amended by some recent results. Much emphasis is on generalizations 
of the model and many connections to other models and descriptions of quantum 
phase transitions, especially in two dimensions. 
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In the following sections of this introductory chapter, we first provide shortly 
some insight into, and understanding of the integer Quantum Hall Effect as it 
has been originally detected. Our considerations will be based on the so-called 
localization model. This is used for describing some fundamental physical as- 
pects of the Integer Quantum Hall Effect, the interplay between a high - 
so-called quantizing — magnetic field, and the disorder in the system due to 
impurities. This interplay can lead to a localization-delocalization transition 
in a two dimensional electronic system. To the best of our present knowledge, 
this appears to be a paradigm of a genuine quantum phase transition. Together 
with the gauge argument first proposed by R. B. Laughlin, this can explain not 
only the very existence of plateaus in the Hall conductivity as a function of the 
electron density, and the simultaneous vanishing of the magneto-conductivity, 
but also the quantization in integer units of e 2 /h. The Chalker-Coddington 
network model is now widely accepted as one possibility for describing the 
fundamental physics behind the quantum Hall phase transition. While the 
latter is nevertheless still waiting for a complete and quantitative theoretical 
description with predictive power, especially including the precision aspect, 
the network model of Chalker and Coddington seems to have acquired more 
fundamental importance also in other fields of solid state physics like super- 
conductivity and magnetism, as will be seen below. 

1.1 The Discovery of the Quantum Hall Effect 

The Quantum Hall Effect has been discovered by Klaus von Klitzing in 1980 
when working as a guest researcher at the High Magnetic Field Laboratory of 
the Max-Planck-Gesellschaft in Grenoble [1] (Fig. 1). He was investigating the 
electronic transport in a Silicon MOSFET subject to a high magnetic field 
of about 18 T flux density at temperatures near 1 K. To his great surprise, 
he found that the Hall resistance -Rr- — the ratio between the Hall voltage 
across the two dimensional electron inversion layer in the transistor, and the 
source-drain current, (J = 1/xA) - - shows extremely well-defined plateaus 
when changing the gate voltage U g . In the gate voltage regions of these Hall 
plateaus, he found the longitudinal magneto-resistance to be vanishingly small. 
Most strikingly, he was able to identify the values of the Hall resistances of 
the plateaus as integer fractions of h/e 2 

Rk = -- 2 (j = 1,2,3,...). (1) 

J e 

Here, h is the Planck constant and e the elementary charge. 

Klaus von Klitzing found that the relative uncertainty of the plateau values 
was much better than 1CT 5 in the very first experiments. Thus, as a surprise, 
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despite of the presence of strong disorder and electron interactions in the 
MOSFET, the effect appeared to be very promising for measuring the Som- 
merfeld constant a = (/i c /2)(e 2 //i) (c = 299 792 458 ms -1 speed of light in 
vacuum, /x = 47r • 1CT 7 NA~ 2 permeability of vacuum) independently from 
optics. Also, von Klitzing immediately realized the importance of his finding 
for metrological applications, i.e. the realization, dissemination and mainte- 
nance of the "Ohm", the unit of the electrical resistance. Eventually, this 
led to the re-definition of the "Ohm" in terms of the von Klitzing constant 
Rk = 25 812.80851] in 1990. Meanwhile, the experimental reproducibility of 
the quantized values of the Hall resistance has been improved to values better 
than 10~ 9 . This is of special importance in view of the metrological applica- 
tions [2]. Some excellent reviews about the Integer Quantum Hall Effect can 
be found in the literature [3,4,5,6]. 

The discovery of the von Klitzing effect stimulated intense experimental and 
theoretical research in many international laboratories. A particularly impor- 
tant discovery was made only two years later by the group of Daniel C. Tsui, 
Horst Stormer and Arthur C. Gossard at Bell Laboratories [7]. In this experi- 
ment, a GaAs/AlGaAs-heterostructure which contained the dopant impurities 
only far from the electron inversion layer was used. As a consequence, these 
samples had a very high electron mobility in the inversion layer at the interface 
of a few hundred thousand cm 2 /Vs. The researchers detected quantization of 
the Hall conductivity at fractional multiples of e 2 /h at temperatures of a few 
0.1K. 

By improving the sample fabrication technology, more and more such addi- 
tional fractional features were uncovered. Hitherto, electron mobilities of more 
than 10,000,000 cm 2 /Vs have been achieved and several dozens of the fraction- 
ally quantized features in the Hall resistance (and their counterparts in the 
magneto- resistance) are known now [8]. They have been associated with hi- 
erarchies of novel correlated electron states induced by the electron-electron 
interaction [9,10,11]. The existence of these states goes far beyond the tradi- 
tional phenomenological Fermi liquid picture for electrons in metallic systems. 
Thus, the Fractional Quantum Hall Effect has opened a new field in the physics 
of correlated electrons, but is not the subject of this review article. 



1.2 The Localization Model 

According to the theory of the classical Hall effect the Hall resistivity of inde- 
pendent particles with the charge e is a monotonic function of the magnetic 
flux density B and the number density n e of the charges. This can be easily 
seen by using the Drude friction model for diffusion of charges in external 
crossed electric and magnetic fields [12]. It yields for the components of the 
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Fig. 1. The Integer Quantum Hall Effect [1]. The Hall voltage Uu (left axis) shows 
plateaus in regions of the gate voltage U g where the source drain voltage U (right 
axis) is vanishingly small. The experiment was done on a Si-MOSFET of length 
400 /xm and width 50 nm at B = 18 T and at a temperature T = 1.5 K. 



conductivity tensor 



a , 



cr 



1 + Kr) 2 ' 



(2) 



Here, <To = e 2 n e r/m* is the Drude conductivity. The mean free time r of the 
particles due to scattering is assumed to contain all microscopic processes. 
The effective mass m* contains the effect of the lattice of atoms and interac- 
tions. The quantity ujb = eB/m* is the cyclotron frequency. The resistance 
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components are obtained by inverting the conductivity tensor 



a. 



^ xx " yx / \ 

Pxx — ~~^> j K~ ' P X V ~ ~~2 i 2~ ' v ' 



xx 1 yx xx 1 ^ yx 



For strong magnetic field (B — > oo) a xx oc 5 2 r 1 — > 0. Correspondingly 
since er ya; 7^ 0, the magneto-resistivity p^ oc cr^ — * 0. The Hall resistivity is 

Pxy := Ph = = • (4) 



As the von Klitzing experiment contradicts this result, it is clear that the 
quantum nature of the two dimensional electron system in the MOSFET sub- 
ject to the high magnetic field must be the reason for the quantization of the 
Hall resistivity at sufficiently low temperatures. 

For the Landau model for a single spinless electron in two dimensions in a 
perpendicular homogeneous magnetic field B = V x A, the vector potential 
given in the Landau gauge A = (0, Bx, 0), the Schrodinger equation is 

1 (h \ 2 

H ipnk{x,y) ■= — ; f -V + eAj ip nk (x,y) = E nk if) nk (x,y) . (5) 



It is solved by the Landau states 

(x, y\nk) := i[> nk (x, y) = e lky <p n (x - X k ) (6) 



where X k = —M\ and £b '■= \Jh/eB the magnetic length. The eigenvalues 
are the Landau levels 

E„k = huj B (n + ^ . (7) 



They are associated with the normalized eigenfunctions 
1 



M*) = 1 = exp (J[pj H n (J- ) (8) 



with the Hermite polynomials H n (n = 0, 1, 2, 3, . . .). The Landau states are 
degenerate with respect to the wave number k with a degree per unit area of 

eB B 

nB = T = * > (9) 
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the density of flux quanta, $o = h/e, in the system. This degeneracy can be 
easily understood by considering the maximum possible number of Landau 
wave functions (Eq. (6)) that can be contained in a system of the size L 2 [12], 



max 



2ir£ 2 B 



n B = L. 



(10) 



k 



L 



An important quantity is the filling factor which is defined by the ratio of the 
electron number density and the density of flux quanta, 



Integer filling factors v B = j correspond then to j completely filled Landau 
levels. At the corresponding electron densities, the Hall resistance is an integer 
fraction of h/e 2 (cf. Eq. (4)). 

However, this does not explain the existence of the wide plateaus in p xy {n e ) 
since upon increasing n e the Fermi level jumps between the Landau levels. In 
order to generate the observed wide plateaus, it is necessary to keep the Fermi 
level continuously varying between the Landau levels when changing n e , but 
without changing the resistivity values of the plateaus. 

A mechanism that can account for this pinning without affecting transport is 
localization of the wave functions due to disorder induced by the presence of 
the impurities in the system [13,14,15]. By introducing randomness into the 
Landau model of Eq. (5), 



the degeneracy of the Landau levels is removed. Here, V(r) is a randomly 
varying potential which is usually defined via its statistical properties. The 
Landau levels are broadened into Landau bands by the disorder, with localized 
eigenstates occurring in the band tails. These localized states can pin the Fermi 
level at the corresponding eigenenergies. A qualitative picture of the density 
of states of such a model Hamiltonian is shown in Fig. 2. 

Localized states correspond to random wave functions with envelopes that are 
exponentially decaying at large distances from some localization center [16,17]. 
This implies exponentially decaying correlation functions which indicates that 
the localized states cannot contribute to dc-transport at zero temperature. 
This means that the magneto-conductivity, and correspondingly the magneto- 
resistivity, vanish for filling factors corresponding to the energy regions of the 
localized states, while the Hall conductivity, and correspondingly the Hall 
resistivity, stay at a constant value. 



n e 

v B = — 

n B 



(11) 



H = H + V(r) 



(12) 
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On the other hand, non-localized — extended — states near the band cen- 
ters, which do not decay exponentially and are spread randomly across the 
entire system, can account for electron transport. For filling factors close to 
half integer numbers, which correspond to the energy regions near the cen- 
ters of the Landau bands, one expects peaks in the magneto-resistivity with 
widths that reflect the widths of the energy regions of the extended states. 
Simultaneously, the Hall resistivity is expected to change from one plateau 
value to the next. Often, these spectral regions are denoted as compressible 
since the Fermi energy hardly changes with increasing electron number while 
the localized states are associated with incompressible spectral regions where 
the Fermi level strongly changes when varying the electron number [18]. 

As long as the Fermi level stays within the region of localized states, the 
zero temperature conductivity components are not changed. Assuming that 
the extended states provide the correctly quantized values of the Hall con- 
ductivity, and simultaneously a vanishing magneto-conductivity, one obtains 
the experimentally observed behavior of the conductivity. However, that the 
Hall conductance remains correctly quantized exactly at integer multiples of 
e 2 /h in the regions of localization, needs closer consideration. Laughlin's gauge 
argument, which we will shortly discuss below, is supposed to provide an ex- 
planation. 

All of the analytical and numerical results obtained until recently are consis- 
tent with the picture that at the absolute zero of temperature the localization 
length diverges only at specific energies close to the centers of the Landau 
bands with a universal exponent v which is independent of the specific form 
of the randomness and the band index [19,20], 

^ E - E ^=w^kr (13) 

where the constant £ n depends on microscopic details of the randomness and 
on the Landau band index n. It will be one of the main tasks of the following 
sections to explain this critical behavior in some detail. 

Table 1 contains a representative selection of the exponents determined nu- 
merically. If the exponent were universal the quantum Hall effect could be 
considered as a paradigm of a genuine universal quantum phase transition. 

An important issue is how to detect this critical exponent in an experiment. 
As in an experiment the size of the two dimensional electron system is fi- 
nite, it is intuitively clear that the singularities of the localization length near 
the band centers cannot be resolved. Furthermore, at non-zero temperature, 
inelastic processes will lead to phase breaking scatterings between different 
eigenstates [16,17]. The latter can be described by a mean, temperature de- 
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Fig. 2. Schematic picture of the density of states (full line, arb. units) and the 
localization lengths (dashed line, arb. units) as a function of the energy (in units of 
the cyclotron energy) in the localization model. When introducing a static random 
potential into the Landau model, Landau levels are broadened into bands consisting 
of "incompressible" spectral regions that correspond to localized states in the band 
tails, and "compressible" regions of effectively delocalized states (shaded) near the 
band centers. In the latter regions, the localization length exceeds the system size 
L. At finite temperature, the system size is replaced by the temperature dependent 
phase coherence length L^. 

pendent phase breaking time r^(T) which can be assumed to increase with 
decreasing temperature according to 



oc 



1 



(14) 



with an exponent p of order 1. During the time interval r^, the electron can be 
considered as moving diffusively under the influence of the impurity scattering. 
This suggests to define a phase coherence length 

L^T) = y/Dr^T) (15) 



with a disorder-induced diffusion constant D that is related to the dc-conduc- 
tivity via the Einstein relation a = e 2 D(Ep)p(Ep) where p(Ep) is the density 
of states at the Fermi level. 

In order to connect the singular behavior of the localization length with the 
results of a transport experiment it is important to note that a localized state 
appears extended if its localization length does exceed the system size L or, 
at non-zero temperature, the above phase coherence length (Fig. 2). This 
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model 


method 


reference 


oo 


short-range impurities 


self-consistent perturbation 


[211 

L J 




Peierls tight binding 


transfer matrix scaling 


[221 

L J 


« 2.0 


short-range impurities 


recursive Green function 


[23,24] 


2.35(3) 


random Landau matrix 


recursive Green function 


[19,20] 

L * J 


2.3(1) 


random Landau matrix 


recursive Green function 


[251 

L J 


2.4(2) 


random Landau matrix 


recursive Green function 


[26] 

L J 


2.4(1) 


finite range impurities 


Chern number scaling 


[27] 

L J 


2.3 


spin orbit scattering 


Thouless number scaling 


[28] 

L J 


2 


double layer system 


Thouless number scaling 


[29] 


Rj 2 


random matrix model 


scaling of level statistics 


[30] 


2.5(5) 


Chalker-Coddington 


transfer matrix scaling 


[31] 


2.4(2) 


random saddle points 


transfer matrix scaling 


[32] 


2.5(5) 


Chalker-Coddington 
type 


real space renormalization 


[33] 

L J 


2.39(1) 


Chalker-Coddington 
type 


real space renormalization 


[34,35,36] 


2.5(4) 


super spin chain 


density matrix renormaliza- 
tion 


obtained 
from [37] 


2.33(3) 


counter-propagating 
chiral Fermions 


Monte Carlo 


[38] 



Table 1 

Critical exponents v of the quantum Hall phase transition in the lowest Landau 
band obtained by various theoretical methods (numbers in parentheses denote the 
uncertainty in the last digit of the exponent). 



implies that the widths of the peaks in the magneto-conductivity at the centers 
of the Landau bands are given by the condition £(AE) = min{L, L^(T)}. The 
resulting characteristic temperature behaviors [19,20] have been found to be 
consistent with experimental data [39,40,41,42,43,44]. 

One still needs an argument as to why the Hall conductance is exactly quan- 
tized at integer multiples of e 2 /h in spite of the presence of the impurity 
potential. Such an argument has been pioneered by Laughlin [45]. The idea 



13 



is to relate the current in an ideally metallic cylinder, subject to a homoge- 
neous magnetic field perpendicular to the cylinder's surface, to the change of 
the total electronic energy E(<&) when adiabatically changing a magnetic flux 
piercing the cylinder along its axis (Fig. 3), 

AE , . 



Introducing a gauge flux $ along the axis of the cylinder corresponds to a 
shift of the wave number in azimuthal direction by e^/hL. This can also be 
considered as a change in the boundary condition. 

A change of the flux by A$ = $o corresponds to a shift of exactly 2n/L 
and leads to a corresponding shift in Xk (cf. Eq. (6)) that is equivalent to 
transfering one electron per Landau band from one edge of the cylinder to the 
other. 3 When occupying j Landau bands, a transfer of j electrons is associated 
with one flux quantum. Due to the presence of a Hall voltage Un across the 
cylinder, the Landau levels are no longer degenerate, but on a straight line 
with a slope given by the Hall electric field Uu/L (Fig. 4). The energy change 
caused by the flux change is then AE = jel/u, and the Hall current 

/ = JT f/ H, (17) 



consistent with integer quantization of the Hall conductance. 

As presented, the argument seems to hold only for the ideal Landau Hamilto- 
nian applied to a system of finite width. However, it can also be used in the 
presence of disorder [22]. To understand this, one has to have in mind that 
on the one hand, by definition, localized states are insensitive to changes in 
the boundary conditions or, equivalently, flux changes. On the other hand, 
the flux sensitivity of the extended states is enhanced, in order to compensate 
for the localized states [3,4,5,47,48,49,50]. This compensation is assured to be 
exact such that the final result for the current is the same as in Eq. (17). 

Despite there exist several experimental features apparently related to electron 
interaction [51,52,53,54,55,56,57,58,59], it seems that the integer quantization 
of the Hall conductance can be understood within the localization model with- 
out taking into account the correlations between the electrons [5]. 

3 Strictly speaking, for a two dimensional system with periodic boundary conditions 
in one direction and fixed boundary conditions in the perpendicular direction, one 
finds a non-degenerate energy spectrum that differs from the Landau spectrum by 
the presence of edge states [12,46]. However, for the present qualitative argument, 
this is not important. 



14 



E A 



E2 




E1 





Fig. 3. Illustrating Laughlin's gauge argument: introducing a magnetic flux <3? along 
the axis of the cylinder corresponds to a shift of the wave number in azimuthal 
direction by e&/hL. A change of the gauge flux by A<1> = $0 corresponds to a shift 
of exactly 2ir/L. This leads to a corresponding shift in that is equivalent to 
transfering one electron per Landau band from one edge of the cylinder to the other 
that corresponds to a change in energy of ellu. 

On the other hand, for explaining the fractionally quantized Hall effect, inter- 
action and correlation effects are generally accepted to constitute the neces- 
sary ingredients for generating the excitation gaps between the many-electron 
ground states associated with rational filling factors. Nevertheless, in order to 
generate the fractional plateaus in the Hall conductivity, disorder is believed to 
be essential in pinning the chemical potential, in close analogy to the Integer 
Quantum Hall Effect. 



1.3 Plan of the Article 

In summary, the Quantum Hall Effect — consistent with all of the presently 
available experiments — seems to be an exact and very fundamental phe- 
nomenon. One needs for its physical understanding, and in particular for 
a theory with predictive power, not only interactions and correlations in a 
two dimensional many-body system but in particular, and somewhat counter- 
intuitively - - in view of the precision aspect - - the presence of a certain 
amount of disorder. It is one of the main purposes of this article to review our 
present understanding, why, in spite of and due to the disorder, the Quantum 
Hall Effect is the most exact phenomenon known in low dimensional elec- 
tron systems. We will restrict ourselves to the Integer Quantum Hall Effect, 
ignoring interactions completely. 

Especially, we will discuss the localization concept by introducing the so-called 
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random network model originally invented by Chalker and Coddington [31] in 
order to describe the localization of quantum states in a random potential 
with long-range spatial correlations in the presence of a strong, "quantizing" 
magnetic field. Later, we will discuss several generalizations of the network 
model which underline the general importance of the phenomenon. 

We will proceed in the following chapters by pursuing a strategy which we hope 
will enable non-specialists in the field to understand the main ideas in a self- 
contained way. We will always start from elementary accessible facts that form 
the background of more complicated relations before we eventually describe 
to some extent the issues that are more at the frontier of present research. As 
the field has been, and still is, extremely rapidly growing, we cannot guarantee 
that all of the available results have been included, especially the most recent 
ones. We have attempted our best to include at least references to all of the 
important results. We apologize in advance to all colleagues, whose research 
efforts might not be sufficiently highlighted according to their importance or 
even not be included at all. 

The plan of the paper is as follows. In the next chapter, the so-called per- 
colation model for an electron in a spatially slowly varying random potential 
landscape is introduced in some detail. It can be viewed as the backbone of 
the random network model. If it is supplemented by quantum effects, it can be 
considered as a precursor of the random network model of Chalker and Cod- 
dington for describing localization in the presence of a quantizing magnetic 
field. 

The latter will be introduced in chapter 3. In chapter 4 the localization proper- 
ties of the model will be discussed in some detail. Strangely enough, technical 
difficulties seem to have prevented a high-precision determination of the criti- 
cal exponent up to now for the original Chalker-Coddington random network 
model. Chapter 5 contains some results for the transport properties, especially 
near the quantum critical point. 

In the chapter 6, we investigate the renormalization group approach for de- 
scribing the critical properties in some detail. Especially, we will emphasize 
here that one can use this as a starting point for the definition of a trun- 
cated Chalker-Coddington model which can be treated numerically exactly 
and allows for determining the critical exponent with high precision. 

Chapters 7, 8, 9 and 10 contain extensions to a number of equivalent Hamil- 
tonians, to systems with other symmetries, and to an equivalent field theo- 
retical formulation of the random network model of Chalker and Coddington, 
respectively. In these chapters, we have attempted to provide some flavor of 
the power and importance of the model in fields other than the quantum Hall 
phase transition. While we will review much numerical and experimental ev- 
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idences for the universality of the quantum Hall transition, any attempts to 
classify it, like other two dimensional phase transitions, based on the conformal 
invariance at the critical point have failed so far. 

In chapter 10 we review a very recent suggestion, namely that the quantum 
Hall critical point may belong to a new class of critical points being described 
by a supersymmetric conformal field theory [60,61]. This is based on the fact 
that an anisotropic version of the Chalker-Coddington model as well as the 
random Landau model at the critical point can both be mapped on the Hamil- 
tonian of a chain of antiferromagnetic superspin chains. We review its deriva- 
tion, and the progress which has been achieved towards the characterization 
of the quantum Hall transition that way, although an analytical calculation of 
its critical exponents is still missing. 

Generalizations to several layers and higher dimensions are briefly discussed 
in the chapter 11 before we conclude by summarizing the status of the field 
and comment on possible future developments. 



2 The Percolation Model 



We start by describing the physics of the so-called percolation model for the 
Quantum Hall Effect. In this model, which is valid in the limit of a very high 
magnetic field, the quantum mechanical wave functions are assumed to perco- 
late along the equipotential lines of a slowly varying random potential land- 
scape, in analogy to a classical percolating fluid in a random system [62]. Using 
this picture of a percolating wave function, one can understand why localized 
wave functions may exist, and, in particular, why there may be isolated crit- 
ical points in the energy spectrum of the Hamiltonian where the localization 
length diverges. This high-field limit is often denoted as the classical perco- 
lation limit, although what is percolating is probability amplitude and not a 
classical fluid. In the quantum mechanical version of the model, tunneling of 
probability amplitude between the equipotential lines is taken into account 
whenever they get close to each other in space. In this version, the necessary 
competing ingredients — quantum tunneling and interference — are included 
such that a generic universal quantum mechanical localization-delocalization 
transition can be described. 

The percolation model provides the physical background of the Chalker-Cod- 
dington network model — a generic model which is assumed to describe the 
universal quantum mechanical properties of non-interacting electrons in two 
dimensions in the presence of a random potential subject to a strong perpen- 
dicular magnetic field. 
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2. 1 Wave Functions in a Spatially Correlated Random Potential in a Strong 
Magnetic Field 



The origin of the random potential in the plane of the inversion layer V(x, y) 
of the electrons in a GaAs/AlGaAs heterostructure or a MOSFET transistor 
is the impurities in the semiconductor material, especially as a consequence 
of the doping. These impurities are distributed randomly at some distance 
from the inversion layer. Thus, only the Coulomb tails of their potentials do 
influence the dynamics of the charges. This implies that the distribution of 
the potential energy V(x, y) is long-range correlated in space. For convenience, 
we assume in the following that the spatial average of the random potential 
vanishes. This can always be achieved by a suitable choice of the zero of energy. 
For the correlator we assume 



with the correlation function C (x, y) having an exponential or Gaussian de- 
cay length £ c much longer than the magnetic length £b- It is also assumed 
that the probability distribution of the potential, P[V(x, y)\ is symmetric and 
homogeneous, i.e. independent of the origin of coordinate system. 

Figure 4 shows an example of such a long-range correlated random potential, 
together with several examples of the corresponding wave functions. They 
were obtained by solving numerically the Schrodinger equation 



in the basis of the Landau states. An important characteristic feature to be 
kept in mind is that the wave function amplitudes are essentially non-zero only 
along equipotential lines of the potential. In addition, when two equipotential 
lines get close to each other near a saddle point of the potential, the amplitude 
becomes high across the saddle point. This can be taken as an indication of 
quantum tunneling as indicated in Fig. 4 by arrows. 

The qualitative features of the wave functions in Fig. 4 may be taken as the 
motivation for introducing the concept of a percolating probability amplitude 
which we will now describe. The resulting model is important because it al- 
lows for qualitative understanding of localization in two dimensional electron 
systems in the presence of a strong magnetic field. It allows also to understand 
why, within each Landau band, there must be an energy where the localization 
length diverges. This energy corresponds to the percolation critical point at 
which the wave function can percolate throughout the entire system. The basic 
physical ingredients of the percolation model may then be used to establish 



V(x>, y>)V(x, y) := W 2 C(x' -x,y'-y), 



(18) 



[H + V(x, y)] i/j u (x, y) = E v ip v (x, y) 



(19) 
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(a) (b) 




(c) (d) 

Fig. 4. Long-range correlated random potential obtained as superposition of ran- 
domly placed Gaussian potentials with a width i c = 21b and some examples of 
eigenfunctions. (a) gray scale plot of the potential landscape with equipotential 
lines indicated; white: high potential, black: low potential; (b)-(d) the moduli of 
a selection of characteristic wave functions corresponding to (b) a low energy in 
the tail of a Landau Band, localized in a deep potential valley; (c) an intermediate 
energy in the Landau band, still localized and following mainly an equipotential 
line; (d) an energy near the center of the Landau band, extending essentially along 
equipotential lines, occasionally inter-connected via tunneling near the saddle points 
of the potential (indicated by arrows). 
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as an idealized version the random quantum network which allows for system- 
atically dealing with certain universal localization and transport phenomena, 
especially in the quantum Hall critical regime. 

In order to introduce the essential idea, we consider the Schrodinger equation 
(19) in the limit l c ^> Ib that can always be achieved for very large magnetic 
field, B — > oo [63]. In this limit, we will eventually also fulfill the condition 
hjjB W such that the disorder-induced mixing of the Landau bands can be 
neglected. 

It is of advantage to introduce new coordinates [5,63,64] 

X = -kP B , Y = -it* B ±:=*jLp. (20) 

They are called center-of-motion coordinates. This is reasonable since X is 
the center of the Gaussian wave packet in the x-direction of the Landau states 
(Eq. (6)). Furthermore, one observes that (id/dk) p plays the role of (y p ) by 
calculating the matrix elements of y p in the Landau states (p integer). At 
finite magnetic field, the electrons can be viewed as describing cyclotron or- 
bits with radius £b around the center of motion. For £ B "C £ c the diameter 
of the cyclotron orbit becomes vanishingly small such that (x, y) can be re- 
placed by center of motion (X,Y). A The center-of-motion coordinates fulfill 
commutation relations similar to position and momentum operators, 

[X,Y]=iP B . (21) 

We start by writing the solutions of the Schrodinger equation for the system 
of the size L in the representation of the Landau states 

*(x, y) = Y, C n (-i 2 B k) ip nk (x, y) := $>y|nX> C n (X) . (22) 

nk nX 

This gives 

£ (E n 5 x , x/ + (nX\V\nX')) C n (X') = EC n (X) (23) 

X' 

with inter-Landau band couplings neglected as a consequence of the high field 
limit, index k on the energy E n of the Landau level omitted due to degeneracy, 



4 Note that this limit has to be taken only after the thermodynamic limit L —> oo 
has been performed which is necessary to make X a continuous variable. 
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and the Kronecker delta 5x,x'- The potential matrix elements are 

(nX\V\nX') = ^J dxdy (f) n (x - X) V(x, y) <f) n {x - X')e iy{x - x ' )/e B . (24) 



Since the potential is slowly varying on the scale of £b it may be expanded 
into a power series near some arbitrary yo(x) in Eq. (24). With this, one easily 
verifies that 

J dy V(x, y ) e Mx-x')/el = y ^ -if B ^j ± J dye ~^~x')/fi B (25) 



Finally, the integration with respect to x may be approximated by assuming 
x ks X in the potential since at high magnetic field \<f> n (x — X)\ 2 can be 
assumed to be very well localized near X within a distance £b- Now Eq. (23) 
becomes 

V(X,Y(X'))C(X')\ X=X , « EC(X). (26) 



Here, E is the energy corresponding to the equipotential line, and the index n 
has been suppressed. Furthermore, E n has been assumed to be the zero of the 
energy. This is justified since the distance between the Landau bands increases 
oc B and the disorder broadening of the bands is only oc \[B such that at high 
field the couplings between the bands due to disorder can be neglected. 

The random Schrodinger equation (26) can approximately be solved by using 
the semi-classical WKB (Wentzel-Kramers-Brillouin) Ansatz 



C(E,X) ocexp 



x 



— i J dX'k(E,X') 



(27) 



neglecting the derivatives of the local wave number k(E,X) with respect to 
X. This is possible for sufficiently smooth randomness such that second order, 
and higher order derivatives in the expansion may be omitted. Under these 
assumptions, the dispersion k(E,X) fulfills the equation 



E — V 



X,£ 2 B k(E,X) 



(28) 



This establishes a most remarkable result, namely that the eigenfunctions can 
be considered as the superpositions of the Landau states that are associated 
with the equipotential line associated with the energy E. 
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2.2 The Form of the Semiclassical Wave Functions 



The general form of these wave functions which propagate on the equipoten- 
tial lines may be conjectured by considering the elementary example of a two 
dimensional harmonic potential V(x, y) = Vq{x 2 + y 2 ). This represents approx- 
imately the situation far from the Landau band center near a local minimum 
or maximum of the potential. Of course, this model can be solved exactly. 
However, it is instructive to look at it using the above percolation viewpoint. 



The equipotential lines are the circles X 2 + £ B k 2 (E ,X) = E/Vq 
applying Eq. (26) and calculating C(X) from Eq. (27) one obtains 



e. By 



C(X) oc exp 



l€r 



X 



2£\ 



arccos 



iX 

W B 



x 5 



(29) 



The energy is given by the condition that C(X) must be a periodic function 
of X. This implies for the closed integral over the equipotential circle 



<f dX' k(E, X') = := 2ixm (m integer). 

J In 



(30) 



This is consistent with the result one obtains for the lowest Landau level in 
the symmetric gauge with the definition Vo := e 2 £>Q/8m*, and m being the 
angular momentum quantum number [65]. 

By inserting the result Eq. (29) into the expression for the wave function 
Eq. (22) one gets for £ 2 B <C Hub/Vo, i-e. neglecting the inter-Landau level cou- 
pling, an expression for the wave function which is a reasonable approximation 
for X/ v /7^< 1, 



ip(x, y)oc J dX exp ^2y ± ^Je m -X 2 



exp 



{x - Xf 
2£l 



x exp 



it. 



2£ 2 B 



^Pm(X) 



(31) 



with <f m (X) = arccos(X/ v /e^). Approximate evaluation of the integral gives 



if)(x, y) oc exp 



ixy 



exp 



(32) 



The wave number corresponds to k m := m/^fe^ and u = v^v^™ ^ s ^ ne az " 
imuthal coordinate on the circle with the radius JT^. The phase factor me- 
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diates a gauge transformation to the symmetric gauge. The Gaussian factor 
localizes the wave function along the circle within an interval of the width 
« 2y/2£ B . 

This suggests that in general the wave functions corresponding to the equipo- 
tential lines can be written in the form [66,67] 

^x,y)^f{v)e^"> (33) 



with a local wave number k(u, v) that depends on coordinates u and v which 
parameterize the distance along and perpendicular to the equipotential line 
V(x, y) = E, respectively. The function f(v) is non-vanishing only within a 
region of the approximate width £b along this line. Locally, the wave func- 
tion has the form of a wave propagating along the equipotential line. This is 
associated with an equilibrium current density 

j oc ^e u \f(v)\ 2 hK(u,v) (34) 

that yields a net current along the equipotential line 

hn{u,v) 

3 oc —e u - 35 

m* 



Perpendicular to the direction of u, the current density obviously must vanish 
(Fig. 5). Classically, the net current is produced by the cyclotron motion of 
the electron which drifts along the equipotential line in the strong magnetic 
field and the local electric field of the potential that is directed perpendicular 
to the equipotential line (Fig. 5). This visualizes that only superpositions of 
those Landau states can contribute to the eigenstates which are located near 
the equipotential line at the corresponding eigenenergy. 

The quantum mechanical problem of solving the Schrodinger equation of a 
particle in a high magnetic field and a random potential has thus been re- 
placed by the task of finding the equipotential lines in the random potential 
landscape. This corresponds to a classical percolation problem [68]. Recently, 
another connection between quantum Hall plateau transitions and percolation, 
based on the classical limit of quantum kinetic equations, has been discussed 
[69,70]. 



2.3 Localization in the Percolation Limit 



The localization properties in this limit of a very strong magnetic field may 
now be easily discussed [66,71]. Consider the smoothly varying landscape of 
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(a) 



(b) 



Fig. 5. (a) Scheme of a potential profile near the local potential minimum, and 
(b) the current density of a state localized along an equipotential line. The current 
density circles within a region of diameter 21b around the equipotential line. Only 
the net current flows along the equipotential line. 

the random potential (Fig. 4). When the energy is very low (or very high), the 
corresponding equipotential lines are closed trajectories captured within the 
minima (or maxima) of the potential (Fig. 4b). When the energy increases, 
these closed trajectories begin to meander around several of the minima (or 
maxima) (Fig. 4c). 

The wave functions that correspond to these equipotential lines are necessar- 
ily localized superpositions of Landau states. They are localized exponentially 
as has been shown in [72]. Only at a certain critical energy, say E c , close to 
the center of a Landau band, an equipotential line can percolate through the 
whole system. Only at this energy the wave function can propagate through- 
out the entire system in the limit of infinite system size (Fig. 4d). The energy 
E c corresponds to the percolation threshold of the classical percolation prob- 
lem. For symmetric distribution of the potential, P(V) = P(—V), the critical 
energy is at the center of the Landau band, E c = 0. 

Thus, the percolation picture allows immediately for a very important conclu- 
sion: since the percolating equipotential lines are closed for all energies except 
the critical energy of the percolation threshold, all of the wave functions must 
be localized except for the one associated with the critical energy. 

Let us define the localization length £ P (E) as the correlation length of a per- 
colating equipotential line. Then, percolation theory says that £ P (\E — E c \) 
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must increase according to a power law as \E — E c \ decreases, 

U£ - £c) oc ^-p- , (36) 



with a universal exponent v v = 4/3 that has been calculated exactly [73]. 
2.4 Tunneling Correction to the Percolating Wave Functions 



When two equipotential lines get close to each other near a saddle point of the 
potential, tunneling processes will occur. In order to determine the localization 
length, these have to be taken into account [74] . For a symmetric distribution 
of the potential, the tunneling, however, will not change the position of the 
critical energy. 

A general definition of the localization length in terms of the Green function 
is [17] 

7(E) 1 = _ lim my-m 

/V ; £(E) |r-r'|-oo \r-r'\ V ; 



with (...) denoting the ensemble average with the above probability distri- 
bution of the potential, P. It is well known [75,76,77,78,79], that in weakly 
disordered systems, defined by kpl ^> 1, where I is the disorder induced mean 
free path and fcp the Fermi wave number, the disorder averaged electron wave 
function amplitude (ip{x : t)) decays on length scales of the order of /, since 
the random phase shifts associated with the scattering at the impurities are 
averaged out. This destroys the information on multiple scattering, and thus 
on localization. In order to describe localization, it is therefore necessary to 
average over disorder functions containing higher moments of the propagator 
G, like the expression used in Eq. (37). Furthermore, it is well-established 
that this quantity is self-averaging, namely its ensemble average coincides 
with its most probable value [17] in the thermodynamic limit. Therefore, one 
can calculate instead of the ensemble average the spatial average for a given 
realization of the random potential. We will now use this definition of the lo- 
calization length for estimating the correction to the critical exponent, u p , of 
the percolation model due to tunneling at the saddle points of the potential. 

Using the above implicit dispersion relation, Eq. (28), we first determine the 
Green function in the mixed representation (X,k(E ,X)) where Eq denotes 
the energy of an equipotential line, 

G(XMEo,Xy,E)= E _ Eo( l x) _. n , (38) 
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with r] a positive infinitesimal real number. By Fourier transforming with 
respect to k one can determine the behavior of the Green function along the 
y-axis which is sufficient for estimating the asymptotic behavior at \y\ — > oo, 

G(X, y;E) = ^-Jdk e~ ik y G(X, k; E) , (39) 



one obtains 

G(X, y;E) ocJ2 e- ik ^ E ' x)y = J2 e~ iRck > y ' e -^(E,x)\ y \ ^ 



where kj(E,X) are the complex roots of Eq. (28) for given X and E = E . 
The asymptotic exponential decay of the Green function for \y\ — > oo is given 
by the spatial average of the smallest lmkj(E, X) > 0, 

l(E) = (imn j lmk j (E,X))p. (41) 



As described above (compare Fig. 4), we need to consider saddle point regions 
of the random potential near E c , since predominantly it will be here where 
the wave functions are connected between different equipotential lines via 
tunneling, and for the critical behavior, energies close to E c are important. 
Near a saddle point (X , fc ), the potential can be expanded (Fig. 6), 

V(X, e B k) = E c -a 2 (X- X ) 2 + b 2 £ 2 B (k - k ) 2 . (42) 

/From E = V(X : £ 2 B k) one obtains the solution 

k - k = ±^-jE -E c + a\X - X ) 2 . (43) 

bin v 



Due to symmetry, it is sufficient to consider, say, E < E c . Then, Irak ^ only 
for E c -E> a 2 (X - X ) 2 

lmk(E, X) = A-jE c -E-a?{X-X Q y. (44) 

01 b 

Averaging this with respect to X yields 

7 (E) oc J dX ^E c -E- a 2 (X - X ) 2 oc E c - E (45) 



with X± = X ± VE C - E/a. 
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Fig. 6. Equipotential trajectories of a saddle point potential for an energy below 
the saddle point potential. Dashed lines: equipotential lines exactly at the energy of 
the saddle point; shaded: region with potential energy higher than the saddle point 
energy E c (after [80]). 

The asymptotic exponential decay for energies near E c can be estimated by us- 
ing the following heuristic argument. The equipotential lines representing the 
wave functions form percolation clusters with an energy dependent correlation 
length £p(E) ~ \E — i? c |~ 4 / 3 . Derealization of the wave functions associated 
with these clusters can occur via tunneling through the saddle points that 
connect different equipotential clusters. The number of saddle points will be 
proportional to the size of the clusters at energy E, namely oc £p(E). Thus, 
the probability of finding the saddle point with the smallest 7 at energy E 
connecting two adjacent clusters will be £b/£p(E). The saddle point with the 
smallest 7 will dominate the asymptotic behavior of the wave functions. All of 
the contributions of the other saddle points will be exponentially suppressed 
with the diameter of the clusters. Therefore, we finally can write for the aver- 
age inverse localization length the power law 

<* TT^\ E - E d <x\E- £ c | 4/3+1 • (46) 



According to this qualitative argument, the critical exponent of the localiza- 
tion length in a quantum Hall system including the effect of quantum tunneling 
should be v = 7/3, without interactions. 

It has taken considerable theoretical [19,20,23,81,82,83] and experimental [39,40,41, 
efforts to determine quantitatively this exponent for two dimensional disor- 
dered Landau system without interaction. The most accurate numerical values 
obtained so far are v = 2.35 ± 0.03 by using a random Landau matrix model 
[20], v = 2.33 ± 0.03 from a Monte Carlo calculation for chiral fermions [38], 
and a renormalization group approach based on a random network model 
v = 2.39 ±0.01 [34,35,36] (cf. Tab. 1). Numerical investigations were done for 
a variety of completely different models, including white noise [19,23] as well 
as long-range correlated [31,82] randomness and including also higher Landau 
bands [20,83]. All of the results are consistent with these values within error 
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bars. This strongly suggests that the critical exponent of the quantum phase 
transition associated with the Quantum Hall Effect is indeed universal, and 
does not depend on the microscopic details of the model used 5 . Although the 
near coincidence of this exponent with the value 7/3, as obtained from the 
simple heuristic argument outlined above, seems compelling we caution that 
an analytical derivation of the critical exponents at the quantum Hall transi- 
tion is still lacking. Recent progress towards achieving this goal is reviewed in 
chapter 10. 



3 The Random Network Model 

The results obtained in the above discussed high-magnetic field limit may 
be used not only for estimating the critical behavior but also as a starting 
point for constructing a model that contains all of the necessary physical 
ingredients - - backscattering and tunneling — for describing the quantum 
critical points near the centers of the Landau levels and the corresponding 
transport quantities. If the critical behavior could be shown to be universal 
such a model should be of great importance. This would be especially true, 
if one could use the model as a starting point for more rigorous theoretical 
formulations. 

As argued above, percolation of probability amplitude along equipotential lines 
as well as tunneling of amplitude between equipotential lines near the saddle 
points of the potential have to be taken into account on an equal footing when 
determining critical properties. This suggests to construct a model consisting 
of a regular lattice of saddle points that are connected via links along which 
probability amplitude can propagate. 



3.1 The Scattering Wave Functions Associated with a Saddle Point 

Before doing so, it is instructive to consider saddle point tunneling more for- 
mally from the scattering point of view, and establish an exact expression for 
the transmission probability. As a side remark, we note that this is yet another 
illuminating example of a quantum wave function that can be treated exactly 
in the high magnetic field limit [80]. 

The most simple saddle point is defined by the quadratic potential considered 



5 For very long range correlated randomness there recently seems to be evidence 
that v changes [84]. 
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in the previous chapter with equal coefficients a 2 = b 2 = U [80] 
V sp (x,y) = E c + U(y 2 -x 2 ). 



(47) 



The classical trajectories in such a potential for energies below E c are shown 
in Fig. 6. For the corresponding quantum mechanical wave functions the con- 
siderations of the previous sections apply when neglecting tunneling between 
the trajectories. This can be justified for energies well below (or well above) 
E c . However, when E E c this is no longer the case. 

The total Hamiltonian including the saddle point 

H = H + V sp (x, y) = -L (- V + eA] + V sp (x, y) (48) 



is quadratic in the variables. It can be diagonalized in analogy with in the 
case of the two dimensional harmonic oscillator. 

However, here we follow the very illuminating method used in [80]. For the 
vector potential, we assume the symmetric gauge, A = (B/2)(—y,x). We 
perform first a unitary transformation that brings the Hamiltonian into a 
separable form. 

By introducing the new variables 



that fulfill the commutation relations 



[a x , 4] = [ay, a\] = 1 , (51) 
[a x ,Oy] = [al,al]=0, (52) 



the Hamiltonian becomes 



H = —{a l x a x + a l y a y + 1) + —(ala y - a y a x ) + 



+l[{a y + a\) 2 -{a x + al) 2 } + E c 



(53) 
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where 


7 = 






' a x 1 


-( 









2 cos <z> sin < 



I 


'I 







(54) 



with tan(20) = —Mje/^I, transforms the Hamiltonian into a sum of two 
independent terms, Hi(bi,b\) and H 2 {b2,b\), 



= &l) 



H2 = (b\b 2 ) 



( *f- - Ml 7 



V -7 ^ + Ml 







) 


;:;) 
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7. 



where Ml = y 7 2 + (M>_b/4) , and 

[&i,&S] = [&2,&5] = 1, 
[61,62] = [6i,4] = 0. 



(55) 



(56) 



(57) 
(58) 



These can be diagonalized by a Bogoliubov transformation, 
/ 



bi 
6j 



cosh#j sinh^j 
sinhft cosh6>- 
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) 


7 



(59) 



with tanh(20i) = 7~ 1 (-M> B /4 + Ml), and tanh(2# 2 ) 
which gives for the Hamiltonian 



7(^/4 + Ml)- 1 , 
H = E l (c\ + cf) + E 2 (c\c 2 + £) + £ c , (60) 



with 



[ci,c{] = [c 2 ,c 2 ] = 1 , 
[ci,c 2 ] = [ci,c|] = 0, 



(61) 
(62) 



and the energy eigenvalues 
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(63) 
(64) 



Finally, we introduce variables 



X = M C ^~ Cl ^ s =T7f( c 2 + c 2 ), 



p = ^(«i + c 1 ), p = ^(c 2 -4), 

and obtain 

H = H 1 + H 2 = E 1 (p 2 - X 2 ) + -E 2 (p 2 + s ; 



(65) 



(66) 



with the commutation relations 



[X,P\ = [s,p\=i, 
[s,X] = [s,P} = [p,X] = [p,P] = 0, 



(67) 
(68) 



Thus, apart from the factor £ 2 B /h'm Eq. (20), H 1 corresponds to the center-of- 
motion part of the total Hamiltonian, and H 2 is a one dimensional harmonic 
oscillator which implies that the wave function is harmonically confined in the 
direction of s. 



The eigenfunctions of the Hamiltonian can now be factorized 
*(X,s)=<j>(X)Ms), 



(69) 



where ip n corresponds to the n-th harmonic oscillator level. 

In order to arrive at a scattering solution of the Schrodinger equation, we 
prepare a initial wave packet at (X 2 ) ps (P 2 ) 3> 1. Since 

^- = y/2( ai X-M, j- = V^i^P + a 2 p) (70) 



with 



COS < 



A = -sin0e^ (i = 1,2), 



(71) 
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we have asymptotically 
Oil Pi 



(72) 



Since Hub 7, Pi/c^i ~ —1; and the wave packet is centered near an asymp- 
tote of the equipotential line of V sp in the upper left quadrant of the (x, y)-plane 
(Fig. 6). 

In order to find the transmission probability, we need now to construct a 
scattering wave function. The eigenfunction satisfies 



H^ftX) = E^P 2 - X 2 )<f>(X) = 



E - [ n + - ) E 2 - E c 



0(X), (73) 



or equivalently 
/ d 2 



dX 2 



+ X 1 + e <j>(X) = 



(74) 



with the energy parameter 

E-{n + 1/2) E 2 - E c 



e = 



E 1 



(75) 



In the limit of high magnetic field, E 2 = huJs and E\ = 7. Therefore, e 
measures the energy deviation from the saddle point energy E c in a Landau 
band normalized by a typical potential strength. 

This Schrodinger equation is discussed in detail in [85]. The eigenfunctions 
even and odd in X are 



$+= e~ ix2 / 2 F 



<$>_=Xe- ix2/2 F 



n + ie 


1 
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where F is a confluent hypergeometric function. For large \X\ they can be 
written in the form 



£(1/2) 
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r(3/2) 



e -i{./2)z*_ f(z*_-X)Xe lX2 l 2 + c.c. 



T{Z_) 

with the definitions f(z;X) = \X\~ 2z , Z_ = (3 + ie)/4, Z 4 
The currents in these states are related to the derivatives 

9 ±ix 2 /2 = ±X e ±ix2/2 . 



(l + ie)/4. 



(78) 



Therefore, each first term in Eqs. (77) corresponds to a currents in the positive 
X-direction. The conjugate complex term represents a current in the negative 
X-direction. 

These solutions are linearly superposed, t = A<3>+ + !?$_, to form the scatter- 
ing wave functions which must fulfill boundary conditions such that for X > 
only a transmitted wave function exists with the current flowing away from 
the origin. This leads to the relation 



T(Zt) 



(79) 



One finds after some algebra for the transmitted wave for large values of X < 



A r(l/2) c -^ /8 | B T(3/2) c _ i37r/8 
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and for the incoming wave 

UX) = f{Z+,X)e-™ / 8 e-^ 2 /2 
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r(z;) r(z*) 



(81) 



This guarantees that for X > no incoming wave exists. The transmission 
probability is defined by 



T(e) = lim 



X^oo|0i(-X)|2- 



(82) 



One gets from the above equations by using the relations for the T-functions 

V2n 



T(zy = T(z*), rQ + iy)r(?-iy 



cosh ny + i sinh 7ry ' 



(83) 
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the final result 



T(e) 



1 



(84) 



1 + exp (— Tie) 



For e — > — oo, well below the saddle point energy, the transmission proba- 
bility vanishes. The incoming wave is completely reflected. For e — > +oo, the 
transmission probability approaches unity. Exactly at the saddle point energy, 
e = 0, one has T(0) = 1/2. For energies near the saddle point, one can expand 



3.2 Parameterizing the Scattering at a Saddle Point 

As we have seen, the percolating wave functions of the electron carry an equi- 
librium current density. Near the points where percolating wave functions 
approach each other closely, quantum tunneling takes place. 

A natural model for the description of the percolating quantum states can 
then be obtained by considering a network of current loops, occasionally inter- 
connected via tunneling near the saddle points of the potential (Fig. 4). The 
simplest choice is a regular network. Disorder can be introduced either by 
randomizing the relative phases of the current loops and/or the heights of the 
saddle points. The distribution of the wave function amplitude is then closely 
related to that of the corresponding quantum currents and the asymptotic 
behavior is essentially determined by the coherent interplay of the quantum 
transmission and scattering near the saddle points. 

Let us consider the region close to a specific saddle point, say E c . The poten- 
tial landscape can be modeled by two neighboring adjacent potential maxima 
and minima (Fig. 7). We assume that the directions of the currents along the 
equipotential lines in the potential minima and maxima are clockwise and 
counter clockwise, respectively. Let us consider energies close E c . There are 
four equipotential lines that enter the region of the saddle point. The total cur- 
rent entering and leaving the saddle point region must be conserved (Fig. 7). 

The effect of the saddle point potential can be considered as a quantum trans- 
mission problem as described above [31]. In more simplified terms, we assume 
that the amplitudes which correspond to the currents on the boundary of the 
saddle point region are ipi . . . tp4 with the amplitudes of incoming and outgo- 
ing currents ipi, ip 3 and ip2, respectively. The formal relationship between 
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(a) 



(b) 



Fig. 7. (a) Potential landscape near a saddle point with adjacent maxima (full 
equipotential lines) and minima (dashed equipotential lines), (b) The corresponding 
idealized potential landscape with incoming (ipi and ^3) and outgoing (ip2 and ^4) 
current amplitudes flowing along equipotential lines. Note that the saddle point 
energy is denoted here by Vq. 



them can be described by a scattering matrix S, 





(86) 



In general, the S matrix can be written as [86,87] 




(_ rt \ i 



t r 



e i<Pi 
e"^ 3 



(87) 



where 0<r<l,0<£<l, r 2 + t 2 = 1, and r 2 and t 2 are the probabilities of 
the current ipi to be reflected to ip2 and to be transmitted to ^4, respectively. 
Note that the current conservation relation implies 



^ir+i^3i 2 = i^r+iv>4i 



(88) 



which is always valid due to the unitarity of the scattering matrix S. 

For iterative numerical calculations, it is useful to introduce the transfer ma- 
trix T that is equivalent to S. It relates the amplitudes on the left, (^i, V ; 2) T , 
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to those on the right hand side of the saddle point, (^4, vp 3 ) T (T denotes the 
transposed). 



(89) 
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Using Eqs. (86) and (87) one finds the result 



T = 
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e^ 3 
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(90) 



For completeness, we only mention here, that due to the current conservation, 
the 2 x 2-matrix T must fulfill the symmetry condition (cf. chapter 9) 



J = T f JT 



with 




(91) 



It is useful to relate the present description to the microscopic saddle point 
model introduced above. By using Eq. (91), and imposing the correct asymp- 
totic behavior from the scattering geometry in Fig. 7, one can write the re- 
flection and transmission amplitudes r and t in terms of a single parameter 



cosh 6 



(92) 



and by the current conservation relation t = \/l — r 2 , 
t = tanh O . 



(93) 



Equation (90) can be written in terms of the variable 0, 
\ / 
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O' and 6 being related via the equation 
sinh O sinh 0' = 1 . 



(95) 



By comparing with Eq. (84) one finds that the dimensionless parameter 
must depend monotonically on the energy of the equipotential line. It charac- 
terizes completely the transmission properties of the saddle point. For E <C E c 
and E ^> E c , <C l(r ps 1) and ^> l(t ps 1), respectively. This is because 
for £ < £ c , f ! ps ^2 for almost all -^3 and ^4, and for E ^> E c , ipx m vp^. 
One can determine the energy dependence of by comparing with the above 
microscopic result near e = where T(0) = 1/2 := \t\ 2 = tanh 2 C and ex- 
panding for Ps C since the model is only suitable for the region near the 
saddle point energy, E c . One finds (cf. Eq. (85)) 

9(e) ^ 0c+ 2^ + '" (96) 



with the saddle point value C = ln(l + a/2). 

The characteristic feature of this model for the transmission through a sad- 
dle point is that incident and transmitted channels are locally separated. The 
interplay between the saddle point potential and the high magnetic field intro- 
duces a spatial separation of the incoming and outgoing channels. This makes 
the model particularly suited for describing the critical localization features 
in the quantum Hall region. 



3.3 Establishing the Random Network of Saddle Points 



Starting from the smooth random landscape of the potential, the network 
model could now be defined from a random system of circular, localized wave 
functions of the type Eq. (33), equivalent to circular equilibrium currents. 
These would correspond to states associated with randomly distributed sites. 
For technical convenience, however, it is preferable to assume the circular 
wave functions to have random phases and being associated with the sites of a 
regular lattice. In addition, it is reasonable to assume that nearest-neighbored 
wave functions are connected by tunneling contacts described by the above 
transfer matrices T (cf. Eq. (90)). This enables the particles to hop from 
one circular state to another. A derealization mechanism for the total wave 
function is introduced in this way. 

The model constructed is the analogue at high magnetic field of a two dimen- 
sional Anderson model for a disordered system [17,88,89]. The site states of 
the latter are replaced by the circular currents of the former, and the hopping 
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amplitudes by the unitary transfer matrices. As in the Anderson model, the 
network contains the necessary competing ingredients — localization within 
the circular current states and tunneling between them — for describing a 
localization-delocalization phase transition. In the case of the Anderson model, 
there is no phase transition in two dimensions. As we have already seen above, 
the present network model must show a singularity in the localization proper- 
ties which represents a quantum phase transition point. This will be discussed 
in more detail in the next chapter. 

But let us first complete the model by providing a more formal description, 
especially suitable for numerical work. We consider a two dimensional rectan- 
gular geometry with 2M x 2L current loops on a square lattice (Fig. 8). The 
structure of the current flow near the nodes described by S' is rotated by n/2 
as compared with that of the nodes described by S (Fig. 8). We have 





(97) 



with the 7r/2 rotated scattering matrix 
S' = 




(98) 



that results in a corresponding transfer matrix T' relating (*0 4 , *0 3 )* to ^2)*, 
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(99) 




cosh O sinh 
sinh cosh 6 




(100) 



The transfer matrix that relates the amplitudes on the left of the sample 
(^0,1, i>o,2, ■ ■ , ^o,2lY to those on the right (ip M ,i, i/>m,2, ■ ■ • , ^m,2lY (L integer) 
is then defined as follows. We divide the system along the direction of M into 
slices, each containing L of subsequent saddle point scatterers and the 2L 
corresponding channels (Fig. 8). The latter are assumed to have completely 
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Fig. 8. Left: The network model in the quantum Hall regime: + and — indicate the 
hills and valleys of the potential, respectively. The scattering at the saddle points 
is characterized by S matrices, S and S'. The arrows indicate the directions of the 
guiding center motion. Right: The scattering described by S' differs from that of 
S by ir/2 rotation in the incoming and outgoing channels; L width of the system, 
M length of the system in the direction of the repeated application of the transfer 
matrix. 

random phases. The amplitudes at the end of the system can then be written 
in terms of a transfer matrix 6 



( Mi ^ 



1 ^0,1 ^ 
^0,2 



(101) 



with the product of the transfer matrices 



T M = r p( M ) , j i ( M - 1 ) . . . r p( 2 ) r j i ( 1 ) 



(102) 



We use here the somewhat unconventional notation that M and L denote the 
length and width of the system, respectively. The reason is that when discussing 
the numerical scaling in the next chapter, we prefer to use L for the length scaling 
instead of M. 
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where 

T(*)=vf ) V 2 V? ) V 4 . 



(103) 



Here the matrix V4 is given by 



(V 4 )y = 



cosh 6' i = j = n, 

i = j + 1 = 2n, 

sinhO' < 

% = j - 1 = 2n - 1, 
otherwise 



n = 1,---,2L 

n = 1,2, 



(104) 



and V 2 by 



(V 2 )i 



cosh 6 i — j — n, 

i = j + l = 2n+l, 
i = j — 1 = 2n, 



sinh < 



n = 2, 3, • • • , 2L - 1 



n = 1,2,---,L- 1 



(105) 



Assuming periodic boundary conditions in the transverse direction gives 



(V 2 )i,i = (V 2 ) 2 l,2L = coshG 
(V 2 )i,2L = (V 2 ) 2 l,i =sinh9 



(106) 



while for fixed boundary conditions one has to use 

(V 2 ) M = (V 2 ) 2L , 2L = 1 
(V 2 )i, 2L = (V 2 ) 2£i i = 0. 



(107) 



All other matrix elements of V 2 are 0. 

When propagating from one node to the other, the probability amplitudes 
gain phase factors. This is contained in the matrix elements of vf^ and V^, 



(V) A ' ! )„, V- (Z = l,3). 



(fe> 



(108) 



Since the distance between the nodes is random, we assume that (j)f^ are 
independent and uniformly distributed between [0, 2ir). So far, this is the only 
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source of randomness in the model. Note that the phases <p and tp' in Eqs. (90) 
and (99) can be included in the 0's in Eq.(108). 

In this form, the transfer matrix method [17,88,89] has been applied to the 
model for estimating the critical behavior of the localization length in the 
quantum Hall regime [31] to be described in more detail in the following 
chapter. 

Attempting some completeness, we mention at this point that some versions 
of the network, which essentially reproduce features of the underlying classical 
percolation problem, and with random saddle point potential values E c have 
been used to obtain detailed information about the longitudinal conductivity 
and several classical percolation properties [32,90,91,92]. A classical version 
of the regular network consisting of coupled metallic wires has been used to 
calculate the Hall conductance. Quantization in integer (positive and negative) 
units of e 2 /h has been predicted consistent with earlier findings for the pure 
quantum Hall case without any disorder [93,94]. The model has been also 
used to study the localization problem in two dimensions in the presence of 
a random magnetic field [95,96]. As our emphasis will be on the universal 
features of the network model in the forthcoming chapters, we will not go into 
the details of these works. 



4 The Localization-Delocalization Transition in the Network Model 

In this chapter we will explore the localization properties of the original 
Chalker-Coddington network introduced in the previous Chapter 3. As this 
has been a key issue for establishing the model, this will be done in some 
detail. Necessarily, for obtaining quantitative information about the critical 
behavior, numerical methods will be used. We attempt to provide a complete 
overview of the statistical properties of the energy spectrum and the wave 
functions near the critical point, as far as it is presently available. 

The model has been the subject of numerous numerical studies which are 
all similar in spirit but different in the details [32,97,98,99,100]. We provide 
here results obtained recently by the numerical scaling method [88,89] which 
have been undertaken in order to improve the treatment of corrections to 
scaling. For obtaining reliable and precise results for the critical exponent, 
the latter has been shown to be decisively important [83,101]. To our great 
disappointment, however, as we will show below, up to now it has not been 
achieved to remove corrections to scaling to such a degree that the precision 
of the exponent of the Chalker-Coddington network model can compete with 
the results obtained for other models in particular the random Landau model, 
see Tab. 1! 
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We will also discuss the critical properties of the wave functions as well as the 
eigenenergy statistics. Due to its simplicity, the network model enables us to 
investigate them in detail, and because of universality, they are supposed to 
be generally valid for two dimensional electron systems in high perpendicular 
magnetic fields. 

To study the localization-delocalization transition, we investigate the scaling 
properties of the localization length in a quasi-one dimensional long strip for 
which the localization length is calculated by the transfer matrix method and 
apply the finite size scaling method developed earlier for estimating the asymp- 
totic value in the thermodynamic limit. The value of the critical exponent v in 
the network model [31,32,102,103] is consistent with, though considerably less 
accurate than, earlier results obtained for several very different models such 
as the random Landau model with a white noise potential [23,24] and the ran- 
dom matrix model including spatial correlations of the randomness [19,20]. 
This supports strongly that the model of Chalker and Coddington belongs to 
the same universality class as the models considered previously, in spite of the 
intriguingly large corrections to scaling which defy the high accuracy in the 
estimate of the critical exponents. 



4-1 Numerical Scaling at the Localization- Delocalization Transition 

To determine the critical behavior quantitatively, we use the numerical finite 
size scaling approach. We define a quantity F L which is a function of a set of 
parameters Xi characterizing the system, such as the Fermi energy, parameters 
characterizing the disorder, as the variance and the correlation length, the 
magnetic field, and the system size L [104], 



We assume the existence of a scaling law such that F L can be expressed as 



with x being the relevant scaling variable and fa denoting the irrelevant ones. 
The latter nevertheless can cause corrections to scaling as long as the system 
size is finite. Therefore, in a numerical calculation they must not be ignored. 
These variables characterize distances from the critical point, v is the critical 
exponent and yi < are the exponents of the irrelevant scales. Eventually, 
in the limit of infinite system size, only the relevant scaling variable survives, 



F L = f({ Xi },L). 



(109) 




(110) 
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and Eq. (110) becomes 

A -*(f). 



mi) 



with £ ~ x v ■ The quantity % as a function of some control parameter, say x, 
can be expanded near the critical point x c 



Which quantity should be used as the scaling variable Fj?. It should be a quan- 
tity that shows a singularity in the limit of infinite system size as one crosses 
the localization-delocalization critical point. Furthermore, it should be a quan- 
tity to be determined numerically easily, precisely and effectively. There are 
several such quantities, such as for instance the level spacing distribution [105] 
or the conductance (see below). Here, we consider the renormalized localiza- 
tion length in the finite system (MacKinnon-Kramer variable), A, introduced 
previously [88] which we will now define. 

Consider a very long strip of width L. This is a quasi-one dimensional sys- 
tem and all the states are expected to be localized. The localization length 
A(L; xi, x 2 , . . .) is a function of the width of the strip, L. Now let us define 



If the system parameters are such that X(L — > oo;xi,x 2 , . . .) remains finite, 
the system is in the localized regime and A(L — > oo; xi, x 2 , . . .) — > 0. The 
localization length is then given by £ = X(L — > oo;xi,x 2 , ■ ■ ■)■ On the other 
hand, if X(L — > oo; x±, x 2 , ■ ■ ■) increases faster than oc L, the system is in 
the delocalized state, and X(L — > oo;xi,x 2 , ■ ■ .), the correlation length in the 
metallic regime, corresponds to the inverse of the dc-conductivity [88,89]. The 
critical point is defined by the condition A(L — > oo; x± c , x 2c , . . .) = const = A c 
(critical MacKinnon-Kramer variable). 

The most efficient way to calculate X(L) is the transfer matrix method [17]. 
We define the product of the transfer matrices Eq. (102) and consider the limit 

T = lim (Tf Tf t ) 1 / 2M . (114) 

M— »oo 

The theorem by Oseledec [106,107,108,109] guarantees that T has always pos- 
itive eigenvalues, which are denoted as exp(±7j) where 7i(> 0) can be inter- 
preted as the exponential change of the wave function for a single slice. The 



X = Xi(x- x c ) + X2(x - x c f H . 



(112) 



A(L;xi,x 2 , . . .) 



X(L;x 1 ,x 2 , ■ ■ ■) 
L 



(113) 
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smallest value of 7$ is the inverse of the quasi-one dimensional localization 
length [17], 

^y = min{ 7i }. (115) 



This definition can be shown to be equivalent to the one in terms of the 
exponential decay of the Green function which could also be used to study 
A(L) [110]. If one is interested only in the localization length the transfer 
matrix method has been shown to be superior. 



4-2 Numerical Results near the Quantum Critical Point 



The two dimensional electron systems in high magnetic fields are character- 
ized by many control parameters such as Fermi energy, magnetic fields as 
well as parameters describing the properties of randomness. In the Chalker- 
Coddington model, however, all the informations are contained in a single 
parameter, the transmission at the saddle point. In the actual simulation, it 
is convenient to use as control parameter 

x = — InsinhO, (116) 



With this choice of x, we obtain from Eq. (92) 

11 . . 

t = - r= — = . 117 

Ve 2x + 1 Ve~ 2x + 1 



Due to the particle-hole symmetry (r <-> t) the localization length is an even 
function of x. Comparing Eqs. (84) and (117) we note that x can be interpreted 
as the energy measured from the center of the Landau band. 

In the actual simulation, we need to simulate not too small systems so that 
only a single relevant scaling variable and at most one irrelevant variable are 
sufficient for fitting the data. In this case, the scaling form Eq. (110) reads 

A(L)=F( X L^,0^), (118) 



where x is related to x via Eq. (112). We consider in the calculation the region 
x < 1 so that x ~ x an d <t> = const. 

As noted above, due to particle-hole symmetry A is an even analytic function 
of x as long as L is finite. By assuming xL l l u to be sufficiently small in order 
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to truncate the expansion after the second order, Eq. (118) can be expressed 
as 

A(L) = A c + a 1 (xL 1 ' v ) 2 + a 2 (xL^) 4 + b L y + c L 2y + ■■■, (119) 

Previously, this non-linear fitting scheme has been working perfectly well for 
the three dimensional Anderson transition [101,111,112,113] as well as in the 
case of two dimensional systems with strong spin-orbit coupling [114,115]. 
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Fig. 9. MacKinnon-Kramer variable A(L) = \{L)/L as a function of x for 
L = 6,8,12,16,24,32 and 48. The relative uncertainty of the data is 0.02%. The 
curves are results of the fitting according to Eq. (119) with A c = 1.23, u = 2.48 and 
y = —0.53. The fitting, however, is unstable, so the curves should be regarded as 
guides to the eye. 



High precision data are required for a precise determination of the exponent. 
As seen in Fig. 9, for such data the corrections to scaling are not negligible 
[102]. Unfortunately, a stable fit that takes into account these corrections has 
not yet been found [116,117]. As a result, there has not been any significant 
improvement of the precision of the estimate of v over the original estimate 
of v = 2.5 ± 0.5 by Chalker and Coddington [31]. 

Fluctuations of the transmission properties of the saddle points cause another 
type of randomness in addition to the phase randomness assigned to the wave 
function when traveling from one node to another. These fluctuations are 
equivalent to randomness in the mass of the Dirac Hamiltonian to be described 
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later in the chapter 7. It has been numerically demonstrated that this effect 
gives larger — but irrelevant — corrections to scaling [102]. 



4-3 The Critical Properties of the Wave Functions 



Exactly at the critical point t — r — l/y/2 (x = — InsinhG = 0), the wave 
function is delocalized. It is now well-established that due to the divergence 
of the length scale, the wave function at the critical point shows multi-fractal 
behavior [118,119,120]. This self-similar structure is in principle reflected in 
an anomalous temperature behavior of the diffusion constant, the energy level 
statistics, the conductance distribution, and many other properties. 

Since the Chalker-Coddington model is characterized by the scattering matri- 
ces at saddle points rather than by the Hamiltonian, the wave function must 
be determined from the scattering matrices. Let the amplitude on the link 
I be ipi. As explained above, the stationary states of the model must satisfy 
[121,122] 

= tmk^k + t m l1pl, (120) 



where t m k and t m i relate the amplitudes ipk and ipi to ip m (Fig. 10). This can 
be rewritten as 

\J(E)V = V, (121) 



with XJ(E)ei = t m i(E)e m + t n i(E)e n , being the unit vector with the ith 
component unity, |t mfc | 2 = T = 1 — \t m i\ 2 := 1/2 at the critical point, and 
\I/ the vector with the components ip m . Again, we assume for simplicity that 
randomness enters only via the phases of the amplitudes. 

Equation (121) has non-zero scattering solutions only for certain discrete en- 
ergies E u , which determine the eigenstates and eigenenergies of the system 
[123]. Solving Eq. (121), however, is in general very complicated for a disor- 
dered system. We will derive an alternative method by establishing equivalent 
Hamiltonians later in the chapter 7. 

For the moment, we proceed by calculating instead the return probability [124] 
from the equation of motion determined by Eq. (121). We start by introducing 
the time to during which the incident wave packet passes the scatterer and 
measure time t in units of to- We can then define the time evolution operator 

U(t) = XJ\E) (122) 
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Fig. 10. The stationary wave function relation characterized by the scattering matrix 
amplitudes. Incoming waves and outgoing waves are ipi,tpk and ipm&n, respectively. 



where t is an integer. 

The probability P(t) of return within the period of time t is determined by 
the probability density of the wave packet at time t near a given scatterer, 
say located at the origin, which was previously located there at t — 0. This is 
defined by 

P(t) = |(0|U(t)|0)| 2 . (123) 



The return probability is related to the fractal dimension D(2) defined by the 
second moment of the probability density, the inverse participation number 
P2(L) 

P2(L):= <|*(r)| 4 ><x^^. (124) 



with (. . .) denoting a configurational average. If we integrate P2{L) over the 
two dimensional system and then take the inverse, we have L D ^ and we 
can interpret this quantity as the portion of space where the wave function 
amplitude is significant. 

Let the wave packet spread to the radius r(t) after time t, the number of 
sites occupied becomes ~ r(t) D ( 2 \ On the other hand, the conductivity at the 
quantum Hall transition is finite, and from the Einstein relation 

a = e 2 pD , (125) 



(p the density of states per volume and D the diffusion constant), we see that 
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D is finite, hence r(t) ~ t 1/2 . One then finds [125,126,127,128,129] 

P(t) oc r^ 2 . (126) 

In general dimension the diffusion radius r is proportional to t 1 / d [130], hence 
P(t) oc r D{2 ^ d . (127) 

For the spatial dimension d — 2, we can estimate .0(2) = 1.52 ± 0.06 from 
numerical simulations [124]. This is in agreement with the results in the con- 
tinuum model [118,131] as well as in the tight binding model [126]. 

A similar idea has been developed to calculate the local density of states [132] 
from which one can calculate the multi- fractal exponents D(q) that are defined 
by the higher moments of the density, 

Pq (L) := <|*(rr> oc oc — ^ , (128) 

The moments of the density characterize the degree of localization of the 
wave function. For plane waves, p q (L) oc L~ 2q — > in the thermodynamic 
limit. For localized states p q (L) oc L 2q — > const ^ 0. Generally, the fractal 
dimensionalities D(q) depend on q and r(q) is non-linear in q. This is the 
celebrated multi-fractal behavior [133]. It means that the complex structure 
of wave function at the critical point can not be described by single dimension 
-0(2) but infinite number of generalized dimensions are required to characterize 
it. 

A convenient quantity to summarize the behavior of the multi-fractal expo- 
nents is to introduce a Legendre transformation on r(q) by defining [134,135,136] 

a = ^ (129) 
aq 

and the /(a)-spectrum 

f(a) = aq-r(q). (130) 



The physical meaning of the quantity f{ct) can be understood by relating 
it to the probability distribution of the random amplitudes of P$(|^| 2 ). 
Starting point is the identity between the latter and the distribution, P a {a), 
of the random variable 

ln|^| 2 . . 
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P a (a)da = P*(|#| 2 )d|#| 2 . (132) 

The inverse participation numbers are obtained in terms of these distributions 
p q (L) = J d|^| 2 l^l^d^l 2 ) = J dae~ aqlnL P a {a) . (133) 



In order to obtain the correct dependence of p q {L) on the system size the 
distribution of a must be of the form 

P a (a) oce [/(a) ~ 2]lnL . (134) 



For very large L, the average in Eq. (133) is dominated by the maximum of 
f(a) at «o- Thus, «o is the exponent of the system size dependence of the 
typical value of the probability, 

m 2 f yP (L) ■= e< ln ('*' 2 » oc -L (135) 



Figure 11 shows some numerical results for f(a) together with the result 
obtained from the Dirac model with random vector potential ([120] see chapter 
7.3 below) which is exactly parabolic 

/(a) = 2 -4(^2) (136) 



with a = 2 + A^/V where is the variance of the random vector potential 
which is in that case equivalent to the random phases associated with the links 
of the network. 

Again, this is in agreement with the diagonalization analysis [131]. Together 
with the analytical results from the Dirac model, which will be described be- 
low, this provides strong evidence for the overall conjecture that the Integer 
Quantum Hall Effect is a quantum critical phenomenon. A graphical represen- 
tation of a multi-fractal wave function at the critical point is shown in Fig. 12. 



The fractal dimensionality of the wave function can be related to the MacKin- 
non-Kramer parameter A c by assuming conformal invariance [31,139,140,141,142]. 
The quantity a , is related to A c by (see Sect. 5.3 for more details) 

«„ = 2 + -L (13 7) 



which is in excellent agreement with the numerical results [102,137,142,143]. 
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Fig. 11. The /(a)-spectrum at the quantum Hall transition determined numerically 
(full line) and the result of an analytical conjecture (dashed line) with ao = 2.262. 
Inset: Spectra of typical eigenfunctions for different system sizes L and extrapolated 
with L -> oo [from Ref. [137]]. 

4-4 Energy Level Statistics near the Critical Point 



Random matrix theory is a very powerful tool to characterize the proper- 
ties of complex systems [144,145]. It can also be applied to the localization- 
delocalization transition. In the metallic regime, the overlap of the wave func- 
tions is strong, which results in strong level repulsion. In the limit of large 
system size, the energy level spacing s (in units of the mean level spacing) is 
given approximately by the distribution function 

P(s) oc s p e~ A ^ . (138) 



This is often denoted as the Wigner surmise. The value of the parameter (3 
(= 1,2,4) is determined by the symmetry of the system. The value of the 
constant A depends on f3. When the system has both time reversal and spin 
rotation symmetry, (3 = 1. This is called the orthogonal symmetry class. If 
the system has only time reversal symmetry but spin rotational symmetry is 
broken by spin-orbit interaction, (5 — 4. This characterizes the symplectic class 
where level repulsion is strongest. If (5 — 2, time reversal symmetry is broken, 
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Fig. 12. Squared amplitude of a critical wave function in a Chalker-Coddington 
network of 256 x 256 saddle points. Darker areas denote lower square amplitude. 
No characteristic lengths scale can be identified [Figure taken from [138]]. 



irrespectively of whether or not spin rotational symmetry is present. This is 
the unitary symmetry class. Symmetries classes are indeed very important 
ingredients for characterizing the universal properties of a quantum phase 
transitions. We will discuss this in more detail later in the chapter 9. 

On the other hand, when the states are localized, the correlations between 
eigenenergies vanish in the limit of large system size. Then, the level spacing 
distribution P(s) is Poissonian, P(s) = exp(— s). When the system size is 
finite, the correlations between energy values are still present and the spac- 
ing distribution deviates from the Poissonian. The deviation from the limit 
of infinite system size can be described by using a single parameter scaling 
assumption Eq. (Ill) [30,146,147,148,105] 




(139) 
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Fig. 13. Scheme of the quasi-energies uj v as a function of E. The condition lo v = 
gives the eigenenergies, E v . 

where £ is the localization length. For the ordinary metal-insulator transition, 
the above equation has two branches corresponding to metallic and insulating 
phases. In the Chalker-Coddington model, we do not have a true metallic 
phase. Therefore, we expect to find only a single branch. 

Exactly at the critical point, £ diverges and P(s) becomes size independent. It 
is neither a Poissonian nor it corresponds to the Wigner surmise, but has both 
characteristics. It grows as s 13 for s <C 1 as in the case of the Wigner surmise 
and it decays according to exp(— As) for s ^> 1 like the Poisson distribution 



We now want to obtain the energy level statistics from Eq. (121). To achieve 
this, we need a method to extract information about energy eigenvalues from 
the transfer matrix approach. This requires the construction of an equivalent 
Hamiltonian which we will do below in great detail. For the present purposes 
it is sufficient to consider the following. 

First we note that the eigenvalue equation of the unitary operator \J(E) is 



The eigenenergies E^s that correspond to the stationary states are obtained 
from the condition exp[iu u (E u )] = 1. 

The statistics of the energy levels obtained from the condition uo{E) = is the 
same as that of levels obtained by uo = Q, since the latter correspond to the 
solution for U'(E) = e~ tn U(E) which belongs to the same universality class. 
This together with the level repulsion leads to the conjecture that 0^(0) obey 



[105,149,150]. 



U(E)* V = e 



(140) 
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the same statistics as E v [151]. The quantities u v (0) are often called quasi- 
energies. This assumption greatly simplifies the numerical calculation. The 
results, for example the form of P(s) at the critical point, agree in fact very 
well with those obtained for the continuum model [152] . Due to the simplicity 
of the network model and the possibility of using the quasi-energy concept, 
the level statistics can be investigated in detail. 

Another important quantity to characterize a random sequence of energy levels 
is the so-called number variance defined by [144] 

£ 2 (A) = ((n - A) 2 ) = (n 2 ) - A 2 (141) 



where n is the number of levels in a randomly chosen interval AA (A average 
level spacing), (• • •) denotes the configurational average, and A = (n). In the 
insulating region, £ 2 (A) = A. In the metallic region, it increases only loga- 
rithmically £ 2 (A) ~ In A due to the level repulsion that makes the spectrum 
rigid. At the metal-insulator transition, 

ton = X, (142) 

/V— >oo J\ 



with < x < 1 [147]. 

The value x — 1 implies that the system is an insulator with no level repulsion, 
and x = is equivalent to metallic behavior with maximal level repulsion. The 
fact that < x means that the level repulsion is weakened as compared to the 
metallic limit. This is due to the fact that the wave function has a very sparse 
multi-fractal structure. Thus, x reflects the multi-fractal behavior of the wave 
functions. In fact, the quantity x an d the fractal dimension D{2) are related 
via [153] 

d-D(2) 

X = -jjU (143) 



Inserting D(2) = 1.52 ± 0.06 and d = 2 yields x = 0.120 ± 0.005. Klesse 
and Metzler have estimated x from the quasi-energies, x = 0.124 ± 0.006, in 
agreement with Eq. (143) [151]. This agreement, however, is not exact since 
Eq. (143) holds only approximately as has been discussed in Refs. [154,155]. 



5 Linear Electrical Transport at Zero- Temperature 



The linear electrical conductance of a quantum system at zero temperature 
is related to quantum mechanical transmission via the Landauer formula 
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[156,157]. Thus, the network model is perfectly designed to provide quanti- 
tative information about the linear conductance. If the model describes the 
physics near the critical point, one can expect also that it is especially suit- 
able for the critical conductance. The study of the linear conductance tells 
us the importance of the conductance distribution instead of the averaged 
conductance, which converges to a form independent of the system size. The 
qualitative behavior of the conductance distribution is expected to be valid 
even in different situations such as 4-terminal conductance measurement. 

In this chapter we want to discuss the quantum conductance and its distribu- 
tion function at the critical point at absolute zero of the temperature. For the 
conductance, we use a slightly modified random network model (Fig. 14 and 
[158]). 




Fig. 14. The network model as modified for calculating the two terminal conduc- 
tance. The amplitudes Vh an d V'oi (* = 1> ' ' ' represent the incoming and the 
outgoing currents on the left hand side, respectively, an d V'oi (* = 1>"""?-^) 
represent the amplitudes on the right hand side of the sample. 



5. 1 The Transfer Matrix and the Conductance 

The general wave function on the left hand side of the system is a superposition 
of the incident and reflected amplitudes, {ip^} and {ip^i}, respectively These 
are related to those on the right hand side, {ip^} and {V^j}, via the transfer 
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matrix of the system of the lengths M and the width L, T^f , 



1 <i ^ 



= T 



M 



<1 



(144) 



By introducing an operator U which changes the order of the components, 



u 



\A7 



Ui R J 



= u 



^i R i 

V ^ J 



(145) 



one can write 



R 



= U f TU 



(146) 



where = R, L) denotes the vector of amplitudes (^0,1 > Vo,2> ' ' ' ■> Vo,l) T an d 
= R, L) is (^i, ipl 2 , ■ • ■ i ^i,L) T ( T denotes the transposed vector). 

The linear zero-temperature conductance G is given by the quantum trans- 
mission probability through the system [156,157,159]. 

e 2 

G=-Trtt t . (147) 



Thus, a relation between the transfer matrix and the transmission matrix tij 
= 1...L) is required. In the following paragraph this relation will be 
established by writing the wave function on the left and right hand sides of 
the scatterer in terms of L x L transmission and reflection matrices t and r 
[160,161]. 
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When a flux of probability amplitude is injected into the ith scattering chan- 
nels from the left and from the right, represented by ^ and ififi, respectively, 
the total scattering wave function in the ith channels is a superposition of the 
incident wave and the reflected and transmitted waves from all of the channels. 
This is described by the scattering matrix S 




(148) 



with 



S = 



t r' 



(149) 



By solving Eq. (148) for ijj^ and ipf and inserting into Eq. (146) one obtains 
straightforwardly 



U f TU 



r't'^r r't'- 1 ^ 
-t'^r t'- 1 




(150) 



In order to determine the (L x L)-matrix t' that one can use equivalently 
for the conductance instead of t, it is sufficient to generate Ti 2 and T 22 by 
multiplying T to the 2L x L matrix which consists of the zero-matrix Ol in 
the upper half and the unit matrix 1l in the lower part, 





-* 











(151) 



Multiplying iteratively (k = 0, 1, 2, 3, . . . M) is numerically unstable. 

Thus, after each m steps one must perform a QR decomposition, 

rj-i (i x m+m.) rj-i (i x m+m— 1) _ _ _ rj-i(ixm+l)-yi yi+l^i (152) 

where V* are (2L x L) orthonormalized matrices and oo 1 are (L x L) upper 
triangular matrices. This eventually yields 



(153) 



56 



with matrices u l generated during the iteration and [L x L)-matrices A and 
B the upper and lower blocks of y M l m . 

Using this result, one finds for the transmission and reflection matrices from 
Eq. (150). Details and generalization can be found in [162]. 

t' = (w 1 )- 1 ^ 2 )- 1 ---^*)-^- 1 , r' = AB 1 . (154) 
5.2 The Critical Conductance and its Statistics 

The conductance G is a strongly fluctuating quantity with a broad distribution 
function. Therefore, not the conductance but its distribution function P(G) 
has to be considered. 

In the metallic regime, the conductance is described by a normal distribu- 
tion with its variance independent of the details of the system as well as of 
the system size. This is known as the phenomenon of reproducible universal 
conductance fluctuations (UCF [163,164]). In the insulating regime, the distri- 
bution function has a log-normal shape, reflecting the exponential localization 
of wave functions. At the critical point of the localization-delocalization tran- 
sition not only the variance but the whole distribution function becomes size 
independent. This reflects again the scale invariance of the quantum critical 
point [111,165,166,110,167,168]. Remarkably, however, the critical distribution 
depends on the boundary conditions [112]. 

In the case of the quantum Hall transition, all the states except those at the 
band center are localized, so that the conductance distribution is log-normal. 
As expected, at the band center the conductance distribution function becomes 
size independent. However, it has the peculiar form shown in Fig. 15. It neither 
is consistent with the log-normal distribution expected in the localized regime 
nor reproduces the non-universal behavior of the moments predicted for the 
metallic region [169]. It also depends strongly on the boundary conditions. 
This is consistent with earlier results obtained for the level statistics at the 
Anderson transition [170,171]. 

iFrom the distribution, one can calculate a configurationally averaged conduc- 
tance (G) and its moments [172]. We consider here as an example only the two 
terminal conductance for periodic boundary conditions. We expect, in addi- 
tion to corrections to scaling, effects of the contacts [173] and of the boundary 
conditions [112]. Both of them can be expected to give rise to contributions 
proportional to L -1 . Therefore, we attempted to fit the average conductance 
to the scaling Ansatz 

(G) = G c + aL y + y. (155) 
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Fig. 15. Conductance distribution at the critical point. The number of incoming 
links is 128. The geometrical shape of the system is a square. Left: for periodic 
boundary conditions; right: for fixed boundary conditions. See also [158] 

This fit yields the bulk (L — > oo) conductance 



G c = (0.570 ±0.02) 



h 



y = -0.56 ±0.05 



(156) 




Fig. 16. The system size dependence of the "bulk conductance" at the critical point. 
Solid line: fit to the scaling law Eq. (155), G c = 0.5702 - 0.2150L 



-0.5587 



0.1486 



The form of the distribution of the critical conductance is similar to that ob- 
tained from simulating the conductance of the tight binding model in magnetic 
fields [174]. However, there are quantitative discrepancies. For example, the 
two terminal conductance in units of e 2 /h is 0.506 for the tight binding model 
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[169], but 0.57 for the Chalker-Coddington model. This might be due to the 
inter-band coupling present in the tight binding model. 

In Fig. 15, we observe a continuous increase of P(G) from near G — 0, and 
a kink near G = e 2 /h. The kink is expected in one dimension [175,176], in 
the two dimensional symplectic ensemble [167,177] and in three dimensional 
disordered systems [168] . It is related to the fact that for one transport channel, 
one can achieve at best unit transmission. Fluctuations of the conductance to 
values larger than e 2 /h are related to the contributions from the exponentially 
smaller contributions of the other transport channels. 

For better understanding the kinks in the distribution functions, the trans- 
mission eigenvalues {rj} obtained by diagonalizing tt^ have been analyzed in 
more detail [172]. It has been found that for the quantum Hall transition the 
critical conductance distribution is well approximated by taking into account 
only the largest transmission eigenvalue as demonstrated in Fig. 17. This re- 




0.2 0.4 0.6 0.8 1 1.2 G(e 2 //))orxmax 



Fig. 17. The distribution of the critical conductance (x) compared with the distri- 
bution of the largest transmission eigenvalue r max (+) at the critical point. (After 
[172]) 

suit suggests that the knowledge of the distribution of r max would enable us 
to approximately predict the transport coefficients at the quantum Hall tran- 
sition. 



In order to further investigate P(r max ), we transform to a new variable u min 
according to 

2 

W = — r — r • (157) 

cosh z/ min + 1 



It has been suggested in [178] that z/^ in is Poissonian distributed, 

P(i£j = PNOe-e™-*** (158) 
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or equivalently, 

P(z/ min ) = 2f3Nnv min e-? Nn ^. 



(159) 



In Fig. 18, we show P c (f m in) and fitted the data to Eq. (159) by assuming 
(3Ntt = 0.233. 

P(Vmin) 
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Fig. 18. The distribution of v m i n compared with the Poissonian distribution 
P(fmin) = 2ai/ min exp(-azA n ) with a=0.233 [After [172]]. 

/,From the very good agreement between -P(zVm) an d the Slevin-Nagao form 
[178] one might conclude that the distribution of {Vj} is given by 

PK u 2 , u N ) =C\{e-^'i T7 \ v \ - v )\\ (160) 

i i<j 

which defines the Laguerre ensemble [178,179]. However, the distribution of 
second and higher values, i.e. z/j {i = 2,3, ..) as calculated from the Laguerre 
ensemble [180] deviate from the result of the present numerical simulation. 
This is not unexpected since the Laguerre ensemble has been proposed to be 
a good approximation for the statistics of the metallic diffusive regime. At the 
quantum Hall transition, the so-called eigenvalue "interaction" term may be 
modified, 

PK u 2 , u N ) =C\{e-^^ 1] f(Wf ~ (161) 

i i<j 




This form leads to the same distribution for z/ mm but different distributions 
for higher v { {i = 2,3,...) [172]. 

Before concluding this section, we mention that the energy correlation function 
of the two terminal conductance has been analyzed in detail by Jovanovic and 
Wang [103]. 
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5. 3 Conformal Invariance at the Critical Point. 



We have presented above strong numerical evidence for the scale invariance 
of the wave functions, energy level statistics and the conductance distribution 
at the critical point. Under a scale transformation, the length scale is rescaled 
globally by some factor b such that the coordinates of a given point transform 
according to 

r' = rV. (162) 



Scale invariance at a critical point means that correlation functions of scaling 
variables <f>j(r) (j = 1, 2, 3 . . .) are invariant under such a transformation which 
generally can include also a rotation and a translation 

Mri)fo(r 2 ) ...} = n^Wr'M) • • •) • (163) 

3 



The exponents hj are called scaling dimensions. 

In analogy to classical phase transitions, one can expect not only that global 
scale invariance holds but also that the more general concept of conformal 
invariance [139,104] applies to the quantum Hall transition [181,31]. Generally, 
a conformal transformation corresponds to local translations, rotations and 
dilatations of the coordinates which preserve angles. 

In order to explain this in more detail, we consider the above general correla- 
tion function in two dimensions, (<f>i(zi, z*)02(<22, z 2 ) • • •), with the scaling vari- 
ables <f>j (j = 1, 2, 3 . . .) and the spatial variables represented as Zj = Xj + iyj 
and its complex conjugate z*. Now let us consider an arbitrary analytic map- 
ping 

z = w(z) . (164) 



Locally, close to some point z , which is equivalent to the infinitesimal map 

z' - z' Q = w'(zq)(z - zq) , (165) 



with w'(zo) = dw(z )/dz . Comparing with Eq. (162) one notes that the global 
scale factor b^ 1 of the scaling transformation corresponds to w'(z ). 

In analogy to Eq. (163), conformal invariance of the correlation functions can 
be defined by 
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(Mzi,4)Mz2,z* 2 )...)= (166) 

n™'(*)*M*n%i&, 002(4, 4*) • • •) • 

i 

The exponents /ij and h[ are real- valued conformal scaling dimensions [104]. 

In order to show how powerful the concept of conformal invariance is for 
extracting properties of the correlation functions we consider the two-point 
correlation function of only one scaling variable. Conformal invariance requires 

(<t>{z u zl)<t>{z 2 ,z* 2 )) = \w\z 1 )\v\w\z 2 )r(<f>(z[,z[*) ( f ) (z' 2 ,z' 2 *)}. (167) 



It is a straightforward exercise to show that this implies a power law decay of 
the correlation function [104], 

(0(^,^)0^,^)) ~ \z x -z 2 \-^. (168) 



As a second example, we consider the mapping between the infinite two di- 
mensional plane and the surface of cylinder with circumference L, 

z' = —\nz. (169) 
2ir 



This is motivated by our above numerical studies in which we have used sys- 
tems of finite width L and length M — > oo. It is easy to get convinced that 
variation of x, y in the complete two dimensional plane implies — oo < x' < oo 
and < y' < L. This corresponds to an infinitely long strip of width L with 
periodic boundary conditions in the direction of L which is equivalent to a 
cylinder. The correlation function on the cylinder reads 

(0(4,4*)0(44*)>~ (170) 

(2tt/L) 2t > 

{2cosh[2vr(x , 1 - x' 2 )/L] - 2 cos[2?r(yi - y 2 )/L}}^ ' 

Therefore, the quasi-one dimensional exponential decay length £ cy i and i] are 
related via 

&, = ^. (171) 



If we assume cf>(z) to be the typical density of states defined via the square 
amplitude of the wave function, 

P tyP = (172) 
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where AE ~ L 2 is the mean level spacing, we obtain from Eq. (135) 

7] = a - 2 . (173) 



The transfer matrix estimates the exponential decay length by averaging 
log |^(0)^(,z) I . Therefore, it is the typical amplitude of the wave function which 
is related to the decay length. Combining Eqs. (171) and (135) 

^M^Jy < 174 > 



which is the same as Eq. (137) since A c = 2£ cy i/L. This example illustrates how 
powerful the requirement of the conformal invariance is. Therefore, whether 
or not this symmetry holds for a phase transition is very important and can 
provide useful information about the properties near the critical point. 

An important transport coefficient which can be used to clarify whether or not 
the conformal invariance applies is the point-contact conductance [182,183]. 
Consider a network with two links / and m cut from the interior, and connected 
to two reservoirs (Fig. 19). Then we can define a 2 x 2 S matrix 

( Su Sim 

o = 




and the point-contact conductance is given by (in units of e 2 /h) T = \Si m \ 2 . 
The S matrix is calculated by 

Si^m-UPtPJ-'Ulj) (i,j = l,m) (176) 



where U is the time evolution operator introduced by Eq. (121) and P; = 
1 — |/)(/| and P m = 1 — \m)(m\ the projection operators which describe the 
nature of reservoirs (Fig. 19). 

The point-contact conductance T is related to the wave function intensity at 
/ and m, \^i\ 2 and |^ m | 2 via [183] 

27rp( J E)(|M/ m | 2 /(|^ m | 2 /l^| 2 )) = (F(T)) , (177) 



where 

2w 



F(T) := I — /(T _1 |l - e^VT^l 2 ) ■ (178) 
J 2ir 
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Fig. 19. Average half-integer moments (T n+1 / 2 ) of the two-point contact conduc- 
tance of the Chalker-Coddington model at the critical point as a function of n. Data 
sets (top to bottom) correspond to distances 5,10,15,20 (in units of the lattice con- 
stant) in a cylinder geometry with total length M = 100 and circumference L = 10. 
Solid lines are analytical results obtained by conformal invariance. For best fit an 
exponent Xt = 0.54 ± 0.01 was found. Inset: schematic view for the model. Two 
links, I and m, are cut and connected to reservoirs I and m (after [183]). 



Here p{E) is the density of states and (• • •) means the average over randomness. 
Assuming f(x) = — lnx, we have 

(lnT) = 27rp( J E)(|M/ m | 2 ln(|^| 2 /l^r))- (179) 



Thus one can evaluate the typical conductance T typ = exp((lnT)) very effi- 
ciently using only the information of the wave functions in an isolated system. 

At the critical point, in the infinite two dimensional plane T typ is expected to 
decay as 

T t yP _ r -x t ^ ( lg0 ) 



where r is the distance between the links I and m and X t a critical exponent 
related to r(q) (cf. Eq. 128) via r'(0) - r'(l) [184]. 

^From Eqs. (171) and (180), the decay of the typical conductance on a long 
cylinder surface becomes [183] 



2vr 



sinh 



/ irx\ 



\LJ 



-x t 



(181) 
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where x is the distance between the links / and m along the cylinder axis. 

Numerical results for the moments of T(x), (T(x) n+1 / 2 ) (Fig. 19) and the 
typical conductance T^f(x) (Fig. 20) coincide with the behavior expected 
from the conformal invariance [182,183]. This strongly supports the conjecture 
that the wave function at the quantum Hall transition is conformally invariant 
[181]. 




Fig. 20. Logarithm of the typical two-point conductance (InT) of the 
Chalker-Coddington network model at the critical point as a function of the log- 
arithm of the distance ln|x| for system widths L = 3,4,6,9,14,20,30 (bottom to 
top). Solid lines: analytical result obtained using conformal invariance, Eq. (181). 
Inset: rescaled curves (InT) + X t lnL versus ln|x/L| with X t = 0.57 (after [183]). 



6 The Re normalization Group Approach 



In this chapter, we discuss the application of the real-space renormalization 
group approach that is widely used in classical percolation and spin systems 
[104,62] to the random network model. A closed set of renormalization group 
equations is derived which describes the universal distribution of the conduc- 
tance. The numerical solution is used to estimate the critical exponent of the 
localization length and the behavior of the moments of the distribution func- 
tion [33,185,186,187]. Amazingly enough, by adopting a slightly different point 
of view, eventually one can interpret the approach as providing an additional, 
independent precision determination of the critical exponent. The agreement 
of this critical exponent with the earlier results obtained for completely dif- 
ferent models confirms that the quantum Hall phase transition is indeed a 
universal critical phenomenon. 
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The general idea is that near the quantum critical point all microscopic details 
of the system eventually must become irrelevant since the correlation length 
diverges with a universal exponent which, however, still can depend on the 
fundamental global symmetries. In order to detect such universal behavior 
one must perform a well-controlled thermodynamic limit. An example how to 
do this in a controlled way has been given above when applying the transfer 
matrix method and using the numerical scaling method. This procedure is 
basically exact and suffers only from technical numerical errors, but it does in 
general not provide analytical insight. 

In classical phase transition theory, the renormalization group transformation 
is a well-established procedure to approach the critical point. Generally, this 
method consists of a sequence of unitary and subsequent scale transformations. 
The former are used to diagonalize relatively small systems. The latter scale 
the system back to the original length scale. From the behavior of the coupling 
parameters under this repeatedly applied transformation conclusions can be 
drawn on the critical behavior. In principle, this method would also be exact, 
apart from the fact that in each renormalization step couplings to states that 
are energetically far away from the critical point are neglected. While this 
seems not to be crucial for showing the very existence of a critical point, 
one needs quite substantial numerical efforts in order to correctly obtain the 
critical behavior quantitatively. 

There are many alternative possibilities for constructing the renormalization 
group transformation which are all similar in spirit but may be different in 
the details. We will discuss here two of them. The random network model of 
Chalker and Coddington is an example in which the renormalization group 
transformation can most advantageously be used. 



6.1 An Illustrative Example: the Tile Lattice 

A particularly instructive example which resembles the procedure invented 
originally by Migdal and Kadanoff is shown in Fig. 22 [186]. For a better 
understanding of this procedure, we first introduce different graphical rep- 
resentations of the scattering matrices of the saddle points which are more 
convenient for the present purpose (Fig. 21). The original scattering matrices 
S and S' relate outgoing and incoming amplitudes, Vw, VVt an d V'm, Vm) 
respectively. This is represented as in the top part of Fig. 21. Rewriting the 
equations such that the amplitudes on the right hand side of the saddle point 
are obtained as functions of the amplitudes on the left, ip in , ip ont and ipi n , i/j ou t, 
respectively, one obtains the usual transfer matrix representation which we 
have used so far (Fig. 21, right). The saddle point vertex is now represented 
by a horizontal dashed line. For S, the corresponding representation is shown 
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Fig. 21. Representations of saddle point scattering. The scattering matrices S and 
S' relate outgoing with incoming flux amplitudes (top). The same information can 
be expressed by transfer matrices that relate top and bottom amplitudes (left) and 
left and right amplitudes (right). 



in the left part of Fig. 21. With this, the original random network of Fig. 8 
can be graphically represented as shown in Fig. 22 a. For simplicity, the sys- 
tem is assumed to consist of independent saddle points with randomly varying 
saddle point energies. The phases associated with the links between them are 
assumed to be completely random and independent. 

Starting from this a renormalization group transformation may be constructed 
by removing every other saddle point line in a each row and column, and 
replacing the remaining single saddle point lines by two. Thereby, one reaches 
the situation shown in Fig. 22 b. One observes that every other closed loop of 
amplitudes now is completely disconnected from the rest of the system and can 
be removed. On the other hand, each of the remaining saddle points is replaced 
by two saddle points in series. By combining these to new scattering centers 
as indicated in Fig. 22 c one arrives at exactly the same lattice structure as in 
Fig. 22 a but with renormalized saddle point scattering centers (Fig. 22 d). 

By iteratively applying this procedure, more and more of the network is incor- 
porated into the scattering properties of the scattering centers until eventually 
one of the latter incorporates the whole system. If the energy is lower than the 
saddle point energy E*, the transmission probability will eventually renormal- 
ize exponentially to zero. The system becomes completely localized since the 
closed loops labeled with E < E* (Fig. 23) will be completely disconnected. 
If on the other hand the energy is higher than the saddle point energy, the 
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(a) 



(b) 





(c) 



(d) 



Fig. 22. Illustration of the renormalization procedure on a tile lattice, (a) Lattice of 
original saddle points (b) Second-generation saddle point lattice where every other 
saddle point in the rows and columns are removed and the remaining saddle points 
are doubled. Isolated rings of amplitudes can be removed, (c) By renormalizing 
the saddle point as indicated one arrives (d) at a lattice that, apart from a scale 
transformation is structurally identical to the original one [after [186]]. 



transmission probability will renormalize to one. The states are again local- 
ized since now the loops labeled with E > E* in Fig. 23 become disconnected. 
Thus, the transmission probability has two stable fixed points T = and 
T = 1 which correspond to the localized phases. The transmission probabil- 
ity becomes independent of the "size" of the scatterer exactly at the critical 
point E*. Here, we have an unstable non-trivial fixed point T* ^ 0. Even an 
infinitesimally small deviation of the energy from the critical point localizes 
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the system completely. 




Fig. 23. Illustration of localization in the renormalized network for energies below 
the critical point, E < E* , where the transmission renormalizes to zero, and above 
the critical point, E > E* , where the transmission renormalizes to one. Extended 
states can exist only at E = E* where the transmission has a non-trivial fixed point 
T* / 0. 



For the tile lattice, the calculation of the scattering properties of a new scatter- 
ing center by combining two original ones is an easy but nevertheless instruc- 
tive exercise. The transmission matrix of the new scattering center (Fig. 22 c) 
can be determined in terms of the transmission matrices of the original scat- 
tering centers by performing simple matrix multiplications. Consider, for ex- 
ample, the diagram at the bottom of Fig. 22 c. First, one calculates the trans- 
mission of the two pairs of saddle points in series at the top (A) and the 
bottom (B) of the graph. Second, the two resulting effective scatterers are 
combined "in parallel" in order to obtain the transfer of amplitudes between 
channels (ifim, Vwt) and V'out)- Eventually, one finds for the typical total 
transmission probability T := exp (InT) (with T = t 2 and (. . .) the ensemble 
average) the renormalization relation [186] 

T' = 2T 2 — T 4 := f(T; 2) (182) 



where T' is the typical transmission probability of the renormalized scattering 
center. 
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This scaling relation has three fixed points defined by T* = IT* 2 — T* 4 , two sta- 
ble ones at T* = and T* = 1 and an unstable one at T* = (y/E-1) /2 = 0.618. 
The stable fixed points correspond to the localized states at energies away from 
the saddle point energy. The unstable fixed point corresponds to the quantum 
critical point. The localization length exponent is found by linearizing around 
the unstable fixed point 



with A = [df/dT] T * = 6 - 2^5. One gets v « 1.635. 

In the general case, where one saddle point vertex is replaced by b original 
vertices, the result is 

T' = 1 — (1 — T b ) b := f(T; b) . (184) 



Assuming b = 1 + e one obtains the infinitesimal Migdal-Kadanoff transforma- 
tion which has a fixed point at T* = 1/2 and an eigenvalue A = 1 + 2(1 — In 2)e. 
This gives an exponent 

In 6 1 , N 

In A 2(1 -In 2) v ; 

and the beta function 
d 



f(T; 1 + e) = TlnT — (1 — T) ln(l - T) . (186) 

e=0 



Although the result for the critical behavior is far from being satisfactory, it is 
nevertheless remarkable that the very existence of a non-trivial quantum criti- 
cal point is correctly predicted by the model even in the crudest approximation 
for the renormalization group transformation. It is therefore worthwhile to in- 
vestigate whether one can obtain more reliable quantitative results by refining 
the approximations. 



6.2 State-of-the-art results for the hierarchical lattice 



A different procedure is outlined pictorially in Fig. 24. This is the model of 
a hierarchical lattice. In the first generation, the system consists again of the 
original, statistically independent saddle points. In a first step, a certain num- 
ber of these original saddle points (five in the example of Fig. 24) are combined 
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to form a new scattering center. The corresponding scattering matrix is cal- 
culated approximately as a function of the original ones. This is repeated: the 
new scattering centers are combined again into new units (Fig. 24 c). Their 
scattering properties are calculated as functions of the previous ones, and so 
on and so forth. 

The crucial approximation of the procedure is that in each step only two 
incoming and outgoing channels are taken into account (Fig. 24 b) such that 
the new unit can be considered as a new, second-generation saddle point. 
The exact scattering matrix of the units in each generation would of course 
contain much more channels, and when repeating the construction, the number 
of channels would explode. Neglecting all of these channels apart from four is 
in fact the most severe approximation of the method — as also made in the 
previous procedure — since it cannot be very well controlled. 

However, instead of viewing the hierarchical lattice as an approximation to 
the random network model one may also consider it as a model in its own 
right for which the critical behavior may be determined exactly. This would 
yield only an approximation to the critical behavior of the network model, but 
if the critical behavior was universal — independent of the microscopic details 
- the exponent of this model should be the same as for any other model in 
the same universality class. We will come back to this point below in more 
detail. 

For the derivation of the renormalization group equations of the hierarchical 
lattice, we closely follow [33]. We consider the five saddle points in Fig. 25. 
The amplitudes must satisfy the five relations 

M-MM • (-LV...5) (187) 

where i and o denote the input and output channels, respectively, and and 
U, (= \J\ — rf) (i = 1, . . . , 5) are real quantities. The first goal is to replace 
the five saddle points by one. 

Some of the amplitudes differ only by a phase factor. For example, 

V>2,i = e^Vi,o, (188) 



where #2,1 is the random phase factor gained when traveling from the saddle 
point 1 to 2. The phase factors are related to the Aharonov-Bohm fluxes 
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(a) 



(b) 



(c) 



Fig. 24. Illustration of the renormalization procedure on a hierarchical lattice con- 
structed iteratively. (a) Saddle point represented by a dot with two incoming and 
outgoing channels, (b) Second-generation saddle point consisting of five original 
(first-generation) saddle points. A total of eight incoming and outgoing channels 
are here neglected, (c) Second level of iteration where five second-generation sad- 
dle points are combined to form a third-generation saddle point consisting of 25 
first-generation saddle points. In contrast to the original network (Fig 8), the hi- 
erarchical network contains considerably fewer saddle points and links [after Ref. 
[186]]. 

4>i, {i = 1, 2, 3, 4) obtained when traveling around a closed loop via 

#5,1 + #3,5 + #1,3 
#3,2 + #4,3 + #2,4 
#2,1 + #3,2 + #1,3 
#3,4 + #5,3 + #4,5 



n 
h 
h 



(189) 
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Since the phases 6 are assumed to be randomly and uniformly distributed 
between [0,2n), so are the fluxes 4>i- 




Fig. 25. Illustration of a transformation that can be used in a renormalization 
group transformation. Five original (first-generation) saddle points indicated here 
by dashed lines (right) are renormalized into one single second-generation saddle 
point. This relates (^2,^5) to (^1,^(4). The quantities (pi... ^4 are random 
independent Aharonov-Bohm-fluxes that penetrate closed loops as indicated. 



/,From the relations (187), one can derive that the total outgoing amplitudes 
(■0o,2, ^0,5) are related to the incoming amplitudes (Vv, $,4) according to for- 
mally the same relation as Eq. (187) 





(190) 



The new reflection and transmission coefficients, f, t, characterize the scatter- 
ing properties of the renormalized "super" saddle point. The new transmis- 
sion and reflection amplitudes can be obtained straightforwardly by solving 
Eq. (187). One obtains after some tedious algebra 



t = 



Ms + t 2 t 4 e^ 3+< W -hhUe^ -t 1 t 3 t 4 e^ 4 



D 



+ 



tir 2 r 3 r 4 t 5 e^ 2 -rit 2 r 3 t 4 r 5 e^***-* 1 ) 



D 



(191) 
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The corresponding reflection amplitude is obtained from r = \Jl — t 2 

r^r 5 — r\r 2 e l ^ 1+ ^ — n^^e^ 1 +r 2 r 3 r 5 e^ 2 



r = 



D 

tit 2 t 3 r 4 r 5 e^ 3 +r 1 r 2 t 3 t 4 t 5 &fa+<h-<l>*) 



D 

where D is the abbreviation for 



(192) 



D — l — nr 3 r 5 e^ 1 -ht 2 t 3 ^ 3 ~t 3 Uh e i4,A +r 2 r 3 r 4 e^ 2 

+ ht 2 UU e i( ^ 3+04) -nr 2 r A r 5 e^ 1+<fe) . (193) 

The transformation (191) allows one to generate the new probability distribu- 
tion of the transmission coefficient P(t) from the distribution P(t). A certain 
distribution is unchanged by this transformation, which corresponds to the 
fixed point distribution. Slight deviations from the fixed point distribution 
increase after the renormalization transformation. This can be used to define 
the critical exponent v. 

The renormalization group transformation of the distribution of the trans- 
mission amplitudes has been determined numerically by using this approach 
[34,188]. From P(t) the distribution of the conductance G = t 2 can be deter- 
mined, 

P(G) = lp(f) . (194) 

This may be transformed to the distribution Q(z) of the heights of the saddle 
points z measured relative to the critical energy e = (cf. Eqs. (75), (84)). 
This can be obtained from Eq. (194) by using the relation 

G = — *— -. (195) 
e z + 1 v ' 



The result is 

dG 



Q(z) = P(G) 



dz 



4cosh / (^/2) Ve 2 + 1 



(196) 



The distributions P(G) and Q(z) are shown in Fig. 26. 

The scaling properties of the distribution Q(z) may be used to determine the 
critical exponent with very high accuracy. The idea is that an initial shift of 
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G 




Fig. 26. Top: The fixed point distribution (o) of the conductance, P(G), at the quan- 
tum Hall transition obtained by starting from different initial distributions denoted 
by squares, diamonds, triangles, and the distribution after the 16th renormaliza- 
tion step (+). Bottom: the corresponding distributions of the heights of the saddle 
points, Q(z) [From [188]]. 

the distribution by an amount z will be amplified during the renormalization 
procedure since the quantum Hall fixed point is unstable. After n steps, the 
maximum of the distribution will be shifted by an amount z max ^ n = \ n zo. After 
a certain number of steps, til, the saddle point will be no longer transparent. 
This limit can be denned by 

W,n L = \ nL z « 1 . (197) 

By denning the localization length £ with the length 2 nL a (a lattice constant 
of the original lattice) one finds that £ diverges for z — > 0, 

£ oc az u (198) 
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with the exponent v = In 2/ In A (cf. Eq. (183)). The critical exponent can 
then be determined from 

ln2 n , . 

(199) 



m (z max ,n)/ z 



with z mSuXtn /zQ = A™. Figure 27 shows results of the numerical calculation 
[34,188]. The numerical analysis of the data yields the result v = 2.39 ± 
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Fig. 27. The critical exponent z/ determined from the scaling behavior of the distri- 
bution of the heights of the saddle points Q(z) with the effective system size 2 n in 
the nth renormalization step. Dashed line: v = 2.39. Inset: maximum z max ,n as a 
function of zq. The slopes yield A™ and thereby the critical exponents v{n). 

0.01, in excellent agreement with previous direct numerical simulations for 
the localization length [34,188]. As mentioned above, by assuming the model 
as an independent realization of the fixed point ensemble, one can interpret 
this result as an independent corroboration of the universality hypothesis for 
the quantum Hall phase transition. 

Though the value of the critical exponent is extremely precise, the form of 
the critical distribution of the two terminal conductance is qualitatively dif- 
ferent from the one of the Chalker-Coddington random network. Instead of 
P(G) — > in the limit G — > (Fig. 15), the renormalization group approach 
indicates divergence P(G) — > oo. Also near G — 1, the present approach sug- 
gests P(G) -> oo, while the actual P(G) is constant. For G > 1, P(G) = 
in the case of renormalization group, but in the simulation for the original 
model P(G) shows a tail in this region. This indicates that the renormaliza- 
tion group approach in fact does not approximate the Chalker-Coddington 
network model although the underlying truncated-network model belongs to 
the same universality class. 
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This renormalization group has recently been extended to the calculation of 
the energy level statistics [35]. The finite size scaling analysis of the energy 
level spacing, Eq. (139), has been applied to estimate v to be 2.37 ± 0.02. 



7 Hamiltonians Related to the Network Model 



The model defined above is not only suitable for quantitative studies of the 
critical behavior. By constructing effective Hamiltonians starting from the 
transfer matrix, one can recognize that it is of much more general importance. 
In the present chapter we will achieve such a construction. We first introduce 
an equivalent tight binding Hamiltonian. In effective mass approximation, this 
will be shown to be equivalent to a Dirac Hamiltonian. For the latter, some 
analytical results are summarized which put the quantum Hall phenomenon 
in a new perspective. 

In this and in the subsequent chapters, we will use units such that h — 1, 
m* = 1, and e = 1, in order to simplify notation. However, in results for 
the conductivity tensors, the correct unit e 2 /h is reintroduced since it is of 
physical importance. 

The above network system has been completely specified in terms of its trans- 
mission or scattering properties at a given energy of the original model of an 
electron moving in a strong magnetic field and a slowly varying random po- 
tential. By using this as a starting point for defining a Hamiltonian, we fix an 
entire spectrum of (quasi-) energies which contains as a parameter the original 
energy value. It is not obvious whether or not there is some arbitrariness in 
such a procedure [189]. 



7.1 A Tight Binding Hamiltonian 

We first want to establish a connection with a nearest neighbor tight binding 
Hamiltonian. We start from the unitary matrices S that describe the scattering 
at the individual saddle points, Eq. (87). We first rewrite S by redefining 
phases, 4>i ~ * 4>i + 71 an d 04 ~ * 04 + n i again in order to avoid unnecessary 
notation. This does not change the physics, as it will be seen below that one 
of the relevant system parameters is the total phase accumulated around a 
unit cell of the lattice, = J2j <f>j, which is the number of flux quanta in the 
unit cell multiplied by 2n. 
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(200) 



Then, we rearrange the network in the way shown in Fig. 28. The loops that 
are connected by the saddle points, described by their scattering matrix So, 
are arranged in the form of a two dimensional square lattice such that the 
centers of the loops are associated with the lattice points 



R j xy ^^x ~t~ IJ^y 



(201) 



with x, y integer numbers. This notation might appear unconventional, but 
is convenient for emphasizing the connection with the real space. The x- and 
^/-directions are assumed to point into the directions of the links that connect 
S and S'. The four links between adjacent saddle points within each of the 
loops (denoted by arrowheads labeled with 1,2,3,4 in Fig. 28 [189]) arc now 
associated with four site states within a unit cell of the lattice. The site states 
are characterized by a lattice vector and a quantum number A in the lattice 
cell, and they are assumed to form a complete set, 

4 

]T \x,y,\){x,y,\\ = l. (202) 

x,y,\=l 



These site states are connected via the tunneling and reflection matrix ele- 
ments through the saddle points within a given unit cell, and between nearest 
neighbor cells. 

In order to determine the matrix elements of the effective Hamiltonian, we 
adopt the method used in [121] and [189]. We write for the vector of the 
amplitudes in the lattice cell at R = xe x + ye y after M + 1 iteration steps 

4 

^ RA (M+1)= ]T CWa'Vva'M- (203) 

R'\' = l 



The unitary operator U can be interpreted as describing the evolution of 
the wave function between steps M and M + 1. The eigenstates of the 4L- 
dimensional matrix U (L number of sites) 

U | ip a ) = e iuJa{E) \ip a ) (a = l AL) (204) 
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Fig. 28. The network model rearranged as a tight binding model on a two dimen- 
sional square lattice with four states per site. The links between saddle points are 
associated with the four directed site states (arrowheads) in the unit cell. These site 
states are connected by tunneling through the saddle points represented by <C>. 

with eigenvalues 1 (u a (E) = 0) are the stationary states of the network with 
an energy parameterized by the energy dependent reflection parameter of the 
saddle points, r := r(E) = Jl — t 2 {E). 

Formally, one can interpret the above Eq. (203) as a "time dependent" Schro- 
dinger equation by writing 

4 

V>*a(M + 1)-V*a(M)= J2 (u-i) rX,r'X' 

Vvy(M). (205) 

R',A'=1 



This suggests 



U = 1 - -H. 



(206) 



As the Hamiltonian H must be self adjoint we define 

H = i (U - Ut) . (207) 



Note that the Hamiltonian is here dimensionless, since it contains a factor 
At/h, with At = 1 the width of the time step in Eq. (205). 
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The Hamiltonian has quasi-energy eigenvalues 



sin u a . 



(208) 



Since each lattice cell is connected via four saddle points to its nearest neigh- 
bors, U must have a 4 x 4 block structure. The midpoints of the links form 
a bipartite lattice. This means that each of the 4x4 blocks can consist only 
of two non-zero 2x2 blocks. By arranging the amplitudes according to the 
sequence (ip R i,ip R 3,ip R 2,ip R 4) T one gets for the 4x4 block connecting with 
each other the site states within a given lattice cell at R, and with the states 
in the nearest-neighbor cells, 



U 




(209) 



such that the 2x2 block M connects the states vp R i,tp R 3 with ip R 2,^ R 4 and 
the block N links the states ip R 2,i> R 4 with ip R i,ip R 3- 



This gives for the 4x4 block of the effective Hamiltonian 



H = 

2i 



( Nt-M\ 
Mt -N 



yU\ . 



(210) 



The matrices 

/ te i<p lT x T v re i0! N 

re ^3 _ te i<t>3 T X T V_ 



M 



N 



re i<t> 2 te i4,2 T x T y + 



(211) 



(212) 



contain the translation operators r± that connect nearest neighbor cells, 

t±^ha = *p R ± ex \ rli) RX = ip R ± ey x ■ (213) 



The Hamiltonian is then given by 



H, = — 



2i 



_ re -i<t>A + te i<i> 3T x T y 



(214) 
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This means that within a given cell, say at R, the tunneling matrix elements 
are proportional to the reflection amplitudes of the saddle points r 



h R i jR 4 = 


2i 


h R 2,Ri = 


r 




h R 3, R 2 = 




y 

2i 


h R 4, R 3 — 


r 


~2i 



(215) 
(216) 
(217) 
(218) 



On the other hand, the eight coupling matrix elements to the nearest neighbor 
cells are proportional to the transmission amplitude t. They are 



h 



h — * J4>2 

ri R 2,[R+(i,i)]3- — e ' 



H3,[H+(1,-1)]4 



2i 



h R 4,[ R+ (-i -i)]i = — e i<t>4 , 

h R i,[ R+ ^i,i)]2 = — el<t>1 > 
and the four remaining matrix elements are the conjugates of these. 



(219) 
(220) 
(221) 
(222) 



It is useful to consider as a starting point for the discussion of the properties of 
this Hamiltonian the periodic limit. This will also yield the Dirac Hamiltonian 
in the limit of small wave number, i.e. in effective mass approximation. 

The periodic case is achieved by assuming the phases in each of the unit cells 
and all of the saddle points as identical. The Hamiltonian can be diagonalized 
exactly with the Bloch Ansatz 



MR) = e lq - R u x (q) . 
The resulting band structure is (Fig. 29) 



(223) 



£,-(*!) = (-l)V2 



1 — 2rt cos <p cos (q x — A x ) cos (q y — A y ) 



n V2 



± sin a/1 — 4rH 2 cos 2 (q x — A x ) cos 2 (q y — A y 



(224) 



with j = 1, 2, A x = (0x - 3 )/2, A y = (04 - 2 )/2 and = El=i 0a the flux 
through the cell multiplied by e/H. Exactly at the energy of the saddle points, 
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r = t = 1/ \/2, there is no gap between the bands. For r^ia gap A = ei — e 2 
is opened. For = it is 



A = 2v / 2^( 1 ~ 2rt ) ~ m ( 225 ) 



for r = l/y/2 + m/A (m < 1). 




(a) (b) 

Fig. 29. The band structure of quasi-energies ej (q — A) of the Hamiltonian repre- 
senting a network with identical saddle points and equal phases in all of the cells of 
the tight binding lattice such that there is half of a flux quantum in each unit cell 
{4> = tt). In this case, the four bands (Eq. (224)) are degenerated with respect to ±. 
(a) The spectrum exactly at the saddle point energy, r = t = 1/V2, (b) for r = 0.5. 



The formation of the gap can be discussed in better detail by expanding the 
band structure near q — A « 0. On finds for small m 

e ^ q ) = m 2 + (q-A) 2 . (226) 



Thus, the existence of the gap is closely related to the deviation from the 
saddle point. 

In the following, we discuss briefly the disordered version of the network. 

If the reflection amplitude is unity (r = 1, t — 0), i.e. if the energy E is well 
below the saddle point energy, the unit cells decouple. The system consists 
then of independent wave functions localized within the grey unit cells in 
Fig. 28. On the other hand, if t = 1 (r = 0), i.e. for an energy well above the 
saddle point, the independent cell states are localized in the complementary 
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lattice of the white unit cells. The quasi-energy eigenvalues in these decoupled 
limits are 



ef(R,ct>) = (-iyV2 



l T -L(l-cos0) 1/2 



1/2 



0" = 1,2) (227) 



and (j) = (f>(R) is the flux through the loop at R. For completely independent 
and random phases, the fluxes in the individual cells, and thus the quasi- 
energies e a are statistically independent and random. When varying the en- 
ergy E, i.e. varying r(E) close to the value of 1, the quasi-energies e(E) will 
perform random walks in the (E, e)-space. The energies of the original prob- 
lem, Ep(e a = 0), will then be statistically independent and random. They 
correspond to the localized states associated with the valleys and hills in the 
potential landscape. 

If we assume that all saddle points have the same parameter r and the phases 
are independently and randomly distributed one can see that the localization 
lengths of the eigenstates of H are finite and independent of e a . As r — ■> r c = 
t c — 1/ y/2, the localization length diverges uniformly throughout the spectrum 
[189]. This reflects the singularity of the localization length as a function of 
the energy E of the original problem. 

We note that with disorder, i.e. assuming the phases to be independently 
and randomly distributed, the matrix elements of the Hamiltonian remain 
nevertheless strongly correlated. These correlations introduce peculiar features 
into the localization behavior of the corresponding eigenstates. In particular, 
they are responsible for the localization-delocalization transitions occurring in 
the sub-bands induced by the magnetic field. A very similar observation has 
been made earlier for the random matrix model where the matrix elements of 
the banded matrix are also strongly correlated [190] . 



7.2 The Dirac Hamiltonian 



One can also relate the network model with the two dimensional Dirac equa- 
tion. This can be done most straightforwardly by considering the effective 
mass approximation of the tight binding Hamiltonian. Alternatively, one can 
consider the two-step unitary operator [190] 



U 2 



I MN ^ 



NM 




(228) 



83 



The system of four equations decouples into pairs, and one can deal with, say, 
only the upper block V. 

Furthermore, by assuming that the displacement operators act on smooth 
functions one can replace 



T%=\±d x , T V ± = l±dy 



(229) 



For small phases and with the definition r = r c + m/4 (t ~ t c — m/4) one can 
expand 



V = e 



-iH 



1 -iH 



(230) 



with the unit matrix 1. This is justified for small H. From the matrix elements 
of M and N, Eqs. (211) and (212), one finds 



V « 1 + 



— d x + %A X d y — iA y + m ^ 



y d y - iA y - m d x - iA x 



i<j>l = 1 - iH. 



(231) 



With the unitary transformation 



R 



1 i -1 



V2 



i 1 



(232) 



one can transform 
H = RHR 1 



(233) 



such that one obtains finally the two dimensional Dirac Hamiltonian 

H = (PX ~ A X )0 X + (Py - Ay)Oy + mo \ + <j)l . (234) 



Here, pj = —idj (j = x, y) are the components of the momentum operator 
and 



1 

1 



-i 

1 



1 
-1 



(235) 



are the Pauli matrices. In this Hamiltonian, randomness can be introduced 
in different ways. Via randomness in the individual phases one can make the 
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components of the vector potential random. Randomizing the total Aharonov- 
Bohm phases in the loops produces randomness in the scalar potential <j). 
Finally, assuming the tunneling parameters of the saddle points to be random 
gives fluctuations in the mass parameter m. 

Introducing disorder in the various parameters of the model leads to breaking 
of symmetries, as noted in [120] and indicated in Tab. 2. The first term in the 
Hamiltonian, H = cr ■ p [cr = (a x ,a y ), p = — iV = —i(d x ,d y )], is invariant 
under an effective time reversal operation 

H = <J z H* a z . (236) 



This means that there is a Kramers degeneracy, ip and cr z ip* are degenerate 
eigenstates as is easily seen. Similarly, it is easy to see that ip and a z ip are 
eigenstates to energies E and —E, respectively. This is equivalent to particle- 
hole symmetry 

— H = a z ii a z . (237) 



Finally, H is parity invariant under x and y-reflections PjjPj = —j, (j = x,y), 
for instance 

Hq = P y o~ x HoO- x P y . (238) 



The Dirac mass term H m = ma z possesses none of these symmetries, as is eas- 
ily seen. The vector potential term H A = —A ■ cr preserves only particle-hole 
symmetry, and the random scalar potential = 01 only time reversal sym- 
metry. These symmetries are very important for understanding the quantum 
phases of the system as will be seen below in more detail. 



7. 3 Some Results for the Two Dimensional Dirac Model 



In this section we discuss briefly some instructive analytical results obtained 
for the Dirac Hamiltonian. 

The correspondences between the quantum Hall problem, certain tight binding 
models and the two dimensional Dirac model have been noticed by several 
authors [120,191,192,193]. Fisher and Fradkin [191] have obtained the Dirac 
model by starting from a two dimensional tight binding model with diagonal 
on-site disorder in a perpendicular magnetic field with half a magnetic flux 
quantum in the unit cell. They constructed a field theory for the diffusive 
modes which was shown to be in the same universality class as the orthogonal 
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symmetry 


cr • p 


ma z 


<t> 


a ■ A 


parity 


yes 


no 


no 


no 


time reversal 


yes 


no 


yes 


no 


particle-hole 


yes 


no 


no 


yes 



Table 2 

Symmetry breaking by disorder in the various parameters of the Dirac model (p 
momentum operator; m Dirac mass parameter; cr vector of the Pauli matrices a x , 
a y ; <j) scalar potential; A vector potential). 



0(2n, 2n)/0(2n) x 0(2n)(n — > 0) non-linear a-model. This implies that all 
states are localized as in the absence of a magnetic field. It suggests that 
if derealization occurs with magnetic field, it must be a direct consequence 
of breaking time-reversal symmetry instead of some other properties of the 
field. Generalizations to the several- channel scattering problem have also been 
discussed [194]. 

Ludwig and collaborators [120] have used a tight binding model on a square 
lattice with nearest and next-nearest neighbor coupling, half a flux quantum 
per unit cell and a staggered potential energy //(— l) x+ f as a starting point. 
At low energies, the model was shown to be equivalent to a Dirac model with 
two Dirac fields. They found that without disorder this model exhibits an 
integer quantum Hall phase transition as a function of some control parameter 
which is essentially the mass m of the lighter Dirac field. This is similar to the 
transitions of the Hall conductance between integer multiples of e 2 /h obtained 
in the original clean Landau model as a function of the Fermi energy, and 
does not mean that a plateau exists when the electron density is varied. The 
transition of the Dirac model has been shown to belong to the two dimensional 
Ising universality class. The associated exponents and the critical transport 
properties were determined. 

The density of states, 

p(£) = M e (|£|- m ), (239) 

vanishes at E = 0. It can readily be obtained from Eq. (226) with A = using 



p{E) = I d2q I 6 ( E ~ + m2 ) + 5 ( E + v^ 2 + m2 ) 



(240) 
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Applying linear response theory to the Dirac system one can determine the 
Hall conductivity by calculating the ratio of, say, the current density in the 
x-direction, j x and the electric field, E y , in the y-direction 

]^ = e^m_ t d 2 qd^ 
xy E y h 2tt 2 J [{iuj - E) 2 -q 2 - m 2 ] 2 ' 1 ' 



It is found that at zero energy, where the above density of states vanishes, the 
Hall conductance jumps by e 2 /h at m = (Fig. 30) 

*.<".) = ^ (242) 



The heavier Dirac field contributes towards the Hall conductance with e 2 /2h 
such that the total Hall conductance jumps from to e 2 /h. Simultaneously, 
the magneto-conductivity a xx is non-zero only at this critical point, 

e 2 



with the Kronecker-symbol 5 m .o equal 1 for m = and for m ^ 0, and 
the constant o"o is 7r/8. The critical point shows time-reversal, particle-hole 
and parity invariance. Thus, the clean two dimensional Dirac model exhibits 
a quantum Hall transition at E — m — 0, i.e. a step in the Hall conductivity. 



m 




E 



Fig. 30. The phase diagram of the ordered Dirac model [120] showing regions 
of non-zero density of states p (\E\ > \m\) and non-zero Hall conductance 

&xy — 

e 2 /2h (\E\ < \m\). 



According to our above derivation of the Hamiltonian, introducing spatial 
randomness in the Dirac mass, m = m(x,y), is equivalent to introducing 
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randomness in the saddle point energies. This breaks all of the above symme- 
tries (cf. Tab. 2). However, this kind of randomness alone, if it is sufficiently 
weak, does not introduce a nonvanishing density of states at zero energy. This 
means that it cannot place the system into the quantum Hall universality class, 
in contrast to what one might suspect. For sufficiently slowly varying mass, 
the zero-energy wave functions are argued to be confined along the contours 
m(x,y) = since for any region with a non- vanishing Dirac mass the energy 
must be non-zero. When the mass varies slowly in space, the transition can 
be interpreted in terms of the percolation of these states. In the absence of 
randomness of the phases (i.e. the vector potential), the corresponding critical 
exponent of the correlation length is that of the classical percolation model. 
The analysis in Ref. [120] for weak randomness using the replica trick and per- 
turbative renormalization group analysis strongly suggests that the random 
Dirac mass term is always marginal for the critical properties of the quantum 
Hall transition. The system scales to vanishing disorder, and the transition is 
described by the free Dirac theory. The replicated effective action is formally 
equivalent to that of the random-bond Ising model (cf. next chapter) with the 
randomness in the Dirac mass corresponding to the randomness in the Ising 
bonds. 

Randomness in the scalar potential is equivalent to assuming random fluxes 
piercing the plaquettes of the network. For small random scalar potential one 
can use the perturbative results obtained for the random Dirac mass. This can 
be suspected from the form of the Hamiltonian 

H = a p + ma z + 01 . (244) 



Specifically, it can be shown that the problem with real disorder in the Dirac 
mass is formally equivalent to a problem with a purely imaginary random 
scalar field which corresponds to a negative disorder strength, (0 2 ) < 0. 
Equivalently, a positive scalar disorder strength maps onto a negative dis- 
order strength for the Dirac mass. Since randomness in the Dirac mass was 
found to be marginal, the opposite is true for the scalar randomness, and the 
renormalization flow of the scalar disorder will now be away from the fixed 
point instead of towards it. Thus, the scalar randomness drives the system 
to some strong coupling regime. Eventually, one then expects the generation 
of a non-vanishing density of states. It is further argued by considering two- 
particle properties that the transition very probably has to be described by a 
symplectic non-linear sigma model due to the time reversal invariance of every 
member of the ensemble [120]. 

The case of only a random vector potential was found to be particularly inter- 
esting since it can be treated to a large extent analytically [120,195,196,197,198,199,200]. 
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This limit has several remarkable properties. The Hamiltonian is 



H A = cr ■ p + <r • A (245) 



and according to what has been said before, the random vector potential is 
due to to the random phases along the links of the network. We assume that A 
satisfies a Gaussian white noise distribution with zero mean and variance A^. 
The corresponding random magnetic field, B = V x A, is then also Gaussian 
distributed with the variance 

(B(k)B(k')) = (2rr) 2 k 2 A A 5(k + k') (246) 



with the Fourier transform B{k) of the magnetic field [perpendicular to the 
(x, y)-plane] . As the Dirac mass is assumed to vanish, this Hamiltonian is 
expected to implement a critical theory. In fact, the model exhibits a fixed 
line corresponding to multi-fractal wave functions at E = [120]. 

For nonvanishing uniform Dirac mass the squared of the Hamiltonian is 

H 2 = H 2 A + m 2 (247) 



for each realization of the random vector potential. Thus, the energy eigen- 
values satisfy E 2 > m 2 . The gap in the pure Dirac model with m ^ is thus 
not closed by the randomness in the vector potential, at least as long as the 
disorder is small. 

For this model, several single particle and two-particle properties can be de- 
termined exactly. As an example, we consider the zero-energy wave function. 
It is useful to use the Coulomb gauge, V • A = 0. This implies that A can be 
written in terms of a scalar field $(:r, y) 

A x = d y $ (248) 
A y = -d x $ . (249) 

By inserting the real Ansatz wave functions 

tt ± oc(l±<7 z ) | | (250) 



into the Schrodinger equation of the Hamiltonian Eq. (245) one easily verifies 
that they are exact nodeless eigenfunctions corresponding to E = 0. Fur- 
thermore, it follows from the symmetry properties of the Hamiltonian and the 
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assumption of zero fluxes in the plaquettes that these are all of the zero-energy 
eigenfunctions. 



These random wave functions are not localized, and not normalizable in the 
thermodynamic limit. In order to quantify their statistical nature it is useful 
to introduce their counterparts normalized to a square of the size L 2 . 



C 



(251) 



with 



C 2 = Je 2 ^dxdy, 

L 2 



(252) 



where the integration is over the square L 2 and to consider the moments of the 
corresponding density, the average inverse participation numbers introduced 
in section 4.3 



Pq (L) = (IVfV 



(253) 



with (. . .) L 2 denoting the configurational average in L 2 . 

For the above normalized zero energy wave function, Eq. (251), with Gaussian 
distributed $, the average inverse participation numbers can be determined 
[120,201,202,203,204]. One finds for q c = J2ir/A A > 1 



r(q) 



2(q-l)(l-q/q 2 ), 



2q(l - sgn q/q c ) , 
For q c < 1 one gets 



\q\ < Qc 

\q\ > Qc, 



(254) 



r(Q) 



-2(1 -q/q c ) 2 , \q\<q c , 

4 (g-|g|)g^ 1 , \q\>q c , 



(255) 



The nonlinearity in q indicates multi-fractal scaling of the extended wave func- 
tion. The resulting /(a)-spectrum is parabolic with 



a 



2(1 + q~ 2 ), q c >l, 



(256) 
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depending on the variance of the random vector potential A a, Eq. (246) (see 
Sect. 4.3). 

The density of states was found to show a power law dependence of the energy- 
near \E\ pa with the exponent varying continuously upon moving along the 
fixed line [205], 

p(E) oc E^/ z , (257) 



with 



l + 2/g r 2 , q c >l, 

/yc ' y ' (258) 

^ 4/q c - 1, q c <l. 

For small disorder, the density of states vanishes. For A^ = ir, p(E — > 0) = 
const, and it diverges for A^ > n. 

The diagonal conductance was found to be e 2 /nh along the fixed line. 

The Dirac model with only randomness in the vector potential belongs to the 
AIII symmetry class to be discussed in the chapter 9 (see also Tab. 5). The 
critical wave functions at E — for random it phase gauge field [206,207,208] 
as well as for arbitrary random gauge field [209] have been investigated in 
detail. 

If the phases along the links of the network are assumed to be independent 
and random, the equivalent Dirac model has randomness both in the vector 
potential and in the fluxes through the plaquettes. Such a system has none 
of the above symmetries (cf. Tab. 2) and obviously belongs to the original 
Chalker-Coddington class with a fixed point that corresponds to the quantum 
Hall phase transition. 

The scaling picture of the Dirac model with all of the three different kinds of 
disorder — random Dirac mass (parameter A m , random saddle point energies), 
random scalar potential (parameter A^, random fluxes through the plaquettes 
of the network) and random vector potential (parameter A^, random link 
phases) — has been qualitatively sketched by Ludwig and collaborators on 
the basis of their analytical results [120] (Fig. 31). The phase diagram consists 
of a critical line (A m = A^ = 0, A^ ^ 0) which is unstable with respect to both 
A m and A^, a two dimensional Ising fixed point at A m = A^ = A^ = which 
is only stable for A m ^ 0, and a fixed point at some A^ ^ which was argued 
to belong to the universality class of the symplectic non-linear sigma model. 
According to this scenario, the fixed point of the genuine integer quantum Hall 
transition is in a strong coupling regime with all of the three different kinds 



91 



of disorder present. Until now, this has not yet been accessible analytically. It 
will be the subject of subsequent chapters to review the approaches to extract 
nevertheless analytical information on the GIQHE, and to shed some light on 
the obstacles which have still prevented an analytical derivation of its critical 
parameters. 



4a 




Fig. 31. The global phase diagram of the Dirac model including randomness in 
the vector potential, A a, the Dirac mass, A m , and the scalar potential, A^ (after 
[120]). The generic fixed point of the Integer Quantum Hall Effect corresponds to all 
of the three different kinds of disorder nonvanishing, and seems to be inaccessible 
analytically. 



8 Relation to the Two-Dimensional Random Bond Ising Model 



In the last chapter we have seen how mapping to different Hamiltonians can 
help to provide new insights into the properties of the quantum Hall phase 
transition. In this chapter we will address in addition the opposite question: 
given we know the network model and know how to use it for understanding 
the quantum Hall critical scenario, can we use this knowledge for understand- 
ing better the nature of phase transitions in two dimensions? As the random 
two dimensional Ising model for a system of interacting spins is a funda- 
mental prototype for phase transitions in two dimensions, for answering this 
question it is useful first to establish a relation between the random network 
and the Ising model. Such a mapping has been performed by various authors 
[210,211,212]. We describe here the elementary route that has been worked 
out in the seminal paper by Merz and Chalker [213]. 
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8. 1 The Ising Model 



The Ising model in two dimensions is denned by the Hamiltonian [214] 




(259) 



The exchange integrals connect the sites % and j of a regular square lattice. 
The variables Si are the z-components of the spin operators associated with 
the lattice sites. 

We consider the simplest case, where connect nearest neighbors only, and 
spin 1/2. In the ordered limit all of the exchange integrals are the same. This 
has been solved exactly by Onsager in classical paper using Lie algebras [215]. 
A comprehensive treatment using the transfer matrix method is described in 
[216]. The basic reason for the Ising model to be an exactly solvable many-body 
problem is that eventually it amounts to nothing else but the diagonalization 
of a quadratic form. 

However, in the disordered case, when the exchange parameters are chosen 
at random, the model cannot be treated exactly, in spite of being quadratic. 
We consider as independent random variables with a distribution P(Jjj). 
Specifically one can assume a two-component distribution, 



This is the Random Bond Ising Model. For J > 0, it corresponds to a lat- 
tice containing p antiferromagnetic and (1 — p) ferromagnetic bonds. In this 
simple case, the model shows a phase transition between an ordered ferro- 
magnetic phase at low temperatures and a paramagnetic high temperature 
phase for small p. At larger p the ferromagnetic phase is destroyed in favor of 
a spin glass phase in high dimensions, d > 3, with a multi-critical point, also 
called Nishimori point [217,218,219], where the three phases coexist (Fig. 32) 
[220,221,222,223,224,225]. In two dimensions the spin glass phase becomes 
unstable for T ^ 0. However, the multi-critical point survives. Along the 
Nishimori line which crosses the phase boundary at the Nishimori point, the 
internal energy can be calculated exactly. 

The disordered Ising system has been studied using Monte Carlo simulations 
and the transfer matrix method in the spin basis [223,224,225,226,227,228,229,230,231]. 
Mapping the random Ising model to fermionic models can have technical ad- 
vantages: one can avoid the random-sampling errors of the Monte Carlo tech- 
nique and it is possible to avoid the exponential growth of the transfer matrix. 
Pioneering work in this direction has been done by starting from the solution of 



P(J ij )=p6(J ij + J) + (l-p)6(J i 



J). 



(260) 
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Fig. 32. Schematic phase diagram of the Random Bond Ising Model. For small 
fraction of antiferromagnetic bonds one observes a transition from a ferromagnetic 
to a paramagnetic phase with increasing temperature. For larger fraction of an- 
tiferromagnetic bonds, the ferromagnetic phase is changed to a spin glass. At the 
multicritical point M the three phases coexist (dashed: Nishimori line; arrows: renor- 
malization flow). 

the two dimensional Ising model using a Pfaffian [232,233,234,235,236,237,238]. 
In this approach the statistical properties of the system are written in terms 
of the spectral properties of the corresponding matrix. The latter is essen- 
tially a tight binding Hamiltonian on the underlying Ising lattice with random 
hopping matrix elements. Thus, a link between the random bond Ising sys- 
tem and the non-interacting localization problem is established. Alternative 
approaches have been formulated by using Dirac fermions which eventually 
have been condensed into a random network model with a special symmetry 
[120,210,211,212,213,239]. 

The mapping of the Ising model to a random network model is done in three 
steps. First, one introduces the conventional description of the partition func- 
tion in terms of transfer matrices [216]. In the second step, the transfer ma- 
trices are written in terms of fermion operators instead of spinors. Finally, 
the network model is introduced by using the equivalence of second- and first- 
quantized forms of linear transformations. In the next sections, we establish 
the random network version of the Ising model following [213]. 

8.2 Transfer Matrix Formulation of the Ising Model 

We consider the two dimensional Ising model on a square lattice of length 
L and width M with random nearest-neighbor exchange couplings J$(l,m) 
(Fig. 33). The pair of integers (l,m) denotes the coordinates of a lattice point 
in L- and M-directions, respectively. The index 5 = A,/i indicates whether 
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Js{l,rn) couples to a nearest neighbor in the directions along L or M, respec- 
tively. 

1 ^ 1^ 1^ 1^ i i 



-4 



-4 — * 



4) + i( |m ) 

|j_ ld (lm)^J +1( j(lm)^ 
J -i('m) 



1-1 



1+1 



-4 t t f it it- 



m+1 



it- 



m 



| f-m-1„ 



M 



Fig. 33. The two dimensional Ising model on a square lattice of length L and width 
M with exchange matrix elements Js(n,m) between nearest neighbors. 



In order to introduce the transfer matrix description we start by considering 
the one dimensional case M = 1 [216]. The Hamiltonian in suitable units is 



L-l 



H(ai . . . a L ) = - J\(l)civi+i , 



(261) 



with independent random J\(l). The partition function is 

Z = Tr e -m<n~*L) = Tr exp £ Kx (l) aiai+1 = Tr Y[T t 

i i 



(262) 



with the "bond strength" K\(l) 
transfer matrices 



(3J\{1) ((3 inverse temperature) and the 



y — e KA(0 CT i°"i+i 



(263) 



Since the spin can take only values ±1 on the lattice sites the transfer matrices 
can be rewritten as 2 x 2 matrices with the matrix elements given by (±, l+l \ 
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T\±,l), 



(264) 



With the definition of the Kramers- Wannier dual k* of the bond strength k 
1 



K 



In tanh k 



(265) 



the transfer matrices may be written in the diagonal form 



Ti 



sinh2^(/) 



1/2 



-k*(0ct 3 



(266) 



Here, crj (j — 0, 1, 2, 3) are the complete set of Pauli matrices denoted earlier 
as 1 2 , <J X , & y , and a z (1 2 two dimensional unit matrix). 

In the two dimensional case the Hamiltonian is 

L-l M L-l M 



Z=l m=l 



=1 m=l 



and we assume periodic boundary conditions in the direction of M such that 
M + 1 — > 1. Also in this case one can write the partition function as a product 
of transfer matrices, Eq. (262), but now these consist of two factors 



Ti = vm . 



(268) 



The matrix W\ is the obvious generalization of Eq. (266), 



M 



w t = n 



m=l 



nl/2 



sinh2/^(/,m) 



exp 



M 



m=l 



(269) 



but with the 2M x 2M matrix 



(T <g) <X 



C"3 



CO , 



(270) 



which is the direct product of M — 1 unit matrices and cr 3 at position to. The 
2M x 2M matrix contains the Hamiltonian of the Ith. column of the system 



Vi = exp 



M-l 

£ 

,m=l 



K 5 (/,TO)cr™cr™+ 1 



(271) 
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The formulation in terms of transfer matrices has two implications which are 
of immense importance for practical purposes. First, one can show that it is 
the largest eigenvalue of the total transfer matrix that determines for L — > oo 
the partition function and thus the physical properties of the system [216]. 
Second, one can write the transfer matrix of a system of length L + 1 as a 
function of the transfer matrix of the system of the length L. We will see 
below, that these two properties are essential for numerical evaluation of the 
phase diagram and critical properties. 



8.3 Transforming to Fermions 



The above transfer matrices involve linear and quadratic forms of the spin 
operators. The raising and lowering operators fulfill mixed commutation and 
anti-commutation relations 



[<,4]=0 (m^l), 
K>-} = 1, (272) 



In one dimension, it is well known how to transform these operators to ones 
obeying fermion anti-commutation relations. This is done by the Jordan- 
Wigner transformation [240,241,242,243]. Annihilation and creation operators 
are introduced by 



m—l 

exp in ^2 °+°"- 



m—l 

exp in ^2 o-\g 3 _ 

3=1 



(273) 



a, . 



It is easy to show that these operators obey fermion statistics and 



(274) 



The inverse transformation is 
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m— 1 

exp | iix CjCj 

3=1 



m— 1 

exp in ^2 

3=1 



CjCj 



cl. 



This gives 



_m 



m—l 



exp in ^2 



c j c j 



2c m c m 1 



(275) 



(276) 



The fermion representation of the above transfer matrices is 



M 



m=l 



sinh2/^(/, m) 



1/2 



exp 



-2 ^ <(Z,m) 4c m - - 

m=l V Z 



Vi =exp 



M-l 



m=l 



X] re^mXcJn - c m )(c^ +1 + c m+ i) 

-^(/,M)e^(4-c M )(c t 1 + c 1 ) 



, (277) 



with the number operator N c = Z^m=i c rrt c m- The last term in (277) represents 
the periodic boundary condition in the direction of M. 



The matrix Vi is biquadratic in the fermion operators. It does not conserve the 
number of fermions since it contains terms that create and annihilate fermions 
in pairs. Thus, diagonalization of the transfer matrix in principle should be 
possible via a unitary Bogoliubov-de-Gennes transformation. 

It is illustrative for the disordered case to consider first the ordered limit in 
which K,* x (l,m) = k,* x and /t M (/,m) = First one notes that evenness and 
oddness of N c commute with W and V [216] since the latter contains only 
quadratic forms of the fermion operators, 

[(-1)^,^] = l(-l) N %V] = 0. (278) 



Therefore one can classify the eigenstates of WV according to whether they 
contain even and odd numbers of fermions. One can then write 



V ± = exp 



M 



X/ ( C m ~~ C m)( c m+1 + c m+l) 



(279) 
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with cm+i = =F c i and c M+1 = =Fc{ where ± denote the operators acting on 
the even (with anti-cyclic boundary condition) and odd states (with cyclic 
boundary condition), respectively. The task is then to find the eigenvalues 
and eigenvectors of the transfer operators 



M/2 



±1 



sinh 2k\ 



exp 



Y M 

9 K M X! ( C m ~ C m)( C m+l + C m+l) 



x exp 



M ( t 1 
— 2/t A ^ f c m c m — - 

m=l V Z 



x exp 



M 



9 K M ( C m C m){c m+ i + C m+ \) 
Z m=l 



. (280) 



The diagonalization can be straightforwardly performed by introducing plane 
wave operators via the canonical transformation 



1 



Cm 



(281) 



where q < it is given by the odd and even multiples of n/M for anti-cyclic 
and cyclic boundary conditions, respectively (M is assumed to be even). The 
resulting quadratic forms in the transfer operator which has four eigenstates 
is very reminiscent of the pair Hamiltonian in the BCS theory of supercon- 
ductivity. It can be diagonalized by transforming to new variables 



£± q = cos 4> q r] ±q ± sin (f> q rf Tq , 



(282) 



with tan 2(p q = 2C q / (B q — A q ) containing 



A q = e~ 2K *^ (cosh + sinh cos q) 2 + e 2K * x (sinh sin q) 2 , 

B q = e~ 2K * (sinh sin q) 2 + e 2Kx (cosh — sinh cos q) 2 , 

C q = (2 sinh sin q) (cosh 2«^ cosh k m — sinh 2k* x sinh cos q) . 
With these operators one obtains the diagonal form of the transfer operators 



sinh 2k\ 



exp 



6 9 2 



(283) 
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The energy dispersion is given implicitly by the solution of 

cosh e q = cosh 2k u cosh 2k x — sinh 2k^ sinh 2k* x cos q . (284) 

The four eigenstates in terms of the ^-operators are the vacuum ip , if)± q = 
£,± q ipo and the pair state ip- qq = C-qCl^Po- 

The partition function is given by the L-th power of the largest eigenvalue of 
T ± . The critical temperature T c is defined by k* x = K n . We do not want to 
discuss here how to use the results for obtaining the critical properties. The 
important lesson to keep in mind at this stage is the importance of the parity 
of the fermion number for the structure of the eigenvalue problem. 

8.4 Mapping to a Localization Problem for Non-Interacting Particles 



In the case of a disordered Ising system, the unitary transformation for achiev- 
ing the diagonalization of the transfer matrix depends on the disorder, and the 
different transfer steps are not independent. Therefore, a different approach 
must be used. This consists of mapping the Ising problem to a localization 
problem for non-interacting particles. 

Starting point is the rewriting of the transfer matrix in terms of Dirac fermions 
fulfilling particle conservation [210,212]. In order to achieve this, an identical 
copy of the Ising model is introduced. The corresponding fermions are denoted 
as (ft m , d m . New Dirac fermions are then defined by (suppressing the indices 
for the sake of simplicity) 

c = kf + f + 9 - 9 ] ) d= l -U ] -f-g-g ] ) (285) 



and the inverse of these, 

1 1 

/ = -[ c + c f + i{d + d r )] g = -[c - J + i(d - d r )] . (286) 



The Hamiltonian of the doubled system is the sum of the two Hamiltonians, 
H = (H c + H d )/2. This implies that one deals eventually with products of the 
corresponding transfer matrices. In terms of the new fermion operators they 
are given by 



WfW? = exp 



M 



2 V K* x (l,m)(glf m + P m g m ) 



m=l 
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V^ = exp 



M-l 



2 ^ K^(l,m)(gij m+1 + fl +1 g m ) + & 



(287) 



and they are particle conserving. The term 



b=-K lt (l,M) 



JttN c 



e tnNd )(9lfi + fl9M) 



+ (e 



inN c 



e™ Nd )(gyi + fi9 M ) 



(2* 



is due to the boundary conditions in the direction of M. This contains the two 
boundary operators 



B 



± 



e inN c _j_ e inN d 



(289) 



In order to perform the mapping to the network model it is required that 
the transfer operator conserves the number of fermions. This implies that the 
Hilbert space for first quantization of the transfer matrix must be constructed 
from the Hilbert space of the two copies of the Ising model in such a way that 
the term involving B~ in Eq. (288) vanishes [213]. 

In particular, one may consider the parity operators R defined by 



M 



R = n 



m=l 



(290) 



which change the sign of a complete column of the spins, 

RaiR = -a* , R 2 = l. (291) 



These obviously commute with the transfer operators and have eigenvalues 
±1. The transfer matrix of the two copies of the Ising system written in the 
basis corresponding to the product space of R c + Rd can then be cast into 
a block-diagonal form where the four block matrices describe the transfer of 
amplitudes in the subspaces corresponding to eigenvalues of the parity oper- 
ators operating in the Hilbert spaces of the two copies of the system, (l c , l d ), 
(l c , — Id), ( — l c , Id) , and (— l c , — 1^). Since the eigenvalues of exp (inN) are 
±1 depending on whether M is even or odd, the boundary operators can be 
written as B ± = R c ± Rd- Then, in the subspaces corresponding to the same 
eigenvalues, B~ = 0. The number of fermions is conserved. 
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8.5 The Equivalent Network Model 



In order to construct the first quantized form of the transfer matrix, we start 
by comparing first and second quantized forms of an operator in a Hilbert 
space, say of the dimension N = 2M with a basis Assume that the 

operator can be written in the second quantized form as exp [a^j'IGI/c)^] 
with the creation and annihilation operators aj, and respectively. The 
first quantized form of this operator is the (2M x 2M)-matrix (j\ expG\k). 

By identifying {ai, ol 2 m) = {/i, • • • , /m, gi, ■ ■ ■ , 9m} one finds that the 
first quantization equivalent of the second quantization transfer operator corre- 
sponding to the A-bonds, exp [-2k* x (1, m)](g m f m + f m g m ), is given by (Fig. 34) 

W lm = exp [-2k* x (1, m)<Ti] = cosh [2k* x (1, m)} - a x sinh [2k* x (L, M)\ . (292) 

By analogy, the first quantized form of the /i-bonds in the transfer matrix is 
Vim = exp [2k m (/, m)ai] = cosh [2k^(1, m)\ + o\ sinh [2k m (/, m)\ . (293) 




Fig. 34. Graphical representation of the scattering nodes of the first quantized form 
of the transfer matrix of the Ising model. Left: A-bond; right: //-bond. 



In terms of right- and left-moving amplitudes f) and gj, respectively, the 
transfer operators of the lattice sites can be written as 



If 

J out 
\ fu 



cosh 2k\ — sinh 2k\ 
— sinh 2k\ cosh 2k 




(294) 



and 



t f~ 

J ir 



out 



cosh2K M sinh2/t M 
sinh2/s; M cosh2«; M 



§out 



\ 9ir 



(295) 
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The relations 



0- 3 W f (73 = W 



1 



(TaWs = V 



i 



(296) 



ensure flux conservation (Fig. 35). 



4 




4 



m+1 



4 



m 
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Fig. 35. The network model equivalent to the random bond Ising model. Arrows 
indicate the propagation of flux. The transfer matrix relates flux amplitudes from 
the left to the right. The Ising lattice of spins (arrows) with the exchange matrix 
elements (dashed) are also shown. 



The network model constructed in this way is completely analogous to the 
U(l) network of saddle point scatterings used in the previous chapter for de- 
scribing the critical localization behavior of the quantum Hall system, with 
plaquettes of definite senses of circulation of the fluxes which can be clockwise 
and counterclockwise. Disorder can be introduced in the former network by 
random quenched phases and random saddle point energies. In the present 
case, randomness enters only via the randomness of the exchange integrals 
represented in the Ising model by the bonds. The latter have been mapped to 
the nodes of the network represented by the parameters k and k*. An anti- 
ferromagnetic /z-bond leads to k < 0. An antiferromagnetic A-bond, however, 
leads to a complex n*, where * denotes the Kramers- Wannier dual defined by 
the function, Eq. (265). It implies the relation 



1-1 



1+1 



(-Mr 




(297) 
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This gives a minus sign in the matrices W representing the scattering at the 
nodes. Starting from an Ising system with ferro- and antiferromagnetic bonds 
distributed randomly, one arrives at an equivalent network model with nodes 
W that have random signs. In the node matrices V, the off-diagonal matrix 
elements sinh (2k) acquire random signs. Thus, scattering at the nodes only 
causes real phase factors ±1. Taking into account the topology of the network 
(Fig. 35), this means that a single antiferromagnetic node introduces phases 
7r in both of the anticlockwise plaquettes that are linked at this node. 

This specific property of the disorder determines the symmetry class of the 
model. In the classification scheme introduced recently by Altland and Zirn- 
bauer [244,245] this network model belongs to the class D (see below, Tab. 4 in 
chapter 9). Hamiltonian matrices of this class are purely imaginary, H* = —H. 
The corresponding time evolution operator exp (—iHt) is represented by a ma- 
trix with real elements and can introduce only scattering phase factors ±1, in 
accordance with the above network model of the random bond Ising model. 

The total transfer matrix of the doubled Ising system is then obtained by 
combining the node matrices W and V . It is a 2M x 2M matrix which can be 
arranged in such a way that propagation in one direction is described by the 
first M rows and propagation in the other direction is contained in the M last 
rows. In the former case of the U(l) network the transfer matrix was unitary. 
In the present case it is real and orthogonal. For a random network, it has 
the form exp7L with a real diagonal matrix 7. The fluctuations are oc \fL. It 
can be calculated numerically recursively for very long strips [17,88,89], L — > 
00, and converges according to the theorem by Oseledec [106] to a diagonal 
matrix with diagonal elements 71, ... , 7m, —71, • • • , —7m with the Lyapunov 
exponents 7^. 

The above random-bond Ising network model has been analyzed exhaustively 
numerically in [213] for systems containing as many as M x L = 256 x 5 • 10 5 
lattice sites including scaling analysis along the phase boundary where systems 
of size up to M x L = 32 x 2 • 10 8 lattice sites were used. From the scaling; 
behavior of the data a critical index of v = 1.50±0.03 has been extracted along 
the Nishimori line at the multicritical point. On the other hand, along the 
phase boundary, the scaling analysis gives a critical exponent v T = 4.0 ± 0.5. 

Also other network models, but belonging to the same symmetry class D, 
have been studied [210,211,239]. Most strikingly, these different choices seem 
to yield different localization properties. This will be reviewed in the following 
chapter. 
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9 Symmetry Classes 



Having established above a network model belonging to a certain symmetry 
class corresponding to the quantum Hall transition, one may ask the question 
if one can generalize the model to cover all existing symmetry classes. 

Conventionally, random electron systems are classified into three universality 
classes [144] according to their symmetry with respect to reversal of time and 
rotation of spin. Systems belonging to the orthogonal class have both time 
reversal and spin rotation invariance. Models belonging to the unitary class 
are not invariant under time reversal. The symplectic class contains models 
that have time reversal symmetry but the spin rotation symmetry is broken. 
This is realized in random electron systems with spin-orbit interaction. The 
notation originates from the fact that a distribution of random matrices H 
of the orthogonal, unitary, and symplectic classes is invariant under trans- 
formations, H — > U^HU, O t HO, W R HW with unitary, orthogonal, and 
symplectic (commuting with the antisymmetric Pauli matrix a y ) matrices, 
U, O, W, respectively. Here, O t is the transpose of O, and W R , the time 
reverse of W [144], as defined by = KW T K~ 1 , with the unitary matrix 
K which can be symmetric or antisymmetric, KK* = ±1, corresponding to 
integer spin and half-integer spin respectively. Any random matrix in each 
of the classes can be diagonalized by a unitary, orthogonal, and symplectic 
matrix, respectively. 

The symmetry properties of a random Hamiltonian have strong impact on 
the distribution and statistical correlations between its eigenvalues and the 
distribution of moments of local eigenfunction amplitudes. The level repulsion 
of extended states increases as the symmetry is changed from the orthogonal 
to unitary to the symplectic class. At the same time the tendency to localiza- 
tion decreases, accompanied by an according decrease of local wave function 
fluctuations [144,246]. As reviewed in chapter 4, these properties have been 
studied in detail for the Chalker-Coddington network model. By construc- 
tion, the model belongs to the unitary symmetry class, since the links are 
directed, and correspond to scattering paths which have a well defined chi- 
rality (handedness). This breaks the time reversal symmetry. The question 
arises, whether or not the integer quantum Hall transition as modeled by the 
Chalker-Coddington network model is uniquely characterized by being iden- 
tified as belonging to the unitary symmetry class. This would lead one to 
the conclusion that all quantum Hall type transitions had the same critical 
exponents. 

Since there can be a strong overlap between spin split Landau bands [247], 
the question is of great experimental importance [248] whether or not the de- 
localization transition is sensitive to the spin rotation symmetry. When the 
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time reversal symmetry is broken, the conventional unitary class makes no 
distinction whether or not the spin symmetry is broken. However, for a mod- 
ified Chalker-Coddington model with two spin channels, fixing the scattering 
phases and introducing random SU(2) mixing between the spin channels - 
which corresponds to spin flip scattering by spin-orbit interaction — the crit- 
ical exponent has been found to be v 1.1 [249,250]. This is close to the one 
of classical percolation, v = v p = 4/3, see section 2.4. On the other hand, one 
recovers the critical exponent v = 2.4(2) close to that of the quantum Hall 
phase transition [97,98,249] (Tab. 3) by choosing in addition the scattering 
phases at random. In the latter case, the unitary matrices mixing the two spin 
channels are randomly chosen from the group U(2) = U(l)x SU(2) where 
U(n) denotes the group of complex unitary (n x n)-matrices A, and S stands 
for the special condition that its determinant is one, detA = 1. 

Recently, it has been realized that there can be two additional discrete sym- 
metries in Hamiltonians of disordered electrons. These symmetries have been 
found to give rise to four distinct transitions of the quantum Hall type with 
different critical exponents and different behavior of the quasiparticle density 
of states. In the next section, we introduce these symmetry classes before we 
review the present knowledge about the properties of the corresponding ran- 
dom models. In particular, we will concentrate on the nature of the quantum 
Hall type transitions. 



9. 1 The Additional Discrete Symmetries 

In addition to the invariance under time reversal and spin rotation, there 
can be at least two more discrete symmetry operations in condensed matter 
systems. One is realized in two-sublattice models [251,252,253,254]. It corre- 
sponds to an interchange of the two sublattices together with a sign change 
of the Hamiltonian. A physical realization is a tight binding lattice Hamil- 
tonian with randomness only in the hopping amplitude. Disordered systems 
with such a symmetry belong to the chiral class. 

The other one is the electron-hole symmetry [255] arising from Andreev scat- 
tering in normal metal wires attached to a superconductor [245] or in su- 
perconductors with a gapless quasiparticle density of states [255]. Disordered 
systems with such a particle-hole symmetry are in the Bogoliubov-de-Gennes- 
Oppermann (BdGO) class. 

Both of the discrete symmetries share the property that the eigenenergies of 
the corresponding Hamiltonians come in pairs (—E n , E n ), where E = is 
the center of the band of eigenvalues for the chiral models, and the Fermi 
energy for the BdGO-models. Random Hamiltonians with either of these two 
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discrete symmetries show very different behavior in their density of states and 
transport properties. Their properties are presently the subject of intensive 
ongoing research. In the tables 3, 4, and 5, we have attempted to summarize 
what is presently known about the properties of the non-directed and directed 
network models for all of the existing symmetry classes. 

There are strong arguments, based on a classification of Lie Algebras by Car- 
tan [256,257], that these symmetries make the classification of possible random 
systems complete. In a perturbative renormalization group study in 2 + e di- 
mensions, the corrections to the conductance in the bosonic and the fermionic 
replica formulations have been obtained to 3-loop order for all of the symmet- 
ric spaces [258]. 

There is an ongoing debate whether or not this symmetry classification is 
sufficient to characterize the quantum critical properties of random systems 
[211,259,204,260]. As we will review below, there are some results which do not 
seem to fit into this scheme, since the critical properties sensitively do depend 
on the spatial correlation of the disorder. As discussed above, the conventional 
quantum Hall transition is insensitive to the correlation length of the disorder 
potential, unless the disorder has an extremely long-range correlation length 



In the following, we first discuss the properties of quasi-one dimensional dis- 
ordered wires, and how their transport properties depend on the symmetry 
properties. Then, we consider the critical properties in two dimensions. We 
present an overview of the network models studied for the various symmetry 
classes. 

For a quasi-one dimensional wire, where the localization length £ is larger than 
the wire width L y , a beautiful formula, valid for all of the ten symmetry classes 
has been derived recently by solving the Fokker-Planck equation [262,268] 



Here, mi is the multiplicity of the root on the symmetric space (see below), (3 
the level repulsion coefficient, N the number of channels proportional to the 
width, and s the spin factor. This formula displays explicitly the importance 
of the symmetry properties of the disordered Hamiltonian for the localization. 
We have seen already that the transfer matrix T of the Chalker-Coddington 
model fulfills the commutation relation Eq. (91), which is a result of current 
conservation. Thereby, transfer matrices form a group. Since they are param- 
eterized by continuous variables, such groups are called Lie groups. For the 
Chalker-Coddington transfer matrix this group is G=U(iV, N), the group of 
noncompact 2N by 2N matrices, which fulfill Eq. (91). For the other ensem- 



[84]. 




(298) 



107 



Wigner-Dyson Class 



property 


orthogonal 


unitary 


unitary 


symplectic 


TRS 


yes 


no 


no 


yes 


SRS 


yes 


yes 


no 


no 


Cartan class 


AI 


A 


A 


All 


P 


1 


2 


2 


4 


mi 


1 


1 


1 


1 


s 


1 


1 


2 


1 


NDNM 


insulator 


insulator (a) 


insulator 


MIT [76] 


DOS (E -► 0) 


constant 


constant 


constant 


constant 


phase factor 




U(l) 


U(2) (b) 




in DNM 










PT in DNM 




IQHT 


IQHT (c) 




V 




2.5(5) [31] 


2.4(2) [97,98] 
[20] (c) 


2.8(1) [Sec. 9.2] 


a 




2.261(3) [137] 




2.174(3) [261] 


A c 




l/vr(a - 2) 


2/vr(a - 2) 


1.83(l)[Sec. 9.2] 


equivalent models 




RSVDM 


RSVDM 





Table 3 

Classification of network models (NM) according to their symmetry: the Wigner- 
Dyson class with no discrete symmetry. Abbreviations: TRS time reversal symme- 
try, SRS, spin rotation symmetry, (3 level repulsion exponent, mi multiplicity of 
long root on symmetric space [257], s spin factor, IQHT Integer Quantum Hall 
Transition, MIT metal-insulator transition, DOS density of states, (N)DNM (non- 
directed network model, PT phase transition, v critical exponent, «o position of 
maximum of f(a), A c MacKinnon-Kramer parameter, RSVDM Dirac model with 
both, random scalar and vector potentials; (a) electrons in a random magnetic field 
belong to this unitary class. The U(2) network model proposed by Chalker and Lee 
[97] describes electrons in random magnetic field with large correlation length, and 
does not show a derealization transition; (b) spin degenerate levels with spin orbit 
interaction [97]; (c) two IQHT's at two distinct energies [97]. 
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Bogoliubov-de-Gennes-Oppermann Class 



property 


orthogonal 


unitary 


unitary 


symplectic 


TRS 


yes 


no 


no 


yes 


SRS 


yes 


yes 


no 


no 


Cartan Class 


CI 


C 


D 


Dili 





2 


4 


1 


2 [262] 


mi 


2 


3 





[262] 


s 


1 


1 


4 


2 


transport 


SI 


SI 


TM 


TM (a) 


DOS (E -> 0) 


\E\ [245,263] 


\E 2 [245] 


\E\ (b) 

\/HV\e\) 

(c) [264] 


|^| [245] 


phase factor 




SU(2) 






PT in DNM 




SQHT 


TQHT 




JJUo at Ur 




\hi\ 1 |ZOO,ZOt)J 






V 




4/3 [265] 


(d) 


? 


a 




2.137(3) [184] 


(d) 


? 


Ac 




l/vr(a - 2) 


(d) 


? 


equivalent models 




CP 


RBIM 





Table 4 

Classification of network models according to their symmetry: the Bogoliubov-de- 
Gennes-Oppermann class with particle-hole discrete symmetry. SI spin insulator, 
TM thermal metal, (a) critical metal for uncorrelated disorder of arbitrary strength, 
MIT for correlated disorder, (b) DOS in the localized phase, (c) DOS in the metallic 
phase, (d) line of critical points, on the Nishimori line v = v v = 4/3 [210], SQHT 
spin quantum Hall transition, TQHT thermal quantum Hall transition, CP classical 
percolation model, RBIM random bond Ising model [210]. 



bles additional conditions are imposed, restricting the corresponding transfer 
matrices to be in different Lie groups. Now, all the powerful tools of represen- 
tation theory of Lie groups [269] can be used to analyze these transfer matrices 
systematically. There is a further simplification, that a transport property like 
the conductance Eq. (147) does depend only on products TT + and T + T 
[160,162]. As a result, there is a subgroup K of matrices k by which T can be 
multiplied from right T — > kT or left T — > Tk, without changing the conduc- 
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Chiral Class 



property orthogonal unitary symplectic 



TRS 


yes 


no 


yes 


SRS 


yes 


yes (no) 


no 


Cartan class 


BDI 


AIII 


CII 




1 


2 


4 


mi 


1 


1 


1 


s 


1 


1(2) 


1 


transport (E = 0) 


metal 


metal 


? 


DOS (E -► 0) 


p(E) 


P{E) 


? 


phase factor 


2 x Z 2 


2 x U(l) 


2 x SU 


PT in DNM 




Yes (a) 





V 



Q 

Ac 

equivalent models RXY [267] RVDM 



Table 5 

Classification of network models according to their symmetry: The chiral class 
with bipartite discrete symmetry (two sublattices). Abbreviations: p(E) = 
exp(-c| ln^l 1 ^)/!^! [253] (x = 3/2 for BDI [204,205] and possibly also for AIII, CII 
[205]), RXY random XY model, RVDM Dirac model with random vector potentials, 
(a) The phase diagram is shown in Fig. 46. 



tance through the system. For the Chalker- Coddington model this subgroup 
is found to be K = U(iV) x U(iV), the group of 2N times 2N matrices which 
consist of two blocks of unitary N by N matrices. Thus, this subgroup is not 
relevant for a physical property like the conductance, and can be divided from 
the group G. The resulting group G/K is for the Chalker-Coddington model 
G/K = U(N,N)/(U(N) x U(N)). This is a so called simple Lie group and 
is also called a symmetric space. There is a full classification of all symmetric 
spaces by Cartan. The way to characterize a symmetric space is to consider 
the corresponding Lie algebra, that is the properties of matrices A, defined by 
T = exp(A). Now, one defines eigenvectors X of A, by the commutation rela- 
tion [A, X] = aX. The degeneracy of the eigenvalue a = defines the rank r 
of the Lie algebra. The set of eigenvectors Hj, i = 1, r with eigenvalue a = 
forms an Abelian subalgebra, the so called Cartan Algebra, since [Hj, Hj] = 
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for all % = l,...,r. Since every matrix A commutes with itself, [A, A] = 0, 
it can be written as a linear superposition A = J2i=i A,Hj. There are N — r 
eigenvectors E$, % = r + 1, N of Hj, % = 1, r with nonzero eigenvalues 
a r . The vector of nonzero eigenvalues a = (a±, a r ) is called the root 
of the algebra. These roots can be characterized by their length, and their de- 
generacy, the so called multiplicity. If there is only one eigenvector to a root a 
it is called a simple root. Now, it can be shown that the parameters entering 
in the above formula for the quasi-one dimensional localization length, Eq. 
(298), (3 and mi are the multiplicities of roots on the respective symmetric 
space G/K [256,257]. 

We have listed their values in tables 3,4, and 5. The factor s accounts for 
the fact that spin flip scattering from magnetic impurities mixes the two spin 
channels, and thereby can double the effective number of channels. Thus, 
s = 1, without spin scattering, and s = 2 with spin scattering. The integer N 
is the number of transverse channels in the quantum wire, and I the elastic 
mean free path. This formula, Eq. (298) coincides for the ordinary symmetry 
classes with the well-known result derived independently by Efetov and Larkin, 
and by Dorokhov, namely that the localization length is proportional to the 
symmetry parameter f3 [87,270,271,272,273,274,275]. For the ordinary and the 
chiral classes (3 takes the well known conventional values f3 — 1, 2, 4, mi — 1 
for the orthogonal, unitary and symplectic class, respectively. 

For the Bogoliubov-de-Gennes-Oppermann class, however, (3 and s have quite 
unexpected values. Moreover, in the classes D and Dili, the multiplicity mi 
vanishes. This implies a diverging localization length. For weak Gaussian dis- 
order the conductance in classes D and Dili has the broad distribution typical 
for a critical state [262]. Accordingly, in a gapless superconductor with bro- 
ken spin symmetry the quasiparticles are delocalized for weak disorder. Since 
both, the charge and the spin, are not good quantum numbers in such a sys- 
tem, the metallic character of the quasiparticles can only be measured due to 
their finite contribution to transport of heat which can be measured via the 
thermal conductance [264]. 

In order to understand this better, it is enlightening to generalize briefly the 
concept of the quantization of the electrical — charge — conductance to the 
thermal and the spin conductance. 

The defining relations of charge, spin and heat conductivities are 




(299) 




(300) 
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(£) V(* B T) (301) 



with electron charge — e, voltage V, Bohr magneton /ib, Boltzmann constant 
and temperature T. The charge, spin and heat conductivities are denoted 
as a c , <j s and k, respectively. One can easily verify that these are consistent 
with the usual linear response definition apart from rewriting the current 
densities as the responses to gradients of energies instead of fields. We also 
recall that thermal and electrical currents are related via the Wiedemann- 
Franz law 

Hi 7T /ltd ,„ , 

(302) 



2 



a c T 3 e 



The quantization of the charge conductance in terms of cr = e 2 /h implies 
that the conductance per charge is given by the ratio of the elementary charge 
and the Planck constant, = e /h- Analogously, the quantum of the spin 
conductance will then be given by the ratio between elementary spin and the 
Planck constant, <7q/ (/is^/ '2) = (h/2)/h. This implies 

^0=^6 — ^ — • (303) 



The quantum of the heat conductance should be given by the ratio between 
the heat per free particle Qo and the Planck constant. The heat per charge 
carrier can be deduced from the thermopower S 



eST = —k B T (304) 




where the last factor accounts for the fact that only electrons in a temper- 
ature window k B T near the surface of the Fermi sea contribute. Note that 
the factor 7r 2 /3 is the first Sommerfeld expansion parameter in the expansion 
of the energy integral around the Fermi surface, which does not depend on 
dimensionality. From this one concludes 

Qo = y*;BT, (305) 



and from 

£ = % (306) 
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follows that 



«o = T i- (307) 



In the following, we apply this to the case of a quantum Hall system. We use 
units such that /ib = 1. 

When a temperature difference smaller than the Landau gap is applied to two 
opposite edges of a quantum Hall bar, there is no thermal transport between 
the edges, if the Fermi energy is pinned to localized bulk states. Still, the 
total heat current perpendicular to the temperature gradient is non- vanishing, 
since the quasiparticles on the hot edge of the sample carry more heat in one 
direction than the quasiparticles on the cooler edge in the opposite direction. 
This results in a finite thermal Hall conductance n xy defined by 

4° = Kxy A H T, (308) 



where AhT is the temperature difference between the edges and I® the total 
heat current. One finds that one edge channel contributes towards the thermal 
Hall conductance with the quantized amount given in Eq. (307) [276], 

7T 2 k 2 T 

**y = Yh- (309) 



As the charge Hall conductance, the thermal Hall conductance is non- zero only 
for broken time reversal symmetry and broken parity. Non-vanishing values 
can thus only be obtained in directed network models. 

There has been some debate [239,264,277,278] if this critical metal phase of 
classes D and Dili persists due to the particular symmetry even at strong 
disorder as argued in [268], or if there could be a regime of localized phases 
characterized by a quantized thermal Hall conductance parameter, in units of 
the thermal conductance quantum k , 

T 3hK xy {airi\ 

**» = ' (310) 



with quantum Hall type transitions between phases with aj y = and cr^ y = 1 
[276]. One could refer to this transition, accordingly, as the thermal quantum 
Hall transition. 

Recently, several investigations have concluded that the phase diagram of 
models belonging to class D critically depends on the spatial correlations of 
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Fig. 36. Scheme of the quantum Hall flow diagram of class D systems in the plane of 
the local thermal conductance and thermal Hall conductance parameters, o xx and 
aj y , respectively, for correlated disorder where an insulating and a quantum Hall 
phase can exist [239,264,278]. For weak disorder (large a xx ), the critical metal phase 
is stable for all values of the Hall conductance. For uncorrelated disorder, this is the 
only possible phase for any value of aj x . Correspondingly, a T * = for uncorrelated 
disorder. For strong correlated disorder there is a flow towards an insulating phase 
for a xy < 1/2, and a flow to a phase with quantized Hall thermal conductance for 
a xy > 1/2) respectively. 



the disorder potential and that the metal phase persists to strong disorder only 
for uncorrelated, white noise disorder [211,259,262,279]. If there is at strong 
disorder a transition to a localized phase, and the time reversal symmetry 
is broken, as it is the case in class D, there can be a finite thermal Hall 
conductance n xy . Figure 36 shows schematically the suggested two-parameter 
flow diagram for correlated disorder [239,264,278] as a function of the thermal 
conductance parameter 

^ = ^|T' (311) 



and the thermal Hall conductance parameter a xy . Such a thermal metal in- 
sulator transition and the corresponding thermal quantum Hall transition are 
expected to have very different critical properties from the integer quantum 
Hall transition (Tab. 4). 

When the time reversal symmetry is broken but global SU(2) spin rotation 
symmetry is not, then the quasiparticles in a superconductor can be studied 
both by their contribution to thermal and to spin transport. In two dimensions 
this class, denoted as C in Tab. 4, can be studied in d-wave superconductors. 
Here, the order parameter vanishes at four nodal points on the Fermi surface. 
This allows to study the peculiar properties of the quasiparticles which occur 
close to the Fermi energy due to the particle- hole symmetry in class C. As seen 
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from Eq. (298), the quasiparticles in a quasi-one dimensional wire belonging 
to class C are localized. From results of perturbative renormalization group 
studies, one can conclude that they stay localized in two dimensions [280]. 

However, in close analogy to the quantum Hall transition, a localization- 
delocalization transition is possible in which the quantized Hall conductance 
for the spin, <r xy , is quantized. 

The spin transport is caused by a gradient in the perpendicular magnetic field 
which results in a gradient in the Zeeman energy. In linear response the spin 
current is [281] (Eq. (300)), 

f = a s VB, (312) 



where a s is now a 2 x 2- matrix with components a s xx and a s xy . When the Fermi 
energy is pinned to localized bulk states, the edge states carry a finite spin 
current perpendicular to the gradient of the magnetic field, 

./•; = o^VyB. (313) 



When the quasiparticles carry the spin 1/2, the spin Hall conductance, o s 
turns out to be quantized in two dimensions in multiples of the quantum of 
the spin conductance a (Eq. (303)), 



(H/2)' 

(T s xy = mflB }U-L i (314) 



where m is found to take only values 0, ±2 [278], very similar to the quan- 
tization of the charge Hall conductance in units of e 2 /h. The spin angular 
momentum h/2 is substituted for the electron charge e. Changing the Fermi 
energy one finds spin quantum Hall plateau transitions where the spin Hall 
conductance changes by two units. This transition is in a new universality 
class, and called the spin quantum Hall transition (SQHT). When a strong 
Zeeman term is introduced which reduces the SU(2) symmetry of the spin to 
U(l), this transition splits into two, each of them being in the usual univer- 
sality class of the integer quantum Hall transition [265,266,280,281,282]. The 
spin quantum Hall transition can also be probed by the thermal conductance, 
which in a metal phase is related to the spin conductance by the analogue of 
the above Wiedemann-Franz law 

K -** kl (315) 



a s T 3 (h/2) 2 ' 

Recently, a model belonging to class C has been mapped on the classical 
percolation model [265,283]. This made it possible to determine the exact 
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critical exponents of the spin quantum Hall transition, like the exponent of 
the divergent localization length, v = v v = 4/3. The mapping is reviewed 
below, when the Class C network model is introduced in more detail. We note 
that disordered d-wave superconductors also have been argued to be in the 
chiral class AIII, yielding delocalized quasiparticles close to the Fermi energy 
[263]. This situation arises for a slowly varying disorder which has negligibly 
small Fourier components for scattering between the four nodal points where 
the gap in the quasiparticle density of states is vanishing. For each node, the 
quasiparticles independently can be modeled by Dirac fermions. In a magnetic 
field, which breaks the time reversal invariance, this leads to the sub-lattice 
class AIII [263]. 

Sub-lattice models exhibit some similarity to a quantum Hall transition, since 
they have one extended state in the middle of a band of localized states. Their 
Hamiltonian is given by a tight binding model on a two dimensional bipartite 
lattice. There are no on-site potentials, and the hopping matrix elements are 
random. Thus, they have perfect particle-hole symmetry. In contrast to the 
models for the integer Quantum Hall Effect, two-sub-lattice models have a 
divergent density of states in the center of the band where the extended states 
appear. For the unitary (AIII) and the orthogonal (BDI) chiral classes, the 
density of states diverges, 

p(E)K^-e-^ 1,x , (316) 



with constant c. By using field theory, Gade obtained x = 2 for both classes, 
AIII and BDI [252,253]. By considering random Gaussian surfaces, and nu- 
merically, the value x = 3/2 was obtained for class BDI [204]. This result for 
BDI was also derived with the supersymmetry method [205] and was argued 
to be a direct consequence of the freezing transition of the dynamical expo- 
nent z [204,205,284,285,286], which is related to the fact that the multi-fractal 
spectrum is bounded [286]. While it has been suggested that this freezing tran- 
sition is not sensitive to time reversal and spin symmetry and therefore the 
value x = 3/2 could be valid for the chiral classes, AIII and CII, as well [205], 
a derivation is pending. 

The model of Dirac fermions in a random vector potential (Sect. 7.2) is another 
member of the symmetry class AIII. Although a divergence of the density of 
states has been found here as well, the corresponding exponent is varying 
continuously with the disorder strength [120]. When moving away from the 
band center, the sub-lattice models acquire the properties of the ordinary 
models, and the states are localized for the orthogonal and unitary models. 
Sub-lattice models with directed bonds have recently been established [260]. 
Their critical properties still have to be studied in detail. So far, the symplectic 
sub-lattice class has been studied only briefly [287]. 
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According to the above description, the derealization transition of the Chalker- 
Coddington network model without particle-hole symmetry should be insensi- 
tive to the spin rotation symmetry, since the conventional unitary class makes 
no distinction whether or not that symmetry is broken. This argument does 
not take into account that the lifting of the spin degeneracy can result in 
a change of the statistically relevant density of states and thereby also a 
change in the localization behavior. This has been first noted for localization 
in quasi-one dimensional disordered wires without time reversal symmetry, 
where the localization length is doubled when the spin degeneracy is broken, 
Eq. (298) [270,275]. Recently, it has been argued within a semi-classical per- 
colation model [288], that spin-orbit scattering can change the universality 
class from the one of the conventional quantum Hall transition to the one 
of classical percolation. In short, the argument is that the Hamiltonian with 
spin-orbit scattering and random scalar potential reduces — in an adiabatic 
approximation — to the Hamiltonian of electrons with opposite spins mov- 
ing in the effective potential which is the sum of the scalar potential and the 
positive/negative locally varying random Zeeman fields for the electrons with 
spin up and spin down, respectively. Thus, the electrons with opposite spins 
do see different randomness. Close to saddle points in the effective potential, 
the electrons have the choice of either to tunnel while keeping their spin or 
to flip their spin using the non-adiabatic part of the spin-orbit interaction, 
and then continue to propagate in the potential landscape corresponding to 
the opposite spin. For a potential with a sufficiently short correlation length, 
the potential landscapes for the electrons with spin up and spin down should 
be sufficiently different that this argument applies. This argument suggests 
that the critical exponent of the localization length in that situation is close 
to the one of classical percolation, v = v p = 4/3 [288]. The obvious similar- 
ity of this model to models of the Bogoliubov-de-Gennes-Oppermann class, 
when identifying the spin flip scattering with the anomalous pairing ampli- 
tude A, coupling electrons with holes may be used as a heuristic explanation 
for the change of the universality class to the one of the classical percola- 
tion transition, in class C models. Indeed, the network model of class C with 
SU(2)-scattering phases at the links has been first studied as a model of spin 
degenerate Landau levels with negligible random potential and strong spin- 
orbit scattering [249]. Another realization would be Landau levels with strong 
magnetic impurities [254]. This was found to have a single quantum critical 
point with v ss 1.1 [249]. 

In order to study quantitatively the localization properties of the metal- 
insulator transitions and the quantum Hall transitions it is necessary to for- 
mulate specific network models for the various symmetry classes. In the next 
section, we introduce non-directed network models for the three ordinary sym- 
metry classes, of which only the symplectic one is expected to have a metal 
insulator transition. Next, we review unitary network models with two spin 
channels that describe the transition of quantum Hall systems with mixing 
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between spin split Landau levels. Then, we review the network models for 
transitions of the quantum Hall type belonging to the Bogoliubov-de-Gennes- 
Oppermann class. Finally, we provide an overview on the network models for 
sub-lattice systems. 



9.2 Non-directed Network Models 



In order to construct network models belonging to the orthogonal and sym- 
plectic classes, one should assume links that are not directed. Following early 
work on the scattering matrix formulation of the scaling theory of Anderson 
localization [289], an example of such a model has been established and stud- 
ied by using the real space renormalization group method by Shapiro in 1982 
[290] . Recently, such network models have been explored numerically for all of 
the three conventional symmetry classes [261,291]. The results of the scaling 
theory of localization have been reproduced: a localized phase exists for the 
unitary and orthogonal classes, and a metal-insulator transition occurs in the 
symplectic class. The corresponding S matrices are unitary, symmetric, and 
symplectic [87]. 

The fact that the links are not directed, but nevertheless the time reversal 
symmetry is broken for non-directed unitary network models, makes these 
models equivalent to two dimensional random fermions in a weak magnetic 
field that does not affect the classical motion of the electrons. Accordingly, 
extended states are not expected to exist. This has been confirmed by extensive 
numerical studies [291]. 

In non-directed networks, the scattering matrices representing the nodes relate 
four incoming to four outgoing currents (Fig. 37), 

















— S ort h 

















(317) 



With time reversal invariance, the S matrix is symmetric [87]. By assuming 
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isotropy of the scattering, it can be cast into the form [261] 



'orth 



I r d t d^ 
d r d t 
t d r d 
d t d r 



(318) 



where t, r and d denote transmission, reflection and deflection amplitudes, 
respectively. It is assumed here that the deflection to the left and to the 
right are equal, as it should be for a network without chirality. The scattering 
parameters are complex and satisfy the conditions 



|r| 2 + |t| 2 + 2|d| 2 = l, 

d*r + dr* + d*t + dt* = 0, (319) 

2|rf| 2 + r*t + rt* = 0. 
^From the 1st and the 3rd condition one finds the relations 



|r| 2 + |t| 2 <l, l<|r| + |t|. (320) 

The second inequality is violated when the time reversal symmetry is broken. 

Due to random phases attached to links, we can always set d to be a real 
positive number. Then, once \r\ and \t\ are given, r,t and d are uniquely 
determined by Eq. (320). 
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Though the transfer matrix, at first sight, seems to be complex, the calculation 
can be performed with only real (or quaternion-real) numbers after a proper 
unitary transformation [292]. 

In order to construct a model belonging to the symplectic class, we first es- 
tablish a scatterer that rotates the direction of the spin [261]. The S matrix 
relates the currents with spin up (|) and spin down (j) according to 



( N 








S S y m p 











(321) 



with 



•>symp 




(322) 



This choice of symplectic scattering matrices mixes only different spin chan- 
nels moving in the same direction, but does not have backscattering matrix 
elements. Backscattering occurs thus in this model at the orthogonal nodes 
only. Here, q is a quaternion-real number, 

3 

Q = J2lkTk, T k = ia k (k = 1, 2, 3) r = 1 . (323) 

k=0 



The coefficients q k are real numbers that satisfy 

E dl = 1 (324) 

k=0 



and the 2x2- matrices a k are the Pauli spin matrices. The quaternion conjugate 
is denoted by q. 

3 

q = qol-^qiTi. (325) 
i=i 



The strength of the spin orbit scattering is characterized by the distribution 
of qi. When q = 1, q 1 = q 2 = q% = there is no spin orbit scattering. The 
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spin rotation is most random when 



go + m = e ia cos(3 

qi + Hi = e 47 sin P 



with a and 7 distributed uniformly between [0, 2ir), and (5 distributed accord- 
ing to the probability density, 

, sin(2/3) < < ir/2 
P{f3) = { V ' ~ 1 . (327) 
otherwise 
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Fig. 38. Schematic view of a scatterer that provides spin rotation. Solid arrows 
correspond to spin-up currents. Broken arrows indicate spin-down currents. 

As before, one attaches random phases between the nodes described by the S 
matrix 



/ e ^ 



'link 







(328) 



The resulting network is shown in Fig. 39. This model is similar in spirit but 
different in the details as the tight binding Hamiltonian on a square lattice 
including spin-orbit scattering proposed in [293,294,114]. A similar model also 
has been proposed in [295]. 

Figure 40 shows some results for the localization length X(L) in a two dimen- 
sional strip of the width L. The renormalized localization length A = X(L)/L 
as a function of |r| is plotted for different L. When |r| is smaller than 0.61, A 
increases with increasing L. It is decreasing function of L for |r| > 0.63. This 
behavior indicates a delocalization-localization transition [17]. For the sym- 
plectic symmetry there occurs a metal-insulator transition in two dimensions 
as has been conjectured before [296]. By analyzing the data using the finite 
size scaling method the critical exponent can be extracted (cf. Tab. 3). The 
estimates with 95% confidence interval, v = 2.81 ± 0.05, A c = 1.836 ± 0.027, 
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Fig. 39. Network of nodes realizing the symplectic class. White boxes describe 
scattering with orthogonal symmetry. Shaded boxes rotate spin, thus breaking spin 
rotation invariance. 

agree within the uncertainty with the result of numerical studies done for the 
tight binding Hamiltonian, v = 2.746 ± 0.009, A c = 1.843 ± 0.001 [114,115]. 

To construct orthogonal and a unitary nondirected network models, we assume 
in S symp q = 1. In the unitary case, S un k is no longer symmetric, 



o unitary 
3 link 



e** 
0' 



(329) 



In both cases, we have spin degeneracy, and the dimension of the transfer 
matrix is halved. 



9.3 Network Models of the Bogoliubov-de-Gennes-Oppermann Class 



Class C 

As mentioned above, systems belonging to the Bogoliubov-de-Gennes-Opper- 
mann class have a particle-hole symmetry, and the spectrum of eigenvalues 
consist of pairs, {—E n) E n ). The zero of energy, E = 0, is the chemical poten- 
tial in the superconductor. Class C describes quasiparticles in a spin-singlet 
superconductor in which time-reversal symmetry is broken, but spin rotation 
symmetry is conserved. This corresponds to systems with negligibly small Zee- 
man splitting [245] . Far away from the Fermi energy, the transport properties 
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Fig. 40. MacKinnon-Kramer variable A as a function of \r\ for \t\ = 0.6; \d\ is 
determined according to Eq. (320); L=6,8, 12, 16,24,32 and 48; precision of the data 
is 0.2% except for L = 48 where it is 0.4%. Curves intersect approximately at 
\r\ = 0.62, which indicates the presence of an Anderson metal-insulator transition 
in two dimensions. 

show a crossover to the conventional unitary class A. Therefore, class C be- 
havior can only be studied in gapless superconductors or in metals closely 
attached to superconductors, where quasiparticles close to the Fermi energy 
exist. As summarized in table 4, the quantum states of quasiparticles in nondi- 
rected models of the class C are localized. This phase is named spin insula- 
tor, because the insulating behavior reveals itself only as a vanishing spin or 
thermal conductance, since the charge is not conserved in a superconductor. 
When parity symmetry is broken, as in directed network models, there can 
be a quantum Hall type transition to a phase with a critical state and an 
integer-quantized spin Hall conductance a s xy . This is called the spin quantum 
Hall transition (SQHT). 

Thus, we are now looking for a generalization of the Chalker-Coddington 
model which satisfies the particular symmetries of class C [265,266,282]. In 
the simplest case, one has two-component wave functions, which propagate on 
directed links through the lattice, and which scatter at nodes between adjacent 
links. All the nodes in the lattice have two incoming and two outgoing links. 
It turns out [265], that the SQHT can be determined from the properties of 
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the perimeters, or hulls, of classical percolation clusters on a two dimensional 
lattice constructed of such nodes and links. This classical problem is much 
simpler than the original quantum problem, and can be solved exactly, for 
example by mapping it on an integrable super-spin chain [265] which will be 
introduced below in the chapter 10. 

In the following, we review the models for class C, as introduced by Beamond 
et al. [266]. A similar network model with SU(2) scattering phases at the links 
has been studied before in the context of the Quantum Hall Effect of spin 
degenerate Landau levels [249]. We consider a set of nodes n connected by 
links I. Each node is assumed to be of the degree four, such that at each 
node two directed links enter and two leave. A two-component wave function 
is assumed to propagate along each link. This propagation may be described 
by a unitary evolution operator U, which describes the evolution of the wave 
function one unit forward in "time" , as the particle moves from a given link to 
a neighboring one. The evolution operator plays a similar role in defining the 
network model as the time evolution operator corresponding to a Hamiltonian, 
for example, of a tight binding model. This procedure has been discussed in 
detail before for the U(l) network model in [121,189] (Sect. 7). 

The model can be constructed by associating each link / with a unitary 2 x 
2 matrix U/. This matrix specifies the phase accumulated when traversing 
the link. Each node n is assumed to be represented by a scattering matrix 



where I2 is the 2x2 unit matrix. The S matrix describes scattering at the 
node from the incoming links to the outgoing ones. If the network has M 
links (and therefore M/2 nodes), then U is an N x N matrix, with N = 2M. 
It consists of M/2 blocks, each associated with a particular node and of size 
4x4. The block at the node n has the form 



where (1, 2) and (3, 4) label the links which are incoming and outgoing at the 
given node. 

So far, we only have assumed the links to be directed, so that time-reversal 
symmetry is broken. To identify a network model to belong to class C, one can 



(Sect. 3.2) 




(330) 




(331) 
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start [282] from the defining property of a Hamiltonian H with this symmetry 
[245] 



H* = -a y Ha y , (332) 



with the Pauli matrix a y acting on the spin variables. The operator H* is the 
complex conjugate of H. With this Hamiltonian, and using that U can be 
written as 

U = e~ iH , (333) 



Eq. (332) implies 

U = <t v U*a y . (334) 

/From this, an equivalent restriction follows for the phases of the links, 

= a y \J*a y (335) 

which are therefore unitary Sp(2)-matrices equivalent to SU(2)-matrices, where 
Sp(2N) is defined to be the symplectic group of 2N by 2N matrices A which 
fulfill the condition A T a y A = a y [266]. The wave function consists of two 
components in each link. The space of states on each link may be viewed as 
consisting of a two-dimensional subspace, within which a y operates. 

One possibility to introduce randomness into the model consists of assuming 
the link phases to be random variables drawn independently from a distribu- 
tion which is uniform on the invariant (Haar-) measure of Sp(2). The ensemble 
average of a physical quantity in the network model, denoted by (...), is the 
mean with respect to this measure. 

As reviewed in Sect. 5.1, the conductance of an open network of a finite length 
can be expressed by the multi-channel Landauer formula in terms of the trans- 
mission probability. Specifically, the spin conductance measured between two 
contacts in units of (h/2) 2 /h is given by Eq. (147), 

# = Trtt f . (336) 



Here, t is the rectangular N x N transmission matrix containing the matrix 
elements for z — 1 of the propagation operator 

GW = ^ (337) 
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between the incoming and outgoing link states k and l Q . 

One can now show that the disorder average Tr (G(z)) taken over the Sp(2)- 
phase factors of the Green function 

G(z; 1,1') = -zU)- 1 ^') (338) 



can be expressed in terms of classical paths [265,266]. 

To illustrate this, we define for the same network a classical scattering problem 
as follows. The scattering at each node may be decomposed into disconnected 
processes in two different ways, (1 — > 3, 2 — > 4) and (1 — > 4, 2 — > 3) (Fig. 41). 




Fig. 41. Decomposition of the scattering at a given node: S matrix elements cos9 n 
and ±sin# ra are associated with the transitions (1,2) — > (3,4) as indicated on 
the left. The two decompositions are weighted with factors p n — cos 6 n and 
1 — Pn = sin 2 9 n [Figure taken from [266]]. 




Fig. 42. The network model for the spin Quantum Hall Effect. Scattering amplitudes 
±cos6> ra and ±sin# n are indicated by ±c and ±s [Figure taken from [266]]. 



A two dimensional model exhibiting the spin Quantum Hall Effect is obtained 
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by taking the network to be the so-called L- lattice shown in Fig. 42. The two 
possible classical decompositions of a node may be associated with the pres- 
ence or absence of a bond, with probabilities p and 1—p between neighboring 
sites on an associated square lattice [265]. The latter is rotated by 45° relative 
to the L-lattice, and has a lattice spacing increased by a factor y/2. In this 
model, closed loops of the classical problem form interior or exterior hulls of 
bond percolation clusters on the larger lattice. It is known that such loops are 
finite with characteristic size < oo except at the critical point, p = p c , 
which for bond percolation on the square lattice occurs at p c = 1/2. When 
approaching the critical point, £ diverges as £ ~ \p — p c \~ v with v = u v = 4/3 
(Chap. 2), while at the critical point the distribution of the hull lengths is 
P(L) ~ L~ y at large L, with y = 8/7. 

The quantum localization length in this class C model diverges with the same 
value of the exponent, v — v v — 4/3, as the plateau transition is approached. 
The density of states 

p(E) — — —TmTr{G(E)) (339) 

7T 



varies for small E as 

p(E) ~ E 2 (340) 

in the localized phase, and as 

p(E) ~ \E\ 1/r (341) 



near the critical point (Tab. 4). 
Class D 

As described in the introduction to this chapter (Tab. 4 and Fig. 36), mod- 
els in class D have particularly rich phase diagrams in two dimensions. The 
symmetry is realized in superconductors with broken time reversal invariance, 
and either broken spin rotation invariance, as in d-wave superconductors with 
spin-orbit scattering, or spin-less or spin-polarized fermions, as in certain p- 
wave states. Since the spin rotation symmetry is broken, and the charge is 
not conserved, the quasiparticles can only be measured by their contribution 
to thermal transport. Equation (298) shows that quasiparticles are delocal- 
ized in quasi-one dimensional quantum wires belonging to class D. In the two 
dimensional case, the flow under the renormalization group is towards larger 
thermal conductance [264,277,278,297]. Thus, there is a phase in which there 
is a diverging density of extended eigenstates at zero excitation energy (Tab. 4 
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[264]). A superconductor described by this model is in a thermal- metal phase. 
A phase with localized quasiparticles is a natural possibility, although it has 
only been found so far for models with correlated disorder [239] which gives 
rise to a phase with a quantized Hall conductance (Fig. 36). 



In analogy with the network models belonging to class C introduced in the 
previous section, one can formulate a directed network model with the sym- 
metries of class D [239]. Disorder appears in the network model in the form 
of random scattering phases, and the symmetries of the class D restrict scat- 
tering phases to the values and n, multiplying the wave function at each 
node at random with phase factors ±1. Remarkably, within this framework, 
different particular forms of disorder result in quite different physical behav- 
iors. Depending on the level of correlations between the random phases three 
cases have been found so far [239,259,298]. 



The first of these was introduced by Cho and Fisher [210] with the intention 
of modeling the two dimensional random bond Ising model (RBIM) (Chap. 8) 
which possesses a fermion representation with the symmetries of class D. Sub- 
sequently, it has been noted [211,299] that a precise mapping of the Ising model 
leads to a second version of the model which accordingly is denoted as random 
bond Ising model. In both of the models, scattering phases n appear in cor- 
related pairs. A third model [239,277] denoted by 0(1), can be established if 
one assumes the scattering phases to be independent random variables. Each 
model has two parameters: the disorder concentration, p (0 < p < 1), and 
a tunneling amplitude sin a (0 < a < n/2) which controls the value of the 
thermal Hall conductance at short distances. 



A phase diagram for the Cho-Fisher model in the (p, a)-plane is shown in 
Fig. 43. It contains three regions: a metallic phase, a localized phase, and a 
thermal quantum Hall phase. While the random bond Ising model has been 
shown to have no metallic phase [211] the 0(1) model has been argued to 
have no insulating and thus no quantum Hall phase [268] . Only the Cho-Fisher 
model has all three phases. It features three potentially different critical points: 
a quantum Hall-type transition, an insulator-to-metal transition, and a multi- 
critical point at which all three phases meet. This phase diagram has the form 
proposed generically for class D [264]. The phase diagram is very complex. 
The corresponding critical parameters remain to be determined quantitatively. 
The network model reviewed here [239] will certainly be central in achieving 
a better understanding of this phase diagram in future investigations. An 
attempt of understanding the implications of the all-orders beta functions has 
been made in [300]. 
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Fig. 43. The phase diagram in the (p, a)-plane of the model of Cho and Fisher 
[210] obtained by numerical calculations [239]. Here, crj y is the thermal Hall con- 
ductance parameter defined in Eq. (310), p the disorder parameter, and sin 2 (a) is 
the tunneling probability (Figure taken from [239]). 

9.4 Chiral Network Models 

In the following, we shortly summarize the network models of the chiral sym- 
metry classes introduced recently [260]. The network models presented here 
are connected in two distinct ways to the models without chiral symmetry 
introduced in the previous sections. Class AIII is connected with the U(l) 
network model, while CII and BDI belong to the SU(2) and 0(1) network 
models, respectively. They describe the plateau transitions in dirty super- 
conductors, as briefly explained above. In each case, the model with chiral 
symmetry is constructed by coupling two copies of the partner model. After 
a transformation, this coupling is equivalent to introducing absorption and 
coherent amplification in the original models. This equivalence parallels the 
established link [301] between the chiral symmetry classes and non-Hermitian 
random operators. 

We formulate the network models for each of the chiral symmetry classes 
by treating in some detail the chiral unitary class (AIII), and outlining the 
equivalent steps for the chiral orthogonal (BDI) and chiral symplectic (CII) 
classes [260]. The strategy is to construct the two dimensional internal space 
associated with chiral symmetry using two related copies of a network model 
without that symmetry. 

The symmetry may be discussed in terms of a Hamiltonian H, a scattering ma- 
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trix S, or a transfer matrix T. We consider systems which conserve probability 
density. The Hamiltonian is Hermitian and the scattering matrix is unitary. 
The equivalent condition for the transfer matrix involves the operator of the 
current J, and is imposed by current conservation as in Eq.(91), 

H f = H, S f = S-\ T^JT+J. (342) 



Chiral symmetry is implemented on a two dimensional internal space in which 
the Pauli operator a x operates. For the Hamiltonian this is equivalent to 

a x Ha x = -H . (343) 



By assuming the scattering matrix to have the symmetry of e * H , one has 
a x Sa x = S- 1 . (344) 



Systems in class AIII have no other relevant discrete symmetry. Those in 
classes BDI and CII are also invariant under time-reversal, in the presence 
and absence of spin rotation symmetry, respectively. 

Class AIII 



We recall first the essential features of the U(l) network model for the integer 
quantum Hall plateau transition (Chap. 3), and put them into the present 
context. The forms of the transfer matrix and of the evolution operator follow 
from the properties of the elementary building units, sketched in Fig. 44. A 
wave function in this model takes complex values ipi on the links / of the 
lattice illustrated by the full lines of Fig. 45. A particle acquires a phase (pi 
on traversing link /, so that in a stationary state amplitudes at either end are 
related by 

tj}' = e^tj}. (345) 



In a similar way, stationary state amplitudes on the four links which meet at 
a node are related by a 2 x 2 transfer matrix 



(cosh a sinh a 
sinh a cosh a 




(346) 



where a is real and all phase factors are associated with links. This equation 
may be rewritten in terms of a scattering matrix as 
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cos a — sm a 



sin a cos a 




(347) 



with sin a = tanha. The transfer matrix that results from assembling these 
units is described in detail in [31], and the time evolution operator in [121,189] 
(cf. Sect. 3). 




Fig. 44. The elements of the network model: a link (left) and a node (right) [Figure 
taken from [260]]. 



Introducing a two dimensional internal space associated with the chiral sym- 
metry results in a doubling of the number of wave function components. Thus, 
starting from two copies of the U(l) model the link amplitudes become two- 
component complex numbers ipi. Instead of Eq. (345), the scattering properties 
of a link now are characterized by a 2 x 2 transfer matrix T link , with 

V>' = T link -0. (348) 



Requiring 

O'xTlink^x = Tii nk , (349) 



the transfer matrix must have the form 



^ cosh b sinh b ^ 



- link 



sinh b cosh b 



(350) 



where is a real phase and b is a real hyperbolic angle. It remains to discuss 
scattering at nodes of the doubled system. We replace Eq. (346) by 




cosh a sinh a 
sinh a cosh a 



(351) 
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where I2 is here the unit matrix in the two-component space introduced on 
links. This choice is the most general one which is compatible with chiral 
symmetry, since scattering within the two-component space may be included 
in the link transfer matrices Ti ink . 

By combining these elements into a two dimensional system, we arrive at the 
model schematically shown in Fig. 45. The transfer matrix for the system as a 
whole acts in the [1, 1] (or [1,1]) direction and may be written as a product of 
factors relating amplitudes on successive slices of the system. Alternate factors 
in the product represent links and nodes, and consist respectively of repeated 
versions of the 2x2 and 4x4 blocks appearing in Eqs. (350) and (351). Disorder 
is introduced by assuming the phase in Eq. (350) to be independently and 
uniformly distributed on the links. The parameter a, characterizing scattering 
at nodes, is assumed to be non-random (Chap. 3). 

Two possibilities have been considered for setting the value of the coupling 
between the chiral subspaces. They are parameterized in terms of either the 
hyperbolic angle in Eq. (350) or the compact angle (3 related to b by sin j3 = 
— tanhfr. Either (5 is assumed to be uniformly distributed, or b is assumed 
to be normally distributed with a variance g. A disorder parameter 7 can be 
defined by 

sin 7 = tanh ^/g . (352) 



It is guaranteed that the system is statistically invariant under 7r/2 rotations 
of the lattice. This constrains the node parameter a, exactly as in the U(l) 
model: nodes lie on two distinct sub-lattices, and the node parameter a on 
one sub-lattice is related to the parameter a' on the other sub-lattice by the 
duality relation (cf. Eq. (95)) 

sinh a sinh a' = 1 . (353) 



As before (Sect. 4.3), the system can alternatively be described using a time 
evolution operator instead of a transfer matrix. In order to specify this unitary 
operator, which has the symmetry of a scattering matrix, it is convenient 
first to consider the limit b = 0. Then, the two copies of the network model 
are uncoupled. One denotes the time evolution operator for one copy by U. 
Since from Eq. (350) link phases are the same in both copies but propagation 
directions are opposite, the evolution operator for the other copy is U*. The 
dimension of U is equal to the number of links in the system. It is useful 
to define a diagonal matrix of the same dimension, with angles f3i for each 
link / as diagonal entries: we use (5 to denote this matrix. The time-evolution 
operator for the system with chiral symmetry can then be written 
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Fig. 45. The structure of the AIII network model. Full and dashed lines indicate the 
two U(l) models from which the system is constructed, with nodes located at the 
vertices of the two lattices. Scattering that couples the sub-systems is represented 
schematically by boxes (upper left plaquette), as parameterized by Tn n k, Eq. (350). 
In one step of time evolution, flux propagates between successive points marked with 
filled circles, in the directions indicated by the arrows [Figure taken from [260]]. 

, Ucos/3 -Usm/mA 
U T =| H \ \. (354) 

sin (3 cos (3 IF I 

One can straightforwardly verify that o x \J^a x = U^ 1 and that = U^ 1 - 

For this model, two regimes of behaviors have been identified [260]. For small 
7, there is a localized phase. For 7 close to 7r/2, one finds a critical phase. 
By identifying the value of 7 which divides the two regimes and studying its 
dependence on a, one arrives at the phase diagram shown in Fig. 46. The 
line a = 7r/4 corresponds to the critical-metal states studied by Gade and 
Wegner [253]. On the line without hopping between the sub-lattices, 7 has a 
single critical point at a = 7r/4 which is in the universality class of the Integer 
Quantum Hall Transition. 



Classes BDI and CII 



A model in class BDI can be obtained from one in class AIII simply by im- 
posing time reversal invariance as an additional symmetry. This condition is 
conventionally written in the form H* = H, but for a discussion based on 
scattering matrices it is more convenient to transform 

H -> QHQ 1 with Q = e i7RT * /4 . (355) 
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Fig. 46. Phase diagram of the AIII network model in the (a,7)-plane [Figure taken 
from [260]]. 

This transformation leaves the chiral symmetry relation a x tla x = — H intact. 
In the transformed basis one has 



H* = -H S = S* , T* = T 



(356) 



To ensure a real time evolution operator, we restrict the link phases to the 
values and n. Choosing these values randomly, the BDI model consists of 
two coupled copies of the class D models reviewed in the previous section. 
Alternatively, one could assume = on all links, and introduce disorder 
only via the chiral couplings (5. 

For class CII, Kramers degeneracy is a defining feature of the class and re- 
quires the introduction of an additional two dimensional space arising from 
the spin. The time reversal operation includes reversal of spin direction. Defin- 
ing C = i(j y , where a y is a Pauli matrix acting in the additional space, it is 
conventionally written in the form 



CH*C 1 = H. 



(357) 



As for the class BDI, it is again convenient to transform according to Eq. (355). 
This gives 



CWC L = -H, CS*C~ L = S, CT*C- L = T 



(358) 



as equivalent expressions of time-reversal invariance. Applying these ideas to a 
network model, four channels propagate along a single link which, generalizing 
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Eq. (350), is transferred according to a 4 x 4 transfer matrix Ti ink with the 
generic form 



where v is an SU(2) matrix and b a real hyperbolic angle. Adopting this 
form, the time evolution operator for class CII has the structure given in 
Eq. (354), but with U representing a network model of class C, as reviewed 
in the preceding section. The links of a such a model of class C carry two co- 
propagating channels, coming from two spin components, and the evolution 
operator satisfies 



10 Supersymmetry and Localization 

In the present chapter, we comment on recent developments which open new 
perspectives for better understanding of the field theoretical models that have 
been developed during the past two decades for describing the quantum Hall 
phase transition. This development is closely related to a mapping to a chain 
of superspins starting from the random network model. 

While the numerical and experimental evidences for the universality of the 
quantum Hall transition appear to be striking, the quantum Hall critical point 
has so far successfully escaped any attempts to be placed into a classification 
of two dimensional phase transitions, based on the conformal invariance at the 
critical point. Recently, it has been suggested that the quantum Hall critical 
point may belong to a new class of critical points being described by a su- 
persymmetric conformal field theory [60,61]. Much progress has been achieved 
towards the analytical derivation and characterization of the quantum Hall 
transition. Nevertheless, an analytical calculation of its critical exponents is 
still missing. 

Soon after the discovery of the Quantum Hall Effect a field theory has been 
derived from the microscopic Hamiltonian of the random Landau model with 
short-range disorder, the correlation length of the disorder potential, £ c , be- 
ing much smaller than the magnetic length £ B [302,303,304,305]. It has been 
shown to have two coupling parameters a\ x and a® y , the longitudinal and Hall 
conductance as defined on small length scales of the order of the elastic mean 
free path I. This field theory is based on the theory of localization of electrons 
in weakly disordered systems. 




(359) 



CWC' 1 = U. 



(360) 
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Subsequently, it has been shown rigorously that this field theory is indeed 
critical at half integer Hall conductance parameters cr® y [302,306,307], and 
that it has a spectral gap to fluctuations at other values of a® . This indicates 
the localization of the electron eigenstates of the random Landau model in 
the tails of the Landau bands [308]. Since the longitudinal conductance at 
the critical point a* is known to be smaller than 1 (in units of e 2 /h), the 
critical point is located, quite unfortunately, in the strong coupling limit of 
the field theory. Thus, it is outside of the validity of the available analytical 
methods which can be used to extract quantitative information on the critical 
exponents. 

Recently, an anisotropic version of the Chalker-Coddington model has been 
mapped directly onto the Hamiltonian of a chain of antiferromagnetically in- 
teracting superspins [37,192,309,310]. It had been shown before that the non- 
linear sigma model for short-range disorder at the critical point, a xy = e 2 /2h, 
and in the strong coupling limit, can also be mapped onto the Hamiltonian 
of a chain of antiferromagnetically interacting superspins [311]. This provides 
strong analytical support for the notion of universality of the quantum Hall 
transition as sketched schematically in Fig. 47. Thereby, the problem of the 
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Fig. 47. Sketch of the development of the analytical proofs for the equivalence of the 
critical theory for the random Landau model and the Chalker-Coddington network 
model. 
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quantum Hall transition has been transformed to the task of finding the ground 
state and the dispersion of the excitations of the chain of antiferromagnetic 
superspins [37,310,312]. 

The model of antiferromagnetically interacting superspins has been shown to 
be critical [310]. Numerically, the critical exponent v has been obtained from 
a finite length scaling of superspin chains, by means of the density matrix 
renormalization group method, and found to be v = 2.4(4) [37]. 

So far, no analytical information could be directly obtained on the critical pa- 
rameters, the localization exponent u, the MacKinnon-Kramer parameter A c 
of the scaling function, and the multi-critical exponents D(q). However, start- 
ing from the model of a superspin chain, a class of supersymmetric conformal 
field theories has been suggested to be relevant for the quantum Hall tran- 
sition, which ultimately should yield the critical parameters of the quantum 
Hall transition [60,61]. In yet another development, the Bethe-Ansatz method 
has been applied to the Hamiltonian of superspin chains. Still, it is not yet 
clear if that model is integrable and a closed solution for its ground state and 
its excitations can be obtained in this way [313,314]. 

In the next section we will provide a derivation of the Hamiltonian of an 
antiferromagnetic superspin chain starting from an anisotropic version of the 
Chalker-Coddington model, following closely the original derivations by Lee 
[192], Zirnbauer [309], and by Kondev and Marston [37]. In the second section, 
the derivation of Pruisken's non-linear sigma model in supersymmetric form is 
summarized [302,303,304,305]. Finally, the connection between that model and 
the superspin chain is discussed following the work of Zirnbauer [60,309,311]. 

10.1 From the Network Model to the Antiferromagnetic Superspin Chain 

We start from a modification of the network model, recently proposed by Lee 
[192], where the random scattering matrices are replaced by a Hamiltonian 
which describes a quasi-one dimensional array of counter-propagating edge 
states in the presence of disorder (Fig. 48). The ^-coordinate is discretized, 
the integer index n enumerates the discrete vertices of the network model. 
The fact that the direction of propagation alternates between sites n will be 
shown in the following to result in an antiferromagnetic interaction between 
superspins. A one dimensional array of edge states, which propagate all in 
the same direction, has accordingly been mapped to a ferromagnetic super 
spin chain [312]. This is an effective model of the surface states of a layered 
quantum Hall system, a so-called chiral metal (Chap. 11). 

In terms of electron creation and annihilation ip operators the Hamiltonian 
may be written as 
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1 2 3 4 n 



Fig. 48. The network model with counter-propagating chiral fermions with random 
on-site potential w njy describing the propagation along edge states around Hall 
droplets with random tunneling amplitudes t rhy between them (dashed). 



where the sum is over the discrete vertices of the network model. This Hamilto- 
nian reflects the chiral nature and the linear dispersion of the edge states, with 
alternating propagation forward and backward in the y-direction with Fermi 
velocities rkv-p. The complex tunneling amplitudes t n>y between the edge states 
are assumed to have random phases. These represent the random Aharonov- 
Bohm phases of the electrons, accounting for the fact that in a network of 
edge states, closed orbits vary randomly in size and thereby encircle randomly 
varying magnetic fluxes. There are also random on-site potentials w n>y . These 
random terms are assumed to be Gaussian distributed with zero means, and 
the variances 



with J n = J[l + (— l) n R\. The staggered modulation in the parameter J n allows 
to trace the differences in tunneling between counter-propagating electrons in 
adjacent columns of plaquettes as depicted in Fig. 48. At R — —1, there 
are disconnected pairs of counter-propagating edge states, such that electrons 
circulate only clockwise around the (+)-plaquettes. In the opposite case, R = 
1, electrons circulate counterclockwise only around the (— )-plaquettes. At 
R = 0, there is tunneling between all sites, and critical quantum percolation is 
expected, corresponding to the transition between the quantum Hall plateaus. 

Critical behavior consistent with this scenario was found from the disorder av- 



(tn,y1pn t y1pn+l,y + tn,yYn+\,y^n,y)\ j 



(361) 



(wn, y «Vy) = 2*7 6 n y 5(y- y') , 




(362) 
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eraged two-particle Green function describing transport corresponding to the 
Hamiltonian Eq. (361) by Wang and Lee [38]. From a Monte-Carlo treatment 
of the replicated Hamiltonian they obtained a correlation length exponent 
v = 2.33(3), in good agreement with other numerical simulations and exper- 
iments. It can furthermore be shown that the Hamiltonian Eq. (361) can be 
identified in the continuum limit with the Hamiltonian of two dimensional 
Dirac fermions with random mass and random vector fields [192], which is 
in the universality class of the integer quantum Hall transition, as discussed 
previously in the chapter 7. 

It is well known [75,76,77,78,79], that in weakly disordered systems, k-pl ^> 1, 
one needs to go beyond lowest order perturbation theory in the disorder po- 
tential in order to describe quantum localization. One also has to take full ad- 
vantage of the symmetries occuring in the calculation of correlation functions 
of disordered systems. This can be traced back to the fact that the disorder 
averaged electron wave function amplitude (ip(r,t)) decays on length scales 
of the order of since the random scattering phase shifts associated with the 
scattering at the impurities are averaged out. This destroys the information on 
quantum coherent multiple scattering, and thus on quantum localization. In 
order to describe quantum localization, one needs to consider higher moments 
of the wave function amplitudes such as the impurity averaged evolution of 
the electron density n{r,t) = {\ip{r,t)\ 2 ). 

The time evolution of the electron amplitudes can be written in terms of the 
retarded propagator G R , 

^(r, t) = J dr'G R (r, t; r', t')^(r', t'), (363) 

with t > t'. The electron density becomes 

(n(r, t)) = J dr' J dr"T(r, t; r', r", t')^(r', t')^*(r", t') , (364) 

where the quantum diffusion propagator is given by 

T(r,t;r',r",t') = (G R (r,t;r' ,t')G A (r" ,t';r,t)). (365) 

After performing a Fourier transformation from time t — t' to energy E, a 
non-perturbative averaging of products of retarded and advanced propaga- 
tors, (G K (E)G A (E')) is needed in order to obtain information on quantum 
localization. 

In a useful analogy to the study of spin systems, the field theoretical approach 
contracts the information on localization into a theory of Goldstone modes Q, 
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arising from the global symmetry of rotations between the functional integral 
representation of the retarded propagator G R ("spin up") and the advanced 
propagator G A ("spin down"). The field theory can either be formulated by 
means of the replica trick, where the N replicas are represented either by iV 
fermionic or N bosonic fields, yielding a bounded or unbounded symmetric 
space, respectively, on which the modes Q are defined [302]. 

Because of the necessity and the difficulty to perform the delicate limit N — > 
at the end of the calculation in the replica formulation, a more rigorous su- 
persymmetric field theory has been formulated. This technique represents the 
product of Green functions G K (E)G A (E') by functional integrals over two 
fermionic and bosonic field components, composing a supersymmetric field 
vector ip. The supersymmetric representation enables one to perform the av- 
eraging over the disorder potential as a simple Gaussian integral [303,304,305]. 

Thus, in order to study the localization-delocalization transition in the net- 
work model as described by the random Hamiltonian, Eq. (361), we consider 
the average of the quantum diffusion propagator 

K(l,2) = (Gi(l,2)G£,(2,l)>, (366) 



where 



^ A (l,2) = (11-^-^-12), (367) 



are the retarded and advanced Green functions. The parameter rj is a positive 
infinitesimal number. We have introduced here the usual shorthand notation 
for the coordinates with 1 := (ni,yi) and 2 := (712,3/2)- For non-zero disorder 
parameter R we expect the quantum diffusion propagator to decay exponen- 
tially 

K(l,2) ~ e~ 2ri2/ t , (368) 



where ryi is the distance between points 1 and 2. This defines the localization 
length £ I, reflecting the finite extent of the electron wave functions near 
energy E. At the critical point R — 0, the localization length is expected to 
diverge like 

f ~ R- v . (369) 
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By introducing a pair of complex scalar fields 

<Kn,y)= | Mn,V) | , (370) 
(f)-(n,y) 



with + and — denoting the retarded and the advanced sector, respectively, 
the advanced and retarded Green functions, Eq. (366), can be rewritten. One 
obtains for the two-particle propagator 

(K(l,2)) = 1 J D[0]D[0*]0 + (l)0* + (2)0_(2)0*_(l)e- 5 ^, (371) 



with the action 



S[<j), 0*]=EE [ ia {Kin, y) (-E a + (-l) n tv F d y + w n>y ) <j> a {n, y) 

n,y a=± 

[t n ,v<f>* a (n, y)<p a (n + l,y) + t^ y <p* a {n + l,y)4>Jn, y)] } 
+r](f) a (n,y)(f)* a (n,y) 



(372) 



Here, we have defined E a=+ _ = E, E' . The choice of the signs guarantees that 
the Gaussian integrals are convergent for r\ > 0. The normalization factor 



/ mm 



*] e -s[4>,r] 



(373) 



is an inverse spectral determinant. Therefore, its inverse can be lifted to the 
numerator by introducing another Gaussian integral, this time over Grass- 
mann fields [315] 



D[ X ]D[x],e- 5 ^ 



(374) 



with the anti-commuting fields 

x{n,y) 



1 x+(n,y)^ 



\X-(n,y) J 



(375) 



that are the supersymmetric partners of the scalar bosonic fields. Using Eqs. (374) 
and (371) the correlation function, Eq. (366), can be written as a combined 
Gaussian integral over scalar and Grassmann fields 
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(K(l,2))= J D[0]D[0*]D[ X ]D[x]0 + (l)0;(2)0_(2)C(l) 

Xe -(S[W1+%,xl). (376) 

Now the average over the disorder can be performed as a simple Gaussian 
integral. 

Before considering the full Hamiltonian, let us look first at the correlation 
function of a single chiral edge state n, which is expected to behave metallic 
for any disorder strength U since its Hamiltonian contains no backscattering, 
J = 0. 

Performing the ensemble averaging as the Gaussian integral over the random 
potential w(n,y) leads to the functional integral 

D[-0, $ exp f dx \i>(Av F d x - rj + - U{ipA^) 2 ] , (377) 



We have here introduced a four-component superspinor ip(x) with components 
ipx where X = KB (retarded Boson), X = RF (retarded fermion), X = AB 
(advanced Boson), and X = AF (advanced fermion). The 4x4 supersymmetric 
A is defined by 

V 5 = (ipRa^Rv ~ i>Aa^Aa) ■ (378) 

cr=B,F 



Since the energies in the retarded and advanced sectors are different, we have 
defined u = E - E' ^ 0. 

One can decouple the interaction term (ipAip) 2 by introducing a Hubbard- 
Stratonovitch field Q coupling to ipipA and then integrate out ip and ip. The 
resulting effective action for Q is 



S[Q] = J 



dx STr 



-ttQ + ln ^ + ^ A + ^ 

4«f V V F J 



(379) 



The next step is to simplify the Q-field functional integral by means of the 
saddle-point approximation. As a result, Q gets restricted to the non-linear 
space Q = gAg^ 1 . It belongs to the supersymmetric space U(1,1|2)/(U(1|1) x 
U(l|l)) which is sketched schematically in Fig. 49. Here the notation is reminis- 
cent of the notation for a 2-dimensional hyperbolic space U(l, 1) and the group 
of 2 by 2 unitary matrices U(2). Accordingly, the supermatrices A e U(l, 1 | 2) 
have the form 
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(380) 



where A BB e U(l, 1) and A FF e U(2) and A FB , A BF are parameterized by 
Grassmann variables, and transform between U(l, 1) and U(2). 



This step, being in general only an approximation, here becomes exact in the 
limit U — > oo. By expanding In [Ag~ 1 (d x + A)g + U/vp] to linear order in 
g^ 1 (d x + A)g, one obtains the action of the Wess-Zumino functional Z WZ [A], 



[u] =JDQ 



exp <^ / da; STr 



^9 1 (v F d x + iujA)g 



(381) 



Higher orders are suppressed by powers of v\jL x l) with L x the system size. 

By interpreting the integral over x as an integral over time, one can rewrite 
this in terms of a Hamiltonian 

Z wz [uj] = STre- LyH ^\ (382) 



where 



H 9 (lu) =icu STr --Ag. (383) 

2 g 



This is identical to the super Hamiltonian which one obtains from the uni- 
tary ensemble of random matrices [303], typical for a random metallic system 
with extended states. Thus, as expected, we recover the result that a single 
disordered chiral line supports only metallic states. 

Let us next discuss the case of many counter-propagating chiral modes that 
are coupled by hopping matrix elements between neighboring modes, with 
variance J. This corresponds to the anisotropic Chalker-Coddington model. 
The Gaussian random hopping matrix elements give rise to an additional term 
in the Lagrangian, 

C C + 2 J (^nA^n+l) (^n+lA^n) 
n 

= C + 2J^2 STr (i) n i) n A) (VwiVWi A) . (384) 

n 

This is obtained by averaging Eq. (376) with Eq. (372) over t using Eq. (362). 
By using the bosonization rule v-pip^A — > Q/2 for U — > oo, the additional term 
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Fig. 49. Schematic view of the supersymmetric space U(l, 1|2)/(U(1|1) x U(l|l)) 
on which the Hamiltonian of supersymmetric operators Q is defined. The point at 
which the compact sphere (corresponding to that part of Q which arises from de- 
coupling the fermionic ^p-term) and the hyperbolic (arising from the decoupling of 
the bosonic ^-term) meet, is the classical point which yields the classical correla- 
tion function. Diagrammatic expansion around that point would miss the curvature, 
and the non-perturbative integral over the whole supersymmetric space is needed to 
describe quantum localization. The compact sphere is parameterized by $ 6 [0, 2ir) 
and A G [—1,1], and the non-compact hyperbolic is parameterized by Ai 6 [— l,oo] 
and 9 £ [0, 2tt). The unitary rotations between the sphere and the hyperbolic are 
parameterized by Grassmann variables £ as indicated by dotted lines. 



can be cast into (J/2t>p) J2 n STr(Q n Qn+i) • The condition for the validity of 
this step is U ^> J. As a result one obtains the Q-field action 

S[Q] = I dx^STr 

J n 



- l ) n ^9n ld x9n + ^QnQn+1 



(385) 



This is the action [312,316] of the coherent-state path integral for a quantum 
superspin Hamiltonian. The resulting effective supersymmetric (SUSY) Hamil- 
tonian describes interacting spin- up and spin-down fermions c a and bosons b a . 
The two spin species formally correspond to the retarded and advanced Green 
functions introduced above [37,310,317], 



L-2 



3=0 



a=l 



a=5 



16 L-1 

+\ H (-iy E 9a Sf +1 ] + V £ [Sj + S] + + S?] . (386) 



a=9 



3=0 
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The signs g a are given by 



,1 for a = 1,2,10,12,14,16 
9a = { (387) 
-1 for a = 3,..., 9, 11, 13,15 . 

Equation (386) contains 16 spin operators, the components of a 4 x 4 superspin 
matrix, 



5 1 ■= bfa + 1/2 S 5 : = c{c T - 1/2 S 9 := c\b l S 13 := b iC] 

5 2 := fcft + 1/2 S G := cjq - 1/2 S 10 := c}& T S 14 := 6 T q 

(388) 

5 3 :=fe}fo{ 5 7 :=c}c{ 5 11 := b\ Cjr S 15 := b\c\ 

S A := bfa S s := qc T S 12 := &{c T S w := b\c\ . 



The boson- valued operators S l ,...,S 8 constitute the symmetric sector of 
the Hamiltonian. The fermion- valued operators S 9 , . . . , S* 16 form the anti- 
symmetric sector. Despite that H is non-Hermitian, it only has real-valued 
eigenvalues. For A# = one has only fermionic spin operators. One does arrive 
at this Fermion model directly from the Chalker-Coddington network model 
when calculating the disorder averaged auto correlation function of spectral 
determinants [60,307]. Another derivation via the Landauer conductance for- 
mula was recently reported in Ref. [318]. This yields the antiferromagnetic 
spin-1/2 Heisenberg chain, which is critical [319]. The crossover to the in- 
tegrable superspin chain of the spin Quantum Hall Effect [265,281], which is 
critical in the universality class of the classical percolation model, can be stud- 
ied by freezing out 8 of the 16 spin components [320] . This raises the question 
whether or not the critical parameters vary continuously, as for the crossover 
from the antiferromagnetic Heisenberg spin chain to the XY model [321,319]. 

The Hamiltonian commutes with four fermion- valued supersymmetry genera- 
tors, 

[H, Q la ] = [H, Q 2a ] = (389) 



with the supersymmetric charges 
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Qia ■= 2 (-l)^ CT c j(7 + cj- ff 6 jff 



(390) 



One can see that the supersymmetric Hamiltonian must have a unique, zero- 
energy ground state. All excited states appear in quartets or larger multiples of 
4, half with odd total fermion content. These cancel out in the partition func- 
tion by virtue of the super-trace, yielding the correct value of 1 that we have 
encountered above when introducing the functional integral over Grassmann 
variables to cancel the normalization factor of the bosonic integral 



STre"^ := T^-l^e"^ = 1 . (391) 



Here, N c is the total number of fermions. The ground state of this supersym- 
metric non-Hermitian antiferromagnetic Hamiltonian is very complicated. The 
ground state of the Hermitian supersymmetric ferromagnet that describes a 
chiral metal with all edge states propagating in the same direction is simply 
the vacuum state [312]. 

The Lieb-Schultz-Mattis theorem proves that for half-odd-integer spin antifer- 

— * — # 

romagnets on a periodic chain of length L sites (Sl '■= So) either the ground 
state is degenerate or there are gapless spin excitations in the thermodynamic 
limit L — > oo [241]. In the supersymmetric problem, Marston and Tsai were 
able to make even a stronger statement [310,317], because its ground state is 
unique by supersymmetry. For rj > they showed that low-energy excitations 
are gapless in the thermodynamic limit. This proves that the antiferromagnetic 
superspin chain is critical. 

In the following section the field theory of the quantum Hall critical point 
is approached starting from a model of electrons in uncorrelated disorder in 
a strong magnetic field. The resulting supersymmetric field theory, with the 
action of a non-linear sigma model was used as a basis of the two parameter 
scaling of the Quantum Hall Effect. It has recently been shown that it can in 
a long wave length limit be mapped onto the Hamiltonian of the superspin 
chain, Eq. (386). This, together with the derivation of the spin chain from the 
opposite limit of a model with long-range potential as sketched above, provides 
very strong support to the surmise that the Chalker-Coddington model is a 
good model for the universal quantum Hall transition 
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10.2 From the Landau Model to Pruisken's Non-linear Sigma Model 



The Hamiltonian of noninteracting electrons in a magnetic field in the presence 
of uncorrelated disorder is (Chap. 1) 

H = ^- (p + eAf + V(r) + V (r) . (392) 
2m 



Here, V(r) is assumed to be a Gaussian distributed random function with a 
distribution 



P(\V])=exp 



dr dr' . 



r')V(r)V(r') 



(393) 



where f2 is the volume of the system. Impurity averaging is thus given by 
(■■■)v = SU T d\y]P(W])- ■ We assume 

J(r-r') = nA-6(r-r'), (394) 

T 



for uncorrelated impurities, where 1/r is the elastic scattering rate and A = 
l/(pQ) the mean level spacing of the mesoscopic sample with volume Q. It 
is related to the variance of the disorder potential V(r) according to W 2 = 
Ah/2irr. The function V (r) is the electrostatic confinement potential defining 
the width of the wire L. 

In order to describe localization, we consider again the correlation function 
K(\, 2) studied in the previous chapter for the Chalker-Coddington model, but 
using this time the Hamiltonian Eq. (392). Formulating the supersymmetric 
field theory by representing the product of Green functions G R (E)G (£") by 
functional integrals over two fermionic and bosonic field components, compos- 
ing a supersymmetric field vector ip, as in the previous section, the averaging 
over the disorder potential can again be performed as a simple Gaussian in- 
tegral [303,304,305]. As a result of the averaging one obtains a locally inter- 
acting theory of the fields ip containing an interaction term oc ip 4 , where the 
interaction strength is proportional to the variance of the disorder potential 
W 2 . This term is next decoupled by introducing another Gaussian integral 
over Q— matrices. Clearly, the field Q should not be a scalar, otherwise we 
would simply reintroduce the Gaussian integral over the random potential V. 
Rather, in order to be able to describe the physics of localization, the field Q 
should capture the full symmetry of the functional integral representation of 
the correlation function. Therefore, the Gaussian integral is chosen to be over 
a 4x4 matrix Q which itself is an element of the symmetric space defined by 
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the matrices A that leave the functional integral invariant under the trans- 
formation ip — > Aip. In the supersymmetric formulation, this matrix consists 
of two blocks of 2x2 matrices whose parameter space consists of a compact 
(bounded) and a non-compact (unbounded) sector, Eq. (380) (Fig. 49). The 
off-diagonal blocks, so to say the rotations between the compact and the non- 
compact sector, are then found to be parameterized by Grassmann (fermionic) 
variables. Now, the spatial variations of Q are governed by the action 

jrh f r\r 1 f 

S[Q] = ^J J2 Tt Q 2 ( r ) + 2 J dr < r ' Tr lnG (^P)|r> ( 395 ) 
where 

G-^r,p) = ^±^-^l-V (r) + ^Qm. (396) 

The 4x4 matrix A3 is the diagonal Pauli matrix in the sub-basis of the retarded 
and advanced propagators, r\ > and uj = E — E' break the symmetry between 
the retarded and advanced sector. 

It turns out that the physics of diffusion and localization, which arises on 
length scales much larger than the elastic mean free path /, is governed by the 
action of the long wavelength modes of Q. Thus, one can simplify and proceed 
with the analysis by expanding around a homogeneous solution of the saddle 
point equation, 5S = 0. For uj = 0, this is 

-1 

\r). (397) 

This is solved by Q = A 3 P, which corresponds to the self consistent Born 
approximation for the self energy of the impurity averaged Green function. 
At uj = 0, rotations U which leave the action invariant yield the complete 
manifold of saddle point solutions as Q = tJA 3 PU, where UU = 1 where U 
is the supersymmetric Hermitian conjugate. 

The modes which leave A 3 invariant can be factorized out, leaving the saddle 
point solutions in this supersymmetric theory to be elements of the semi- 
simple supersymmetric space G1(2|2)/(G1(1|1) x Gl(l]l)) [304,305]. Here the 
notation is chosen in analogy to Gl(n), denoting invertible n by n matrices, 
that is detA ^ for all A e Gl(n). Thus, Gl(n | n) denote invertible 2n by 2n 
supermatrices 

/ A FF A FB\ 

A = , (398) 

\A BF A BB 



Q 



1 1 1 

— (r 

TTO 



E-H -V (r) + —Q 
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with 



SdetA = det(A** - A A A )detA ^ , (399) 



where Sdet is the superdeterminant. See [244] for a more detailed definition. 
In the semi-simple space the subgroup Gl(l|l) x Gl(l|l) consisting of matrices 



(h u 

h = 

^ h 22 

with h n ,h 22 e Gl(l|l), are factorized out. 



(400) 



In addition to these gapless modes there are massive longitudinal modes with 
Q 2 ^ 1 which only change the short distance physics, and not the physics of 
localization. They can be integrated out [302,303]. Thus, the partition function 
reduces to a functional integral over the transverse modes U. 

The action at finite frequency oo and slow spatial fluctuations of Q around the 
saddle point solution can be found by an expansion of the action S. Inserting 
Q = UA3PU into Eq. (395) and performing the cyclic permutation of U under 
the trace Tr allows a simple expansion to first order in the energy difference u 
and to second order in the commutator U[iJ , U] [302]. The first order term 
in XJ[H Q , tj] is proportional to the local current. It is found to be finite only at 
the edge of the wire in a strong magnetic field, due to the chiral edge currents. 
It can be rewritten as 

1 r a 0II (r) 

S xy ii = --J dxdy -^-± STr(Qd x Qd y Q) , (401) 



where the pre-factor is the non-dissipative term in the Hall conductivity in 
self-consistent Born approximation [303] 

1 he 1 

W = — — (r\(xn y - yn x )lmG*(E)\r), (402) 



where 7r = (h/i)V + eA. This field theory has now the advantage that one can 
treat the physics on different length scales separately: the physics of diffusion 
and localization is governed by the action of spatial variations of U on length 
scales larger than the mean free path /. That is why this field theory is often 
called diffusive non-linear sigma model. 

The physics on smaller length scales is included in the coupling parameters of 
the theory, which is identified in the above derivation as correlation functions 
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of Green functions in self consistent Born approximation, being related to the 
conductivity by the Kubo-Greenwood formula, 



he 2 

r) = — 2 (rKG*(E)npGt(E + co)\r). (403) 



The remaining averaged correlators, involve products G^(E)G^(E + uo) and 
Gq(E)G^(E + uj) and are therefore by a factor h/rE smaller than the con- 
ductivity, and can be disregarded for small disorder. Using the Kubo formula 
(403), the action of Q simplifies to 



(V,Q(r)) 2 

-A J dra y{uj = Q) r) Tr [Q dxQdyQ ] , (404) 

where a xy (u = 0, r) = o xy {oj — 0,r) + cr^(u; = 0, r) and al y (u = 0) is the 
dissipative part of the Hall conductivity in self consistent Born approximation 
Eq. (403). 

The first term in this action yields localization in two dimensional electron 
systems, signaled by the presence of a gap in the field theory. The second term 
could not be obtained by any order in perturbation theory. It is of topological 
nature. 

In two dimensions, and for a homogeneous Hall conductance it can be shown 
that this term can take only discrete purely imaginary values, 



h 

S Top = 27ri-a° xy n, (405) 



where the integers n count how often the field Q(r) is winding around its 
symmetric space as it varies spatially in two dimensions. Thus, disregarding 
the spatial variation of the coupling functions cr^r) in Eq. (405), and assum- 
ing isotropy, one finds in the two dimensional limit that there are instantons 
with nonzero topological charge q, which are identical to the skyrmions of the 
compact 0(3) non-linear sigma model, as obtained form the compact part of 
the supersymmetric non-linear sigma model [302,303]. Their action is given 
by 

F q = 2ir\q\cr xx + 2mqa xy , (406) 



where a xx = a yy and a xy are the spatially averaged conductivities. 
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Now, we can repeat the derivation of the scaling function, by integrating out 
Gaussian fluctuations around these instantons. It is clear, however that the 
contribution from instantons with q ^ is negligible, as long as a xx > e 2 /h. 
Within the validity of the 1/ ^-expansion one does not find a sizable influence 
of the topological term on the scaling function A = £/L y . Still, the tendency is 
seen that at o xy = —e 2 /2h the renormalization of the longitudinal conductance 
is slowed down and one may conclude from this observation the two parameter 
scaling diagram with a critical state of finite conductance < a xx < e 2 /h, 
which is an attractive fixed point for a xy = —e 2 /2h as shown in Fig. 50 
[302,322]. Considering the beta function of the N-th fermionic replica one 
notes [302,306] that only the last one N — 1, becomes critical and touches 
like (3 ~ — (a xx — c xx ) 2 , which is the beta function of the antiferromagnetic 
Heisenberg chain [307]. From that, one can conjecture that for N — > the 
maximum of the beta function moves to positive values, which yields the 
form drawn in Fig. 50 with another, repulsive fixed point at smaller values 
of the conductance a xx , below which there is a flow to an insulating phase 
for all values of a xy . One can also come to that conclusion noting that, when 
freezing out the bosonic degrees of freedom from the supersymmetric nonlinear 
sigma model at a xy = —1/2, one arrives again at the beta function of the 
antiferromagnetic Heisenberg chain, given above. It remains still to be shown 
explicitly that the beta function evolves indeed continuously from that beta 
function, as one adds the bosonic degrees of freedom, and to see if it results 
in the beta function drawn in Fig. 50. 




XX 



Fig. 50. The conjectured two parameter flow diagram of the integer Quantum Hall 
Effect (left), and the corresponding beta function, (3 XX = d\na xx /dlnL (right) at 
Qyx = —&xy = 1/2 (dashed line) and at a xy = (full line, units of conductivities 
e 2 /h). 

Subsequently, it has been argued by other means that this field theory is indeed 
critical at half integer Hall conductance parameters a xy [302,306,307], and that 
it has a spectral gap to fluctuations at other values of a xy . This indicates the 
localization of the electron eigenstates of the random Landau model in the 
tails of the Landau bands [308]. Since the longitudinal conductance at the 
critical point a* x is known to be smaller than 1, the critical point is located in 
the strong coupling limit of the field theory. Thus, it is outside of the validity 
of available analytical methods which can be used to extract quantitative 
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information on the critical exponents. Furthermore, it is seen explicitly that 
in order that the instanton solutions with nonzero topological charge do exist 
the system must exceed the noncritical localization length ^Dunit ~ exp(-7r 2 g 2 ), 
where g = h<j xx /e 2 , when the assumption of uniform coupling parameters 
is made. Accordingly, it has been shown that the Hamiltonian of a chain 
of antiferromagnetically interacting superspins can be derived from the non- 
linear sigma model for short-ranged disorder at the critical point a xy = —e 2 /2h 
on length scales larger than ^Dunit [311]- At criticality, a xy = —e 2 /2h, and 
at strong coupling, a xx < e 2 /h, Zirnbauer [311] showed by discretization, 
following Shankar and Read [323], that the supersymmetric nonlinear sigma 
model can be mapped on the chain of antiferromagnetic superspin chains in 
the low energy, long wavelength limit, meaning that the lattice spacing of the 
super spin chain is on the order of the noncritical localization length ^Dunit- 

So far, no analytical information has been obtained for the critical parame- 
ters, such as the localization exponent, v and the critical value A c . However, 
building on the model of a superspin chain, supersymmetric conformal field 
theories have been suggested, which ultimately are supposed to yield the crit- 
ical parameters of the quantum Hall transition [60,61,183]. 

Restricting this theory to quasi-one dimension, by assuming a finite width L y 
of the quantum Hall bar, of the order of the unitary noncritical localization 
length £ 2 Dunit = ^ exp^V 2 ^), which serves as the ultraviolet cutoff of the con- 
formal field theory, one finds that the critical value of the scaling function 
A c « 1.2, the ratio of the localization length in a quantum Hall wire and 
its finite width L y when the energy is in the center of the Landau band, see 
Section 2.3, is fixed by the eigenvalues of the Laplace-Beltrami operator of 
this supersymmetric conformal field theory [182,183,61]. This is a character- 
istic invariant of the theory, arising from the conformal symmetry, just as the 
quantization of angular momentum arises from the rotational symmetry of a 
Hamiltonian. Furthermore, based on the properties of this constrained class 
of supersymmetric conformal field theories, it has been predicted that the dis- 
tribution function of local wave function amplitudes is very broadly, namely 
log-normally, distributed. This prediction has recently been confirmed by high 
accuracy numerical calculations [137]. 

The quest to derive the critical exponent of the localization length at the 
quantum Hall transition from a critical supersymmetric theory has thus re- 
cently gained much progress, by mapping both the Chalker-Coddington model 
and the random Landau model on the Hamiltonian of the superspin chain, 
Eq. (386), which has been shown to be critical at it! = [37,310,317] as de- 
fined in Eq. (362). Thus, it has been proven to be a good starting point to 
continue in the quest for an analytical derivation of the critical exponent v 
and the multi-critical exponents at the quantum Hall critical point. 
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11 Extension to Higher Dimensions 

A natural extension of the two dimensional quantum Hall system is to con- 
sider layers stacked in parallel, where electrons can tunnel from one layer 
to another. Such layered systems can be fabricated in experiments and the 
magneto-transport properties have been investigated [324,325,326,327]. 

11.1 Double Layer Network Model 

The simplest case is a double layer system. In this case, there appear two 
energies where electron states are delocalized [29,328,329], instead of a single 
energy for delocalized states in case of a single layer system. If the tunneling 
integral t! between layers is vanishing, the positions of the delocalized states 
are both at the Landau band center, while if the random potential is the same 
for both layers, delocalized states appear at E — ±t'. The really interesting 
situation, however, is none of the above, i.e., uncorrelated disorder and finite 
tunneling. The position of the delocalized states is analyzed in detail in [329], 
while the critical behavior has been discussed numerically in [29]. 

To incorporate interlayer tunneling in the network model, we introduce scat- 
tering between links in different layers. In Fig. 51, a schematic view for double 
layer system is presented. The dashed lines are the saddle points that describe 
intra-layer scattering described by the S matrix, Eq. (87), and dotted lines 
represent tunnelings between the layer. Detailed expressions are presented in 
Section 11.3. Note that the situation is similar to the class AIII discussed in 
Section 9.4, but the discrete symmetry (Eq. (349)) is not present in this double 
layer system. 

The density of states of such a system is sketched in Fig. 52. The Landau 
band splits due to inter-layer tunneling. Whether or not the localization length 
exhibits two singularities with the same critical exponent as a single layer is 
under debate, but it is likely to be so. 

The application of the renormalization group process (Chap. 6) to two channel 
network model has been discussed in detail[330,331], where the limitation of 
the hierarchy model is clarified. 

11.2 Localization- Derealization Transition 

With the increase of the number of layers, the positions of energies where 
delocalized states appear increase. The number of positions coincides with 
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Fig. 51. Schematic view of a double layer network model. The dashed lines indicate 
intra-layer saddle points. Dotted lines represent tunneling between the layers. 

p(E) 




E 



Fig. 52. Schematic view of the density of states p(E) and the positions of delocalized 
states for a double layer quantum Hall system. Delocalized states are indicated by 
solid lines. 

the number of layers. In the limit of infinite number of layers, the delocalized 
states form a band of energies [332]. 



p(E) 




E 

Fig. 53. Schematic view of the density of states p(E) and the positions of delocalized 
states for a multi-layer system. Delocalized states are indicated by the shaded region. 
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As long as the number of layers is finite, the critical behavior of the localization- 
delocalization transition is conjectured to be the same as for a single layer. 
In the limit of infinite number of layers, the localization-delocalization tran- 
sition should become that of the three dimensional unitary system. For this 
limit, the localization length exponent has been estimated numerically in the 
multi-layer Landau model [332] and the anisotropic tight binding model [333]. 
In all cases, it has been found to be consistent with the value of the three 
dimensional unitary class [111,334], v = 1.43±0.03. This value is not far from 
4/3 obtained by mapping the three dimensional layered system onto a spin 
Hamiltonian [335]. 

The numerical estimate of the exponent v in the multi-layer network model 
via finite size scaling analysis of the MacKinnon-Kramer scaling variable gives 
v = 1.45 ± 0.2 [100], which is again consistent with the result for three di- 
mensional unitary class. The quasi-energy spectral properties of the multi- 
layer network model [132,336] as well as wave packet dynamics [124] have 
been studied numerically, which are also consistent with that for three dimen- 
sional unitary class [130,170]. These results suggest that the bulk localization- 
delocalization properties of three dimensional layered Chalker-Coddington 
model (Chalker-Dohmen model) can be described within the conventional 
three dimensional unitary class. 



11.3 Chiral Metal 



Most of the discussions so far have assumed the periodic boundary condition in 
the transverse direction, Eq. (106). If we impose the fixed boundary condition, 
Eq. (107), the electrons travel along the boundary of the systems and form 
edge states. 

If we consider a multi-layer system where each layer has edge states circulating 
around the plane, the edge states are connected due to tunneling. When the 
Fermi energy is set between the center of the Landau band, bulk states are 
localized and the edge states are decoupled. In this situation, they form a 
sheath [100,337,338,339] in which the motion of electrons is directed (Fig. 54). 
The transport properties along the magnetic field (z-direction) is determined 
by this sheath, since the bulk states are all localized. We therefore can consider 
as independent the electron states in this two dimensional sheath which can 
be modeled by a two dimensional directed network model, Fig. 55. 

In this directed network model, at a saddle point, the incoming and outgoing 
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Fig. 54. The electrons circulate about the two dimensional plane in the case of fbc. 
They can hop from one layer to another, and form an "edge sheath" . 

waves are related via an S matrix as 





(407) 



In terms of the transfer matrix along ^-direction, the wave functions are related 
via 




-r'/t' l/t' 



{ l/t' 




(408) 



The resulting two dimensional network model is described in Fig. 56. It is 
similar to Fig. 8, but in this case the arrows are always rightward and S=S'. 

An interesting property of the transport in this edge sheath is that one can 
estimate the sheet conductivity along the z-axis exactly [100,339]. In the limit 
of large circumference C, the path connecting the bottom edge to the top can- 
not circulate along the system. This and the chiral nature of edge states lead 
to the fact that the interference between the paths with different wrapping is 
absent, i.e., a path does not self- intersect, and we can estimate the conductiv- 
ity classically. Let T = \t'\ 2 be the single layer transmission probability, and 
let T N the transmission probability for N layers (7\ = T). Then the following 
recurrence equation holds 

T N+1 = T N - 1 ——T , R= 1-T , R N = 1-T N (409) 

1 — rCriM 
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Fig. 55. Schematic view of electron motion in the two dimensional sheath. The 
arrows indicate the edge states circulating around the two dimensional plane. It 
tunnels to that in a different layer through interlayer coupling indicated by dotted 
lines. 



*x 



Fig. 56. The directed network model. It is similar to Fig. 8, but different in that 
the arrows are always towards the right in the x-direction, and only one type of the 
S matrix appear, not S and S' as before. 



or 



1 1 1-T 1 — T 1 
= 1 = N h — . 

T N+l T N T T T 



(410) 



This gives 



1 T 

Tn = N1-T(1-1/NY (411) 
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The conductance along the z-axis, G z , is given by 
e 2 e 2 N T 

G > = I N < T » = TTf i-rd-i/Af) (412) 



where iV c is the number of channels. Setting the distance of the saddle point 
along the x-direction to be b and the layer distance to be a, the expression for 
the conductivity a zz becomes 

L 2 T 

° zz = C° z = h°bl-T{l-l/N) ■ (413) 



where L is the length of the system along z-axis. For a sufficiently large system, 
one finds 

*" = -hbl^T- (414) 




Fig. 57. Numerical transfer matrix calculation of conductance G z (e 2 /h) as a func- 
tion of the layer transmission probability T. N c = 96 while N = 16, 32 and 96. Bars 
around the points indicate the conductance fluctuation. Solid lines correspond to 
Eq. (412). 



The results indicate than a zz can be much smaller than e 2 /h, but the system 
is still metallic. This is in contrast to the conventional two dimensional metal- 
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insulator [340] or two dimensional superconductor- insulator transitions [341], 
where the transition occurs near e 2 /h. This peculiar property is due to the 
chiral nature of the states in the surface, so that the system is called chiral 
metal. 

The conductivity is related to the localization length along the z-axis. By 
parameterizing the transmission eigenvalues, i.e., the eigenvalues of tv , r n via 
the Lyapunov exponent v n as (cf. Eq. (157)) 

(415) 



cosh^ (u n N/N c ) 
we have 



lim lim ~V 2 , Nc , (416) 

i-oo c-oo C h cosh 2 (u n N/N c ) y ' 



By noting that v n nv\ due to spectral rigidity, the expression of the con- 
ductivity becomes 



Na e 2 ^ N c . . 

a zz = lim lim — -— > ^— ; — - (417) 

n->oo n c ^oo N c b h ^ cosh 2 (z/ n iV/iV c ) v ; 

Nae 2 



hm hm — — — / cosh (uiNx) dx 

jV^ooAT c -+oc b h J 




e^a 1_ 

Since the smallest Lyapunov exponent is related to the one dimensional local- 
ization length £ 2 via 

N L ^ aN c . . 

^ = 6 ' (418) 



in terms of localization length the conductivity is expressed as 

o« = £§. (419) 



/From Eqs. (414) and (419) the MacKinnon-Kramer parameter is obtained 
[100] 

Ac= k = ^_I_, (420) 
C bl-T v ; 
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It is independent of the size of the system. This means that the wave function 
is critical. 

With the increase of the system length L, the paths of electrons from the bot- 
tom to the top begin to self-intersect. In this case, the transport properties are 
conjectured to be metallic. On the basis of mapping onto a ferromagnetic su- 
per spin chain [342], this conjecture was quantitatively discussed in [343,344], 
and then numerically verified [345]. 

In Fig. 58, the qualitative phase diagram of a layered quantum Hall system is 
shown. The width of the metallic region W is expected to increase with the 
tunneling amplitude like W(t) ~ ^Mqht^ where z/iqht is the critical exponent 
of the two dimensional quantum Hall transition. This behavior is obtained 
from the following argument for the derealization [100]: When the level spac- 
ing in a localized region A c ~ l/£ 2 with localization length £ is smaller than 
the tunneling matrix elements between the localized wave functions in adja- 
cent layers, the state becomes delocalized. This tunneling matrix elements is 
estimated to be ~ hence t ~ l/£ should hold at the mobility edge. Since 
in the two dimensional system the localization length diverges at E — like 
£ ~ \E\~ VIC * HT , one obtains the above dependence of the width W(t). 
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Fig. 58. Phase Diagram of the layered network model, where E = denotes the 
center of the lowest Landau band [100]. The width of the metallic region for small 
tunneling amplitude t is expected to increase like W(T) ~ tViQHT. 

The chiral metal has been theoretically predicted in 1995 [100,337]. Soon, it 
has been verified experimentally in an organic conductor [346] and in semicon- 
ductor quantum well structures [347,348,349,350,351]. Magneto-conductance 
[352], conductance fluctuations and the effect of interactions [353,354,355] have 
been discussed in relation to the experiments. 

Recently, layered network models have been considered for other symmetry 
classes as well. Since, as reviewed in the previous chapter 9, classes C and D 
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(disordered superconductor with broken time reversal symmetry, see Tab. 4) 
exhibit a rich phase diagram already in two dimensions, the corresponding 
phase diagram of the layered network can have additional phase boundaries 
to the one shown in Fig. 58 [356,357]. Still, the width of the metallic region is 
expected to be for small tunneling amplitudes t, W(t) ~ t 1 ^, where v is now 
the critical exponent of the corresponding quantum Hall transition at t — 0. 



12 Conclusion 

We have attempted to describe the development during the past decade of 
the random network model originally designed by Chalker and Coddington 
for the critical behavior of the quantum Hall phase transition. As the field is 
presently in a transient state with new ideas and developments appearing very 
rapidly, we cannot hope to have covered all of the different facettes completely. 
We can only hope that we have been successful in sketching at least the most 
important aspects such that a newcomer to the field can get an idea about 
what is going on. 

Two distinct and characteristic features of the model have been very important 
during the development. The first is of great practical importance. Similar 
to the tight binding Anderson model for localization, the network model is 
perfectly suited for numerical studies since the defining scattering operator is 
represented by a sparse matrix. Thus, quantitative numerical studies of the 
fundamental quantum critical properties of the model have been the subject 
of uncountably many works. These include not only the critical exponents 
but also the quantum fluctuations of the multi-fractal wave functions and the 
statistics of the conductances. 

The second property is perhaps of more fundamental nature. The network 
model can be mapped onto a great variety of Hamiltonians ranging from a 
bipartite tight binding Hamiltonian — of which the Dirac model is a limiting 
case, the Ising model to an antiferromagnetic chain of superspins. The versa- 
tility of the model invented by Chalker and Coddington for combining results 
from different areas seems to be truly unique. Using this, and the already men- 
tioned remarkable practical flexibility, reliable quantitative information about 
quantum phase transitions in very different kinds of disordered systems can 
be obtained which include all of the ten presently known universality classes 
of disordered quantum systems in two dimensions. 

Moreover, the mapping to the antiferromagnetic superspin chain and applying 
field theoretical methods has opened novel possibilities of putting the quantum 
Hall phase transition in a much wider context. Several new phenomena, such 
as the thermal Quantum Hall Effect, the spin Quantum Hall Effect, and the 
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chiral metal, have been predicted and are waiting for more theoretical and 
experimental efforts. Eventually, this also may contribute to explaining the 
universality of the quantum Hall phenomena which forms the underlying basis 
for the exactness of the quantization of the Hall conductance. 

Thus, the model can be considered as paradigmatic. It seems to us that the 
development has not yet come to an end. Many of the questions that have 
been raised during the development are still waiting for answers (Tabs. 3 to 5). 
These are especially the quantitative investigations of the critical properties 
at the boundaries of the various novel quantum phases that are predicted to 
occur in the models belonging to the different symmetry classes. 

As a major challenge, it remains to be explored how the model can be gen- 
eralized to include eventually interactions and correlation effects. There is 
evidence that electron-electron interaction is of great importance for under- 
standing the properties of the two dimensional electron system in the region 
of the integer Quantum Hall Effect [58,358,359]. Also, generalizations to the 
regime of the Fractional Quantum Hall Effect should be desirable [360] . There- 
fore, the generalization of the model towards including Coulomb interaction 
will eventually be crucial for getting insight into the physics behind the quan- 
tum critical phenomenon which seems to be of central importance in modern 
condensed matter physics. 
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